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Vill. SOLAR ROTATION AND PHOTOSPHERIC CIRCULATION 
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Summary 

From spectra exposed 1953 April 21-23, sight-line velocities have been 
measured at 332 points over the visible surface of the Sun. For any point the 
sight-line velocity can be expressed as a linear function of velocity com- 
ponents directed to the west and the north. Least-squares solutions of such 
observation equations for various latitudes reveal an asymmetry of rotational 
velocities in the equatorial zone and the presence of meridional currents. 
‘These motions suggest the existencé of a circulation in the photosphere 
similar to that shown in the upper atmosphere of the Earth. ‘This hypothesis 
is examined in the light of an investigation by M. J. Lighthill, and an obser- 
vational test is proposed. 





Introduction.—The rate of solar rotation is normally determined as an 
average in time and over the surface of the Sun. ‘This is necessarily so when 
successive meridian passages of sunspots are timed, but is equally true in the 
differential spectroscopic observations introduced by Adams (3) where the rate 
of rotation for a given latitude is found from half the difference of velocity between 
the northern and southern hemispheres. Provided photospheric motions are 
random their effect is minimized by such averages and the true rate of rotation is 
obtained. 

It is not, however, certain that photospheric motions are random. ‘Thus, 
when the systematic errors in Adams’ differential method are removed, a time 
variation in rotation is found (16) which may be due to long-lived local velocity 
fields, the existence of which has been established by A. B. Hart in Papers IV (8) 
and VI (g) of this series. These variations in velocity are not, as has been sug- 
gested, observational effects originating from the blending of sky and solar spectra. 
The suggestion arises from a misunderstanding of the origin and nature of 
scattered light, and is disproved by the simple experiment described in Paper LV. 

This paper describes an investigation of systematic photospheric motions. 
From observations made in 1953 April sight-line velocities, corrected for the 
observer’s motion, are found for 332 points on the solar surface, and are resolved 
at each point into east-west and south-north components. ‘The latter component 
is independent of rotation, and its amount makes possible an estimate of the 
contribution made to th. observed angular velocity by the east-west photospheric 
motion. ‘The observations and measurements are described in the first section 
of the paper which also contains a sample of the 332 sight-line velocities as functions 
of the heliographic coordinates. These velocities are broken down into their 
components in the second section, and the paper concludes with a discussion of 
the ‘ equatorial acceleration ’. 
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Summary 

From spectra exposed 1953 April 21-23, sight-line velocities have been 
measured at 332 points over the visible surface of the Sun. For any point the 
sight-line velocity can be expressed as a linear function of velocity com- 
ponents directed to the west and the north. Least-squares solutions of such 
observation equations for various latitudes reveal an asymmetry of rotational 
velocities in the equatorial zone and the presence of meridional currents. 
These motions suggest the existence of a circulation in the photosphere 
similar to that shown in the upper atmosphere of the Earth. This hypothesis 
is examined in the light of an investigation by M. J. Lighthill, and an obser- 
vational test is proposed. 





Introduction.—The rate of solar rotation is normally determined as an 
average in time and over the surface of the Sun. ‘This is necessarily so when 
successive meridian passages of sunspots are timed, but is equally true in the 
differential spectroscopic observations introduced by Adams (3) where the rate 


of rotation for a given latitude is found from half the difference of velocity between 
the northern and southern hemispheres. Provided photospheric motions are 
random their effect is minimized by such averages and the true rate of rotation is 
obtained. 

It is not, however, certain that photospheric motions are random. Thus, 
when the systematic errors in Adams’ differential method are removed, a time 
variation in rotation is found (16) which may be due to long-lived local velocity 
fields, the existence of which has been established by A. B. Hart in Papers IV (8) 
and VI (g) of this series. These variations in velocity are not, as has been sug- 
gested, observational effects originating from the blending of sky and solar spectra. 
The suggestion arises from a misunderstanding of the origin and nature of 
scattered light, and is disproved by the simple experiment described in Paper IV. 

This paper describes an investigation of systematic photospheric motions. 
From observations made in 1953 April sight-line velocities, corrected for the 
observer’s motion, are found for 332 points on the solar surface, and are resolved 
at each point into east-west and south-north components. The latter component 
is independent of rotation, and its amount makes possible an estimate of the 
contribution made to the observed angular velocity by the east-west photospheric 
motion. ‘The observations and measurements are described in the first section 
of the paper which also contains a sample of the 332 sight-line velocities as functions 
of the heliographic coordinates. These velocities are broken down into their 
components in the second section, and the paper concludes with a discussion of 
the ‘ equatorial acceleration ’. 
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1. Measured sight-line velocities 

In order to get as nearly simultaneous a survey as possible of motions over the 
whole visible surface a small image of the Sun is formed on the slit of the spectro- 
scope. The slit, oriented east and west, cuts the image in a chord. In Fig. 1 
let the x-axis be the slit with origin at O, marked by a wire crossing it, and let the 
y-axis be directed to the north. Then the equation to the limb of the solar image, 
centre C(x», y)) and radius pp, is given by 


(x ze xy)? + (y — Yo)" = Po 
The limb cuts the slit at A(x,,, 0) and B(x,, 0) so that 
Xo = }(x,, Uy Kes Po? — Vo" =} (s,,— =,). 

If now the south limb is tangent to a line parallel to the slit and distant a from it, 
then 

Yo =~ es 3 i(x,, x,)? 4a sae a) 

rfo= i [(x,, — %, y 4a + a). 

The coordinates wx,,, x, can be found from the distances measured along a spectrum 
line (corrected for curvature if necessary) from the shadow of the wire to the 


FIG. 1. 


bottom and top respectively of the spectrum, the distances being multiplied by a 
factor for magnification of the spectroscope. We can therefore find the polar 
coordinates p and ¢ (origin C, and ¢ measured counter-clockwise from the west) 
of any height on the spectrum, and consequently, as set out in Paper I (19), its 
heliographic and equatorial coordinates. 

In practice the distance a of the limb from the slit is not known. However 
two or more spectra, where the limb is displaced by known amounts from an 
arbitrary fixed origin, may be used to find p, and the y-coordinate of this arbitrary 
origin. We can then proceed as above to find the polar coordinates of any 
height in any spectrum. The advantages of the small image are thus not only 
that it permits a quick survey of motions over the Sun, but equally important that 
the data required for finding the heliographic coordinates of any height in a 
spectrum are given by the spectrum itself. The only quantities to be accurately 
determined when exposing the spectra are the U.T. of mid-exposure, and the 
displacement of the limb from a fixed but otherwise arbitrary zero. 
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1.1. Observations.—The prismatic spectroscope (18) of the 19-8m tower 
was used in the investigation. In order to have the strong telluric water-vapour 
lines as wavelength standards minimum deviation was set at 5900A, the dis- 
persion being approximately 1:5sAmm~!, The slit was covered with a diaphragm 
20 mm high crossed by three fine wires perpendicular to the direction of the slit. 
The central wire, where it crossed the slit (the x-axis), formed the origin O, and 
the two outer wires were used to find the magnification of the spectroscope. The 
slit was opened to 0:036mm. Ilford Thin-Film Half-Tone (R.50) plates were 
used because of their high photographic resolution. 

A specially corrected Wray objective of 5cm aperture and minimum focal 
length at sgooA of 148 cm formed a 13-7 mm diameter solar image on the spectro- 
scope slit. ‘The objective, mounted on an optical bench in the axis of the spectro- 
scope, was fed with parallel light from the Sun by the coelostat, secondary flat 
and Coudé mirrors of the 19-8 m telescope (17), the Cassegrain and the main 
mirrors of the telescope proper having been removed. ‘The objective was 
squared on and diaphragmed to 2-5 cm aperture to give a focal ratio comparable 
to that of the spectroscope. An Ilford filter, No. 202 in optical flats, before the 
objective admitted light on the long wavelength side of 5700A, while a mask 
before the filter and objective screened off unwanted solar light. 

The Guiding Plate, described in Paper I (19), was mounted in front of the 
spectroscope slit. The rocking plate was fixed normal to the incident beam from 
the objective, and setting was effected solely by the controls on the coelostat and 
secondary flat. A filar micrometer in the focal plane of the objective and eye- 
piece permitted measurement of distances and position angles. The coelostat 
and secondary flat were adjusted so as to make the east-west line in the sky parallel 
to the direction of the spectroscope slit, and a subsequent careful determination 
showed that this had been achieved within 0°-1. Correction for this and for a 
o°-3 lack of parallelism of the movable wire of the filar micrometer to the east- 
west line was made in the calculation of the polar coordinates of heights in the 
spectra. 


TABLE I 


Exposures 
No. Exp. time 
Series Spectra U.T. 1953 sec Remarks 
I April 21.541—21.621 8-9 No spots 

II April 22.378—22.461 75-11 Small spot on east limb on 
Ill April 22.576—22.624 9°5 22nd and 23rd. Three 
IV April 23.564—23.624 8-5 cloudless days after dis- 

persal of morning fog. 


With these arrangements four series of spectra, each series complete in itself, 
were taken as set outin Table I. As already noted each of the 14 spectra in any 
one series was a chordal section of the solar image. They were exposed in the 
following order: diametral section, 2mm S of diameter, 2mm N, 4mm N, 
4mm S, 6mm S, 6mm N, diametral section, 1mm N, 1mm S,3mm S,3mmN, 
5mm N,and5mm 8S. The whole image was thus covered in any one series by 
chords one mm apart. The displacement of the image on the slit by an exactly 
measured amount was secured by setting the south or north limb of the image 


16* 








200 H. H. Plaskett Vol. i19 


tangent to the movable wire, pre-set by the screw of the Guiding-Plate filar 
micrometer. 

After A. B. Hart’s investigation, described in Paper VI (9g), the spectroscope 
adjustments and temperature remaining unchanged, the main and Cassegrain 
mirrors of the 19-8 m telescope were restored. Spectra of the centre of the 18-2 cm 
diameter image, de-focussed, were obtained and one of these with an exposure 
of 22 seconds on an Ilford R.50 plate at U.'T. 1953 May 24.50g05 was subsequently 
measured for a standard. 

1.2. Measurement of spectra.—The spectra are measured on a Hilger photo- 
measuring micrometer with a mm pitch screw of 150mm length. Settings can 
be read directly to one micron. ‘The moving stage of the micrometer is provided 
with a subsidiary stage driven by a Moore & Wright micrometer (read directly to 
two microns) at right angles to the motion of the main stage. By measurement of 
a glass scale conversion factors are found and applied to reduce ‘mm’ measured 
with this micrometer, and also with the micrometer of the Guiding Plate, to mm 
of the long screw of the Hilger Micrometer. In the focal plane of the microscope 
of the Micrometer is placed a mask limiting the height of spectrum to be measured 
to o-1 mm; the opening in this mask is crossed by a single fine measuring wire. 

1.21. Standard spectrum.— he wavelengths of the lines to be measured on all 
spectra are listed in Table I]. For the water-vapour lines these are the wave- 
lengths recommended by the 1.A.U. (Transactions 3, 93, 1928; 4, 61, 1932). 
For the solar lines they are 1.A.U.—o-003A, as found from M. G, Adam’s inter- 
ferometric investigation (2) at 6100A and as partially confirmed by the measures 
at 5gooA quoted in Paper V (20). 


TaBLe I] 
Wavelengths 
‘Telluric water vapour Solar lines 


Wavelength R. Wavelength  Rowl. El Velocity 
I.A.U. No. M.G.A. No. = Factor 


'mm-! 


km sec 
5985°977 
5891°661 5892'880 : Ni 
5901°468 5905°677 Fe 
5919°054 E 5916°254 Fe 
5919°644 5927°794 Fe 
5924°272 5929°679 Fe 
5932°092 ‘ 5930°188 Fe 


w 
as 
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i as 
© os 


SIN sa“ sus 
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‘The settings for these lines on the standard spectrum of 1953 May 24 are 
measured direct and reverse on the Hilger Micrometer. ‘The measures are made 
at 12 heights within 1 mm of and symmetrically about the centre of the spectrum, 
corresponding to approximately the same distances east and west of the centre of 
the 18-2cm diameter solar image. ‘The setting for each line referred to an 
arbitrary origin, the same for all heights, and averaged for the 12 heights is then the 
standard setting to which measures of the small-image spectra are referred. 
From these standard settings and the wavelengths of Table II the sight-line 
velocity ly of the standard spectrum is +0°38kmsec~!. If 1.A.U. wavelengths 
had been adopted for the solar as well as the water-vapour lines this velocity 
would have been +0°23 kmsec~}. 
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1.22. Height measurement of small-image spectra.—Each of the four series of 
small-image spectra listed in Table I was measured at pre-computed heights in 
each spectrum and completely reduced before starting on a new series. The pre- 
computed heights corresponded to one or more of the heliographic latitudes 0°, 
+75, +1570, ..., +60°-0, +67°:5, and were found in the following way. 
Each spectrum of the series was ruled in a dividing engine across an unused spec- 
trum line with scratches separated at nominal 1 mm intervals from the shadow 
of the central wire (the origin O). These scratches were numbered and their 
distances, from the origin, and those of the top and bottom of the spectrum, were 
measured in mm of the long screw of the Hilger Micrometer. For this purpose a 
low-power objective was used in the microscope and an unmasked single mea- 
suring line in the focal plane. No correction for the curvature of the spectrum 
line was necessary. After multiplication by a factor for spectroscope magni- 
fication these distances gave the x-coordinates on the slit of the scratches and of 
the east and west limbs of the solar image. In accordance with the procedure 
outlined in the opening paragraphs of this section the polar coordinates and the 
corresponding heliographic latitudes of these scratches were found; from these 
the heights in the spectrum (relative to the scratches) of the required heliographic 
latitudes were obtained by graphical interpolation. 

1.23. Line measurement of small-image spectra.—With the Hilger Micrometer 
arranged as described in Section 1.2 to measure line positions over a spectrum 
height of a tenth of a mm, each spectrum of the series is measured direct and 
reverse at each of the pre-computed heights. For any one height the spectrum 
dispersion is reduced to that of the standard spectrum (Section 1.21) by means 
of the settings of the water vapour lines. ‘The ‘curve’ of errors is assumed linear 
and found by least squares. The difference between the reduced settings of the 
six solar lines and their settings in the standard spectrum, converted by the 
velocity factors of Table II to km sec! and averaged, gives us V, — Vo, the sight- 
line velocity of the particular latitude relative to that of the standard spectrum. 
For positions on the solar image where sin @= p/p, >0°775 the solar wavelengths 
of Table II must be increased by limb effect. From the interferometric deter- 
mination of limb wavelengths by M. G. Adam (2) the smoothed corrections of 
Table III, expressed in km sec~!, are derived and applied to the relative velocities 
V’,— J, of such positions on the image. 


Tasze III 
Correction to V,— Vo for limb effect 


P/Po A, P/Po A, 
km sec™! km sec”! 
0000 07000 6(«0°875) S—- — 0031 
0-775 07000 0°900 = — 0057 
0800 —0'002 0925 —Oo'r02 
0825 —0o°007 o'950 —0'182 
0850 —o'016 0975 —0°322 


The standard error s in the velocity of a single line is found as usual from the 
sum of the squares of the residuals (line velocity — mean velocity for six lines). 
The mean standard error for the 332 measured spectra is §=0-28(5)kmsec~. 
The individual values of s which contribute to this mean form a sampling distri- 
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bution which is satisfactorily represented by the theoretical sampling distribution 
for standard errors. ‘The true standard error o for a single line is thus given by 


pao, [2 I'(n/2) 
n—11'{(n—1)/2}’ 


or since our samples each have n=6, we find g=o0-30kmsec™!. ‘The standard 
error of the mean of six lines, that is of our relative velocity V,— V4, is thus 
0°30/,/6=0'17kmsec~!. The greater part of the error due to the measurement 
of the telluric water-vapour lines is carried forward into the solar lines in the 
objective adjustment by least squares of the dispersion of the measured spectrum 
to that of the standard. We therefore accept 0-17 kmsec™! as the standard error 
of V, — V4, and in view of the narrow spectral height of 0-1 mm available for fixing 
line positions may regard it as satisfactorily small. 

1.24. Final computations.—From the actual heights in each spectrum at 
which V, — V, is measured are computed afresh the polar coordinates of the cor- 
responding points on the image, and their heliographic coordinates B and L — Lp. 
From the polar coordinates of each point are also found the appropriate equatorial 
coordinates which are required for the correction C, for the Earth’s orbital and 
rotational motion. ‘This correction C, is found from eqn. (8) of Paper I (1g) in 
which we have taken the numerical constant for the orbital correction to be 1730-32, 
corresponding to a solar parallax of 8”-80, and have correctly found 6,’ from 
tan 5,’ =Z/V X?+ Y?. 

1.3. Results—The sight-line velocity V for an observer motionless with 
respect to the Sun is given by 


V =(V,— Vy) +0°38+ C, +A). 


In this equation +0:38kmsec~! (= V,) is the measured velocity of the standard 
spectrum (Section 1.21), C, the correction for the observer’s motion for the 
particular point on the Sun to which V refers, and A, the correction for limb 
effect (Sec. 1.23). 

The total number of values of V so obtained is 332, distributed in latitude as 
follows: 


No. 


~~ 
ee 
a 
° 


No. 


19 

; 13 

26 ° II 

21 - 7° 8 
16 


~ 
Nw Ss] 


Ww nN 
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This distribution is shown in an orthographic projection in Fig. 2 where open 
circles show the position of the Series I observations and filled circles the Series II 
observations. Since the Guiding Plate micrometer was centred with respect to 
the slit almost identically for Series II, III and IV the filled circles correspond 
to three determinations of V, the open circles to one. The position angle of the 
Sun’s axis of rotation at the time of observation was approximately — 26° (mea- 
sured from the north) and the slit ran from east to west. The successive positions 
of the slit on the image can therefore be seen by looking along Fig. 2 from east 
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south-east to west north-west when the observed points line up in a succession of 
slightly curved lines. 





4-2-0 + po 4 i 


ae es 


Fic. 2.—Orthographic projection of the Sun with lines of latitude and longitude spaced at 
15° intervals. Open circles indicate the positions where V has been measured in Series I, filled 
circles in Series IT. 

Representative of the complete* table of values of V is Table IV. The 
values of V are given as a function of the heliographic coordinates; an asterisk 
attached to V indicates that the correction for limb effect (‘Table III) is not zero. 
The roman numeral in the first column shows the series number of the obser- 
vations from which the velocity has been obtained. The standard error of V is 
o-17kmsec~! (Section 1.23), and the errors in the equatorial coordinates of the 
point on the Sun where V has been determined are not likely to exceed 

dbB=1°%2: 8(L—Ly)=2°:28ec B. 
These latter errors have been derived on the assumption that the error in fixing 
the x-coordinates of the solar limb where it crosses the slit does not exceed 0-2 mm. 

The sight-line velocity V at any point on the Sun can be resolved into com- 
ponents wu (directed to the west), v (directed along a meridian to the north) and w 
(directed radially outwards). We thus have 

U COS y; + VU COS 72+ W COS y3= V, (1) 
where the direction cosines of the sight-line with respect to these axés are from 
Paper I (19) 

cos y, = sin (L—L,) cos B, sin (@ + 8,) cosec 8 
cos y, = [sin B cos B, cos (L — Ly) — cos B sin By] sin (6 + 8,) cosec é 
cos y3= —cos (8+6,). (2) 


* The complete table will be issued as a supplement to Communications from the University 
Observatory, Oxford, No. 70. 
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These direction cosines are computed for each point of observation, and are 


shown for example in the representative Table IV. 


Tasie IV 


Sight-line velocities V for B= +37°°5 


~ 
Ww 


L-I, cosy cos V2 COs 3 V 
km sec™! 
+1972 +0329 +0640 —0'695 =+0°81 
or) be —O'519 +0°587 -—o'621 —0°59* 
+ 46° +O°719 +0°492 0492 +1°23" 
—0o'066 +0'677 = 0°732 —O'ol 
+ 0°093 +0°676 —o'°73!1 +o 
—0°836 +0402 —0375 —1'43* 
+ 0°483 +.0°599 —0°639 +0°25 
-0°343 +0640 —0o0'687 —O0O'S: 
to852 'org83. —0O357 +%° 
+ O°107 + 0°672 —o732 +0°9 
—0'S20 + O°412 — 0°397 —o'd 
511 +O°591 -0°624 +0 
"321 +0°645 —0'693 
39 +o'671 —0°728 
‘782 + O°445 0°437 
532 +0582 -0°616 


~J J <) 
Aa Au Ww 


~J J <J ~J « 


SWWW WwW WwW WwW Ww 
| i | 


297 +0647 —0°702 


2. Analysis of sight-line velocities 

Two points about the analysis of the measured velocities V are to be noted. 
First it is assumed that all points in the photosphere at any one latitude have a 
common motion. ‘This is almost certainly not true, but there are not sufficient 
observations to justify an analysis in longitude as well as latitude. Secondly any 
radially directed velocity component w, whether due to a general motion or, as 
originally suggested by St. John and recently revived by Schréter, due to local 
ascending and descending currents, will show itself as a change of measured solar 
wavelength from centre to limb. Since our wavelengths have been corrected for 
this limb effect by Table III, our measured V’s will contain no component 
velocity w, irrespective of whether such motion exists or not. Hence for the 
analysis of measured velocities eqn. (1) becomes 


UCOS y, + COS y,= V, (3) 


and this is the observation equation to be solved for the westerly u and northerly v 
components of velocity for each latitude. 

2.1 The westerly u-component.—As a preliminary to finding u and v for each 
latitude, a least-squares solution is made of eqn. (3) combining north and south 
hemispheres. ‘This involves reversing the sign of cosy, in the southern hemi- 
sphere so that the v-component in both hemispheres is directed to the nearer pole. 
The resulting mean values of 7 bring out nothing additional to the v’s for indi- 
vidual latitudes discussed in Section 2.2, but the mean values of u form a standard 
of reference and their smoothed values are plotted as a necessarily symmetrical 
curve in the upper half of Fig. 3. The ordinates in this figure are u— up, where 
Uy = 1°52 cos B is the linear velocity which a uniformly rotating Sun would have 
if its angular velocity were that actually found at |B|=67°-5. The curve shows 
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therefore the equatorial acceleration over and above that existing at |B] =67°-5. 
The smoothed value of u at the equator is 1-92kmsec~', while the numerically 
found value is 1-89. Both are within the range of equatorial velocities which have 
been found by other observers since 1915. 

The interest in the westerly u-component of velocity lies, however, in its 
values for individual latitudes. A least-squares solution is therefore made of the 
observation equation (3) at each latitude, north and south, for which |B|< 52°-5 ; 
individual latitudes poleward of this limit have with one exception less than 11 
determinations of V and are unlikely to yield reliable values of u and v. The 
results are given in Table V where n, in the second column, is the number of 
observation equations on which the solution is based, and where the values of u 
and v with their standard errors are in the last four columns. A point of minor 
interest attaches to the standard error o, of a single observation equation, given in 
the third column of the table. The mean of o, is 0-17kmsec~, the value pre- 
dicted in Section 1.23 from the internal consistency of the line measurements. 
This suggests that the local velocity-fields, described by A. B. Hart in Papers IV 
(8) and VI (9), and covering a linear distance greater (at least near the centre of 
the disk) than the o-1 mm height measured in the 13-7 mm diameter image, have 
not made any substantial contribution to the dispersion of the observed velocities. 


TaBLe V 
Velocity components as a function of latitude 
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The individual values of u with their standard errors from Table V are plotted 
as circles in the upper half of Fig. 3. The ordinates are u— up, that is the excess 
of equatorial acceleration at the latitude in question over that prevailing at 
|B|=67°-5. The interesting feature of the plot is the asymmetry of wu in the 
equatorial zones; for |B|<22°-5 the difference between southern and northern 
latitudes is approximately twice the standard error of the difference, that is on the 
5 per cent level of significance. Least-squares solutions have also been made for 
the individual latitudes for each of the four series separately, though the number 
of observation equations is dangerously low. Series I, Ii, and IV each show the 
asymmetry of Fig. 3, but Series III shows a reversed asymmetry—presumably 
the result of chance accumulation of errors in the small number of available 
observations. Asymmetry about the equator is not a new phenomenon in 
observations of solar rotation. Hulbrecht (10), from observations at Cambridge, 
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somewhat indirectly inferred that the northern hemisphere gave higher velocities 
than the southern during the first half of 1911 June. 
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Fic. 3.—Velocity components u and v as a function of latitude. 


2.2. The meridional v-component.—The v-component of velocity is signi- 
ficantly different from zero, and at low latitudes is larger than u—wu,, the excess 
equatorial acceleration. ‘| he individual values of v with their standard errors from 
Table V are plotted as circles in the lower half of Fig. 3, the ordinate scale being 
half that of u—wu, in the upper half of the figure. The value for B= —7°:5 is 
not shown since for this latitude the mean value of |cos y,|= 0-030 and the derived 
value for v is probably meaningless. ‘The figure shows well the character of the 
meridional currents which are directed in each hemisphere towards the nearer 
pole and which increase towards the equator. By and large these characteristics 
are repeated in the solutions for each of the four series, though the equatorial 
strengthening is more marked in Series land 1V thanin I] and III. Once again 
these differences between series, each with relatively few observations per latitude, 
are not to be regarded as significant. 

Since for a given latitude the coefficient of v, namely cos ys, is roughly constant 
(cf. ‘Table IV), any systematic error 6 in V will appear at that latitude as a 
‘meridional current’ of magnitude 5secy,. Two examples of this effect are 
given in Table VI, which is arranged in the same way as Table V. The upper 
half of Table VI gives the results of solutions of eqn. (3) for V — 0-15, that is the 
values of u and vw which would have been obtained if the I.A.U. wavelengths of 
solar lines had been used (Section 1.21). While the values of u are almost 
identical with those of Table V, as might have been anticipated, the meridional 
currents are now directed equatorwards. 
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Comparing Table V and the upper half of Table VI it might be supposed that 
with an appropriate choice of solar wavelengths the meridional currents would 
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disappear. To test this hypothesis a ‘systematic error’ 5 was found for each 
latitude from a least-squares solution of the equation 


ucosy,+d=V. 


The mean of these ‘systematic errors’, weighted according to the number of 
observations per latitude, was 5= +0-10okmsec™!. This is the ‘ systematic error’ 
which will most nearly make the meridional currents disappear. Solving now 
eqn. (3) for V—o-10 we get the values of u and v shown in the lower half of 
Table VI. It will be seen that though, in accordance with our intention, v is effect- 
ively zero everywhere, nevertheless at latitudes + 45°, + 30°, +15° and — 45° the 
values of v are approximately twice their standard errors. In short, no matter 
what we do, meridional currents persist, and it is only the latitude of their place 
of incidence and their direction which can be changed by applying arbitrary 
systematic corrections to the observed sight-line velocities V of Table IV. 

It would appear, therefore, that there are meridional currents in the photo- 
sphere with velocities comparable with those of the excess equatorial acceleration. 
Adopting the solar wavelengths of M. G. Adam, the system of meridional currents 
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is that shown in Table V and Fig. 3. That these currents have not been found 
before is probably the result of confining spectroscopic observations to the limb 
and the equatorial zone where, from eqns. (2), cos y, is small and where, as a con- 
sequence, the contribution of meridional currents to the observed sight-line 
velocities has likewise been small. 

3. The equatorial acceleration 

The existence of motions perpendicular to the direction of rotation suggests 
the presence of large-scale currents in the photosphere. ‘The fact that these 
meridional currents have speeds comparable with the excess equatorial acceler- 
ation raises the question whether the equatorial acceleration itself is not a mani- 
festation of photospheric circulation. If beneath the photosphere the Sun had 
a uniform rotation, then the equatorial acceleration and its asymmetry and the 
meridional currents would bear some resemblance to the general circulation of 
the Earth’s atmosphere. Here in the upper troposphere at all latitudes there is 
a mean westerly circulation, that is a greater angular velocity than that of the 
uniformly rotating Earth beneath. In a stimulating paper Lighthill (12) has 
suggested that the mean westerly circulation of the Earth’s upper atmosphere is 
due to a thermal wind, arising from the decrease in ground temperature from 
equator to poles, and deviated by Coriolis force to blow parallel to the isotherms 
and isobars. It is tempting to speculate whether a similar mechanism may not 
be operating in the Sun, and if so what the observational consequences would be. 

3.1. The hydrogen convection zone and polar cooling.—Assume, as is customary, 
the existence of Unsild’s hydrogen convection zone in which a not negligible 
fraction of the heat is transported by convective processes. Assume further that 
the top of this zone is a level surface rotating with a constant angular velocity w. 
Then if overstability is avoided, convection is inhibited by rotation and the 
critical Rayleigh number for its onset varies according to Chandrasekhar (5, 6, 7) 
as the two-thirds power of the Taylor number 


where @2 is the angular velocity at the place, v is the kinematic viscosity and h the 
depth of the convection zone. But Q=wsinB so that the Taylor number and 
the critical Rayleigh number steadily increase from the equator to the poles. 
Correspondingly we may expect an increasing deficiency of heat transfer by con- 
vection and therefore a cooling of the top of the hydrogen-convection zone as we 
proceed polewards. On these grounds we might expect that conditions on the 
Sun are such as to permit an application of Lighthill’s theorem. 

Application of Lighthill’s theorem.—Following Lighthill (12) let us 
gine the pressure gradient for an isothermal photosphere lying above the 
top of the hydrogen-convection zone. With our usual axes, x directed to the 
west, y to the north and z radially outwards, the pressure p at some level x’ =z —2, 
in the photosphere is 

p= poexp (—ng2'/RT»), 
where pf, and 7), are the pressure and temperature at z,, the top of the hydrogen- 
convection zone. From Section 3-1, 7, and possibly p, change with increasing 
y, and the pressure gradient at level 2’ i 


Cp _ dp, _ H&P oF ‘dT, ee 
ey ve (F RT,2 dy exp ( MSs RT,). 





No. 3, 1959 Solar rotation and photospheric circulation 209 


This pressure gradient attains a maximum at some height where its value is 


op Po aT, { ( TodPo/ >) Po aT 
= = * —°exp, — (1- =) pr 045 S.- 
3 Ty dy podT,/dy ’ To dy (4) 
On the right hand side we have assumed that the index of the exponential is unity, 
that is that the temperature 7, varies with y by a greater fraction of its absolute 
value than pp. 

According to Jeffreys (11) large-scale atmospheric motions fall into two main 
classes. First there are those in which the pressure and rotational terms in the 
equations of motion far exceed the inertial and frictional terms so that 

A ) 

—2wvsinB= — 12 
pox 
2wusin B= — +P ; (5) 

poy 
pcs J 
These may be called ‘ heliostrophic’ motions in analogy with the meteorological 
term ‘geostrophic’, and are such that the pressure gradient is balanced by 
Coriolis force. Secondly there are those in which the pressure and frictional 
terms exceed the inertial and rotational terms in which case the equations of 

motion take the form 
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u lo } 
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-e| 
The motion in this case is down the pressure gradient, viscosity, represented by 
the terms on the left hand side of the equations, preventing the velocity from 
increasing steadily. Such motions Jeffreys terms ‘antitriptic ’. 
The equatorial acceleration is thus to be regarded as a heliostrophic motion 
with dp/dx=o and op/dy as given by eqn. (4). We thus have from eqns. (5) 
Po AT, __ 4.4 RedTo 


2wu sin B= —o0'4 +> — = ’ 
+07, dy wp dy 


dy R Po Po 


at &, the top of the hydrogen-convection zone. ‘The numerical value holds for 
w=2'18 x 10-*sec~! (corresponding to a linear velocity of o-58kmsec~! at 
|B|=67°-5), for B=30° and for u=0-38kmsec~ (excess equatorial acceleration 
at this latitude), and for a mean molecular weight of unity. Since p/p)<1, the 
numerical value of 2-1 x 10-*degcm~", or equivalently some three degrees per 
degree of heliographic latitude, is an upper limit to the temperature gradient at 
the top of the hydrogen-convection zone which is required to produce the excess 
equatorial acceleration found at B=30°. It should be noted that the observed 


_ 41 =(55 cusin B) £ =(21 x 10-*) deg cm~ 
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usin B is not a constant, as would be the case if the temperature gradient and the 
density ratio were independent of latitude. 

In the equatorial zone where the Coriolis force becomes small we are passing 
from the heliostrophic regime of eqns. (5) to the antitriptic conditions of eqns. (6). 
It is therefore not surprising that we have meridional currents down the pressure 
gradient, that the excess equatorial acceleration changes slowly with latitude and 
that there are from time to time asymmetries in this acceleration. 

‘The temperature gradient of 3° per degree of latitude can only be regarded as 
a guide to the order of magnitude. A more exact calculation for a photosphere, 
not isothermal, but with a temperature gradient corresponding to radiative 
equilibrium would be justified if the limb-darkening test about to be mentioned 
were met. In any event such a calculation will demand highly accurate values of 
the total absorption coefficient for small changes of temperature and pressure, 
for in this problem we are concerned with the pressure gradient at a constant 
level, not a constant optical depth. 

3.3. Limb darkening.—lf Lighthill’s theorem is to be applicable to the Sun 
there must be a ditterence of temperature at the top of the hydrogen-convection 
zone between the pole and the equator, and correspondingly a difference of limb 
darkening along a polar and equatorial diameter. ‘To estimate its amount we 
assume the optical depth to the top of the hydrogen-convection zone always to be 
large, that is we make a distinction between this zone and the higher-lying zone 
of photosphere granulation (21). ‘Then for total radiation and a photosphere in 
radiative equilibrium on the Milne-Eddington approximation we have 


Polar diameter Equatorial diameter 
‘Temperature at 2, i i 
Optical depth at 2, 


‘0 


+ 


Temp. distribution T4=3F (7 


, ew ‘ ; © ns oa” 
Net flux F = i T3/ (7) +3) Ty'*/(t) +3) 
wr 
Surface brightness I(o, = B)=}F,(cos6+ 3%) = I'(0, @). = } Fy’ (cos @+ §) 
Hence a Fo = Toto +5 
Fy Ty ToT 3 
Since the temperature and pressure are higher at @ along an equatorial diameter 
than at 6=8B along a polar diameter, the equatorial absorption coefficient is 
higher than the polar coefficient, or ty)’> 7). That is 
I(o, @=B) _ T,* _ i T, —To 
rie, ~ I," 4 le 
If at the pole we take 7,’ — 7,=270° (corresponding to 3° per degree of latitude), 
and 7,’ = 10,000", then 
I(o, 6= B=7/2) 
ae (0, d=7 2) 
Hence, if Lighthill’s theorem is applicable to the Sun, there should be a 
measurable ditference of limb darkening between a polar and an equatorial dia- 
meter. With one exception the standard limb-darkening determinations have 
been made along an east-west or an approximately equatorial diameter. The 


> 0°89. 





No. 3, 1959 Solar rotation and photospheric circulation 211 


exception is that of Canavaggia and Chalonge (4) which, Professor Chalonge 
kindly informs me, was made along a diameter between east-west and north- 
south. This limb darkening is certainly different from that of the east-west 
determinations of Abbot (1), Raudenbusch (23) and of Pierce, McMath, 
Goldberg and Mohler (15), but is similar to the east-west determination of 
Moll, Burger and van der Bilt (14). Colour-temperature measurements at the 
limb are no more helpful. Thus Proisy and Gauzit (22) from an accurate 
visual determination find the pole to be 100° cooler than the equator, while photo- 
electrically Minnaert (13) finds no difference greater than his small error of 
measurement. Further measures, preferably of surface brightness, are clearly 
necessary, and on their result must stand or fall the applicability of Lighthill’s 
theorem to the photosphere. 

Conclusion.—The observations described in this paper show the presence 
of meridional currents and asymmetry in equatorial rotation. If, as in Section 
3, these motions are regarded as a purely photospheric circulation it follows, 
since sunspots also show equatorial acceleration, that the visible part of a spot 
is not anchored to the top of the hydrogen-convection zone, assumed on this 
interpretation to be in uniform rotation. ‘The motions, however, may be inter- 
preted quite otherwise. For example the spectroscopic abundance of beryllium, 
deuterium and lithium shows that there is a material connection between the 
photosphere and the solar interior ; it is, therefore, not impossible that the photo- 
sphere is revealing in a modified way motions which have their origin in the 
deep interior. Whatever the interpretation, however, the immediate task is by 
other and improved methods to confirm or disprove the photospheric motions 
found in the present investigation. 

Work in this and related fields is planned at Oxford. Since the methods will 
ditfer from those hitherto employed the present paper will conclude the exploratory 
series Motions in the Sun at the Photospheric Level. For invaluable discussions 
and in particular for their own contributions to the series I wish to express my 
thanks to Dr Hart and to Dr Kinman. 


University Observatory 
Oxford: 


1959 February 27. 
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Summary 


Observations were made during spring 1958 of the central emission 
intensity of the Sun at a wave-length of 11-104 (g01 cm™'), by direct 
comparison with the emission of a black body at 1300°K. ‘The solar 
intensity was reduced to approximately that of the black body by the use 
of rotating sector disks, and the ratio of intensities was then measured by 
a calibrated electrical attenuator so that the system was independent of 
the detector linearity. The only subsidiary measurements needed were the 
absolute reflectivity of one mirror and the variation of reflectivity with angle 
of reflection of another. 

The slit width of 0-031 or 2°5 cm~' was small enough to enter an 
atmospheric ‘“‘ window "’ where the terrestrial absorption obeys Lambert’s 
law, thereby allowing extrapolation to extra-terrestrial intensity. 

The weighted mean of five measurements of the central intensity was 
IQ (11°10) = 2°408 + 0°016 x 10" erg cm~ ster~! sec~! (per cm wave-length 
range). This is equivalent to a central temperature 7G (11°10 4) = 5036 + 30°K. 





1. Introduction.—'The astrophysical importance of an accurate determination 
of the energy—wave-length curve of the Sun is well known (1). In recent years 
attention has been directed towards the longer wave-lengths, ranging from the 
near infra-red (2, 3) to the radio region (4). 

No measurements of central radiative intensity seem to have been made in 
the atmospheric window around 1oyz. This is rather surprising, as both solar 
and terrestrial conditions are here somewhat favourable; there are virtually 
no Fraunhofer lines, so that a true photospheric continuum may be observed, 
and the Earth’s atmosphere is relatively transparent between the intense bands 
of H,O at 6-3 and of CO, at 15. Optical techniques in this region are now 
well developed, and—a point of great practical convenience—the energy ratio 
between the Sun and any possible terrestrial comparison source is very much 
less than at shorter wave-lengths. For example, Peyturaux (2), using a tungsten 
ribbon lamp at 2900 °K, had to measure an energy ratio of 20 at 2:3; in the 
present determination the reference source was at 1300 °K and the ratio was 6. 
At this lower temperature an almost-perfect black-body source may be constructed 
fairly easily, enabling a direct comparison with the Sun to be made with a very 
simple optical system. 

This paper describes a series of measurements made during spring, 1958, 
at the Ascot (Berkshire) Field Station of the Department of Meteorology, 
Imperial College. A grating spectrometer was used, enabling the observations 
to be made in a small ‘‘ window”’ at 11-10 (go1 cm~") (5, 6). Although there 
are no known molecular lines in this region, the Earth’s atmosphere is by no 
means transparent, and the study of its opacity forms another aspect of the 
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observations which will be reported elsewhere. It has been studied at lower 
resolution (7) showing that there is a continuous absorption obeying Lambert's 
law over most of the region in the spectrum from 8 to 13, so that extrapolation 
to zero air mass can be made with confidence. 


2. Apparatus. (a) General layout.—The equipment is arranged as shown 
in Fig. 1. The siderostat mirror M was electrically driven and maintained a 
beam of solar radiation on the paraboloid M, (focal length 180cm). When 
M was not pointed at the Sun it was fixed so as to reflect a collimated beam 
of radiation from the black body by M, on to M,. The paraboloid mirror M, 
then forms an image of the Sun or of the black body on the entrance slit of the 
spectrometer, placed just below the incoming beam. It will be seen that 
radiation from the Sun and from the black body follow identical paths except 
that (i) the black-body radiation suffers one additional reflection at the collimating 
mirror M,, and (ii) the angles of reflection at the siderostat mirror differ. It 
therefore becomes necessary to determine the absolute reflectivity of M,, and 
the variation of the reflectivity of M with angle of reflection; the methods used 
are described in Section 3. 

The aperture of the system is defined by a metal diaphragm placed at 5, 
the calculated position of the entrance pupil. 


780cm 
is 


8.8 
Fic. 1.—Layout of the observing equipment. The siderostat mirror M is shown in position for 
taking measurements on the black body B.B. The sector disks are placed at S.D in front of the 
entrance slit of the spectrometer. The dotted plane mirror M, (see Section 3(b)) was placed in this 
position only when the variation of the reflectivity of M was measured. 


(b) The black body.—The black body is shown in Fig. 2. It consists of a 
sphere of sillimanite incorporating a heater winding of Kanthal Al wire 
(m.p. 1520 K) which can be run directly from a 240 V a.c. line. An additional 
winding on the exit tube maintains the part adjoining the sphere close to the 
sphere temperature and so reduces the disturbance caused by the circular exit 
hole. ‘The inside wall of the sphere opposite the exit tube was roughened to 
prevent any possibility of specular reflection and the whole of the inside was 
coated with cobalt oxide (applied as a paste mixed with dextrine and water). 
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A careful exploration was made of the internal temperature distribution by 
means of silica-encased platinum, platinum—-rhodium thermojunctions passing 
vertically and horizontally through the sphere, from which it was established 
that the temperature was uniform within +1 °C. The emissivity of cobalt oxide 
at 10 and room temperature is 0-8 (8). Assuming this value to be true at 1300 °K, 
the normal operating temperature, it has been calculated that with the given 
geometry the departure of the emergent beam from perfect black-body radiation 
is less than o-1 percent. This is about the same error as is introduced in measuring 
the temperature of the black body. 





Ad 





Fic. 2.—Black body. The heated sphere (A) is supported by three refractory bricks (B). The 
small tube (C) inside the exit tube (ID) is separately heated to the same temperature as the sphere. 
The box is filled with “ filter cell’? powder (F). The thermocouple slides inside the projecting 
sheath of quartz (E). 


The black body was run continuously from a constant-voltage transformer, 
a temperature-controlled cut-out being set at a somewhat higher temperature 
than the normal equilibrium temperature as a safety precaution. The temperature 
was measured using Pt-Pt 13 per cent Rd thermocouples calibrated by the 
National Physical Laboratory to an accuracy of +2°C at 1300°K. The 
reference thermocouple was placed in a Dewar vessel of water whose temperature 
was measured with a mercury thermometer, while the other thermocouple was 
placed inside the heated sphere in a selected position. A Tinsley potentiometer 
in conjunction with a Cambridge Spot Galvanometer and a Weston Standard 
Cell was used to measure the e.m.f. The system was capable of detecting 2V, 
while the total e.m.f. produced was 11004V. The temperature varied by less 
than 1 “C during any run. 


(c) The spectrometer.—The infra-red spectrometer was a small double 
monochromator using two 7x5cm Merton N.P.L. replica gratings having 
1000 lines per cm and blazed for 10 in the first-order spectrum. The elimination 
of the higher order spectra of shorter wave-lengths was achieved by the use of an 
indium antimonide filter with an anti-reflection coating of lead chloride (9). 
This remarkable filter has an extremely sharp cut-off at 7, with maximum 
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transmission of 98-5 per cent at 10-2, falling to go per cent at 8-2 and 13-5. 
‘The detector was a Golay cell (Unicam) with a phase-sensitive rectifier and a 
Honeywell-Brown potentiometer recorder. ‘The radiation falling on the Golay 
cell was chopped at 11 c/s. Fig. 3 shows the resolution obtained in the solar 
spectrum and illustrates the efficacy of the combination of gratings with 
semiconducting filters in a small and relatively inexpensive instrument. 

‘Two careful tests for stray light from shorter wave-lengths showed this to 
be an entirely negligible factor. One test was to intercept the solar radiation 
with metal and glass shutters alternately, the glass shutter admitting radiation 
below 2:5 only. ‘The other test was to compare the solar radiation directly 
with that from a metal shutter in a region of very strong atmospheric absorption. 
In both tests the scattered radiation was found to be less than 0-05 per cent. 

It is advantageous in photometry to employ a null method, making the 
brighter beam equal to the weaker by some form of calibrated attenuator, so 
that the measurements are independent of the properties of the detector. In 
practice, however, a continuously variable and accurately calibrated attenuator 
is notoriously difficult to realize (10). In the present work a compromise has 
been made. ‘The intensity of the solar beam was brought into approximate 
equality (+1 db) with that of the black body by the use of one of a set of four 
rotating sector disks with transmissions of 0°6227, 0°3991, 0°2504 and 0-1568. 
These sector disks were accurately machined and calibrated by geometric 
measurement. ‘They were placed in front of the chopper and interrupted the 
beam at a frequency of 100 c/s. Since the time constant of the Golay cell is of 
the order of o-1 sec, the sectors should produce the theoretical attenuation 
independently of the 11 ¢s chopper; this has been confirmed experimentally 
by showing that there was no detectable variation in output with sector speeds 
between 50 and 100c s. The measurement of the ratio of intensities of the two 
nearly-equal beams now places far less stringent demands on the detector linearity. 
As a further precaution a calibrated electrical attenuator, which can be switched 
into the circuit to replace the normal gain control of the amplifier, was used 
to eliminate non-linearity in subsequent stages. 

A careful test of the detector linearity was carried out using a constant signal 
from the black body and the fast rotating choppers as a means of changing 
the applied signal in a known way. This tested the linearity of the Golay cell 
and the amplitier combined. The linearity of the amplifier alone was found 
using the electrical attenuator as a means of changing the applied signal. From 
these tests it was concluded that errors in the measurements due to lack of 
linearity were entirely negligible. 


3. Reflectivity measurements. (a) Collimator mirror for black body (M,). 
he absolute reflectivity of mirror M, was measured by setting it up on a special 
mount in the laboratory so as to form an image of a large vertical Nernst filament F 
at I (Fig. 4). Images of F and I were formed in the plane of the spectrometer 
entrance slit by M, and M,, and by adjustment of the latter either of these may 
be made to fall on the slit. A good quality flat was used for M, to avoid errors 
arising from the use of slightly different parts of its surface for the two images. 
An aperture stop S was placed at the image of M, by M,. Possible non-uniformity 
of temperature over the Nernst filament was allowed for by repeating the 
measurements after rotating it through 180 . 
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‘The holder of M, was attached to its stand by a kinematic device, enabling 
the reflectivity measurement to be made at intervals with minimum disturbances 
to the observing system. 

The reflectivity of this mirror appears to have remained effectively constant 
over the observing period, as shown by the results in ‘Table I. 
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Fic, 4.—System used for measurement of reflectivity of collimator mirror M, 


TABLE | 


Reflectivity of collimator mirror M, 
Date of aluminizing: 1958 February 20 
Date Reflectivity 

1958 


March 21 0°979 + 0°004 
May 7 0°979 +0°005 
° 
° 


June 23 ‘976 +0°004 
June 30 ‘983 +0°005 
\dopted mean: 0°979 +0°003 


(6) Stderostat mirror M.-—Only the variation of reflectivity with angle of 
reflection is required. For this measurement a subsidiary plane mirror M, 
was mounted on an arm which could be rotated about the vertical axis of the 
siderostat (see Fig. 1) he siderostat mirror was also constrained to rotate 
about this axis and the azimuth angle could be read from a divided circle 

\ Nernst filament placed close to the spectrometer slit sent light to M, after 
retlection from M, and MM. M, was adjusted to return this light along an almost 
dentical path, so that an image of the filament was on the slit Phe sensitivity 
of the method was improved by the double reflection at M Lhe angles of 
reflection at the mirrors other than M were kept constant L ntortunately 


he reflecting surface of \I wa - om off axis. so that lateral motion of the image 
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of M, occurred for different reflecting angles. A small stop was therefore placed 
at the image of M, by M, to ensure that the same reflecting area of M was used 
for any angular setting (this stop, inside the spectrometer, is not shown in Fig. 1). 

Results are shown in Fig. 5 from which a table was drawn up for correcting 
observational results to a constant angle of reflection, namely that of the black 
body beam: 37°:5. ‘The observed decrease of reflectivity with increasing angle 
of reflection may be expected from the presence of a thin layer of dust, which will 
obscure an increasing fraction of the aluminium surface as grazing incidence is 
approached. 


Reflectivity in 
ost Arbitrary units 


06 + 





95r 





94 aro ae j j 
0 10 20 30 40 50 60 


Reflection Angle in degrees 





Fic. 5.—Variation of reflectivity of siderostat mirror M with angle of reflection. 
Dates of measurements : 1958 May 24, © to East, @ to West. 
1958 Fuly 20, (_) to East, §§ to West. 
No change appears to have occurred over this period and the curve shown 
zvas used to correct all the observations. 


4. Observations.—The observations reported here were made in the window 
at 11:10u (gor cm~'); the slit was 7mm high and 1-9 mm wide, the latter 
corresponding to a spectral slit width of 0-031 (2:5 cm~*). The diameter of 
the solar image was 17mm. Further series of observations in windows at 
12°02 4 (832 cm~') and 8-63 1 (1159 cm~*) are planned when weather conditions 
permit. 

Observations were made from sunrise until deep convection started at 8" 
or g" G.M.T.; this covered an adequate range of air masses for an accurate 
extrapolation to zero air mass. Alternate recordings were taken on the Sun 
and on the black body, using the sector disks and electrical attenuator to make 
the two deflections as equal as possible. The balancing of the two beams was 
made to within 2 per cent. When this condition was obtained, the electrical 
attenuator readings and the ratio of the recorder deflection for the Sun and the 
black body, together with the reflectivity measurements discussed above, the 
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transmission of the sector disk used and the temperature of the black body 
enable the mean solar intensity across the slit (1) to be found for each air mass. 

A small correction, at most 0-6 per cent of energy, was applied to allow for the 
finite emission of the chopper, which was at room temperature. The gain of 
the electronics changes smoothly by approximately 10 per cent during a four 
hour run, but the change between any two observations of the black body did not 
exceed 0-3 per cent and a linear interpolation was used to obtain the black-body 
deflection at the time of the Sun observation. 

The air-mass values calculated from the time were corrected for curvature 
of the Earth’s surface and atmospheric refraction ; for this purpose it was assumed 
from the work of Roach (11) that the atmospheric absorption in this part of the 
spectrum can be assigned to a thin layer at a height of 1-5 km (850 mb). The 
refractive index of air at 11-1 was taken to be 1-0; 273. 

Log (J) was then plotted against air mass and the best straight line was 
obtained by the method of least squares. This line was then used to find log (J,) 
by extrapolation to zero air mass. A correction for limb darkening must be 
applied to give the central intensity />. At this wave-length, and for the portion 
of the solar disk sampled, the correction is very small, and from limb darkening 
curves of Pierce, McMath, Goldberg and Mohler (12) at 10 and results of our 
own at 11°10 we find a correcting factor of 1-002. 


5. Results.—The weather in the spring and summer of 1958 was so poor 
that a daily watch at dawn yielded only five useful runs. The plots for individual 
runs are given in Fig. 6 and the resulting intensities and temperatures in Table II. 


9.9 


‘ Corrected air masses 
9. 1 mn oe L i i aa 
0 2 4 6 8 10 12 4 


Fic. 6.—Results of observations : @ 1958 April 2. O 1958 April 12, [_] 1958 April 30, 
@ 1958 May 1, x 1958 Fune 18, 
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Tase II 


Results of observations of solar emission at 11°10 
I, 
Date S logie 7) w k erg cm~? ster~! sec-! o 
(1958) percent (J/B), °K gmcm-* (air mass)! (per cm wave-length per cent 
range) 

April 2 o's O°755 1321 0°95 0°0344 2°408 x 10" o’81 
April 12 or) 0°757 1308 0°55 0°0471 2°384 x 10% 0-74 
April 30 o°3 0760 1312 18 0°169 2°411 x 10 1°46 
May 1 1° 0°764 1313 1°s5 0°138 2°437 x 10° 1°83 
June 18 ae) 0°773 (1307 1°6 0°179 2°472 x 10° 1°63 


S = Maximum fraction of slit area covered by sunspots. 

k= Extinction coefficient of atmosphere, base e. 
T, = Temperature of black body. 

B= Intensity of black body at 7). 

I, = Average solar intensity over the slit outside the Earth’s atmosphere. 

o= Standard deviation of Jy. 

w= Atmospheric water-vapour content in the zenith. 


An evaluation of the standard deviation of the regression lines enabled the 
standard deviation of log J, to be found. Hence the standard deviation of J, 
was deduced and is included in Table II. 

The values of the sunspot area (S) in Table II have been estimated from 
photographs kindly supplied by the Astronomer Royal. Their effect on the 
measurements is probably negligible, since the brightness contrast between spot 
and photosphere is likely to be small at this wave-length. The atmospheric 
extinction coefficient (k), deduced from the slope of the lines in Fig. 6, is also 
included in Table II, and so is an estimate from Meteorological Office Daily 
Aerological Records of the water-vapour content of the atmosphere, with which 
it shows a rough, not necessarily direct, correlation. The origin of this extinction 
is of considerable meteorological interest and is being further investigated. 

Taking a mean of the intensities in Table II with weights proportional to 
the inverse square of the standard deviation, and applying corrections for limb 
darkening and deviation of the black body from a perfect enclosure, we obtain 
I¢(11-10p) = 2-408 x 10” ergcm~*ster~! sec! (per cm wave-length range). 


6. Discussion.—The standard deviation of the weighted mean of the five 
measurements of Jo is 052 per cent. This includes indirectly the effect of 
errors in the measurements of the records, and any rapid changes of the 
atmospheric transmission. There are also possible systematic errors of 
measurement as follows: reflectivity of M (0-3 per cent); temperature of the 
black body (0-2 per cent); variation of reflectivity of M (o-2 per cent); 
attenuation of the sector disks (0-04 per cent). Combining these, an estimate 
of the error of 1g was found to be +0°67 per cent. Hence, 
1¢ (11-102) = 2-408 + 0-016 x 10” ergem~*ster—!sec~!(percm wave-length range). 
Using the data of Allen (14) this intensity corresponds to a black-body 


temperature of 


Te(11'10 p) = 5036 + 30 °K. 
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This result taken together with those of Peytureux (2) at 1°55, 2°03, 2°31 
where To was found to be 6330, 6045 and 5880°K respectively (the results 
of Pierce (3) are in excellent agreement) is in line with the expected decrease 
of temperature with increasing wave-length. It also fits in with an effective 
temperature of 4500 °K at 6 mm as measured by Whitehurst et al. (4). There 
is however a possible discrepancy with the measurements of Coates (13) who 
gives a temperature of 6962 °K at 4-5 mm. 
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his interest and advice in the astrophysical questions involved; and to 
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THE MAGNETIC AND DYNAMICAL FIELDS 
OUTSIDE A PROTO-STAR 


L. Mestel 
(Received 1959 March 7)* 


Summary 


A contracting proto-star formed in a magnetic cloud conserves its magnetic 
flux provided the plasma density does not become too low. In this paper the 
structure of the external field is studied, together with its effect on the flow of 
matter. Centrifugal force is assumed always negligible, and the star is 
supposed at rest in the cloud. The fields in the star and at infinity are both 
taken as uniform, and the flow is assumed isothermal. 

Possible equilibrium fields include : 

(i) A poloidal field, with all the lines of force extending to infinity, but 
with a small toroidal component strong enough to keep the centrifugal field 
small in spite of the inflow of matter. The field is weak at the equator, but 
the lines of force make sharp hairpin bends there, yielding a pinched 
equatorial discharge. 

(ii) A poloidal field, with some of the stellar field lines finite, separated 
by a ring of X-type neutral points from the infinite field lines, which have a 
structure similar to field (i). 

(iii) A twisted, torque-free, and approximately force-free field confined 
within a finite surface S, and an approximately curl-free field beyond S. 

If the cloud field initially contains toroidal flux threading poloidal loops, 
then type (iii) will result; otherwise, types (i) and (ii) are the relevant 
equilibrium fields, applying respectively to earlier and later stages in the star’s 
contraction. 

The presence of the field alters the accretion problem somewhat. Solutions 
exist with supersonic and subsonic velocities at infinity, yielding respectively 
an increased or diminished accretion rate. In the absence of an external 
strong pressure field at infinity, the system is likely to approach the state with 
the lower rate, so that the presence of the field reduces the accretion rate 
substantially below Bondi’s value. 





Introduction.—The problem of star formation in a strongly magnetic gas 
cloud was discussed in a recent paper (1). If the cloud field is assumed to have 
the strength required by current theories of interstellar polarization, cosmic-ray 
acceleration, and synchrotron radiation, the cloud must have a mass at least of 
order 10°© if it is to withstand the magnetic and thermal pressures. If, further, 
the field is frozen into the cloud, this lower limit is hardly reduced if the cloud 
contracts, and it is difficult to see how masses of stellar order could spontaneously 
form. The difficulty is resolved when it is recognized that the magnetic field 
in an Hi cloud is tied to the plasma rather than to the gas as a whole, so that 
if the frictional coupling between plasma and neutral gas is weak enough, the 
bulk of the cloud can contract without dragging the field with it. In a dense 
dust cloud, obscured from galactic ionizing radiation, the plasma density can 
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decay rapidly to a value low enough—approximately one per cent of the 
equilibrium value—for the two gases to be effectively uncoupled; neutral matter 
condenses across the magnetic field, high densities are built up without a 
corresponding high magnetic energy, and stars can form without magnetic 
impedance. 

We now suppose that a proto-star has been formed, and that effective 
freezing-in of the magnetic field into the stellar matter has been re-established. 
This is certainly the case when the internal temperatures exceed 10‘°K, the 
ionization temperature of hydrogen, while at lower temperatures the plasma— 
hydrogen friction for normal plasma densities is enough to prevent any sensible 
drift of plasma plus field through the neutral gas (1). In this paper attention 
is focused on the magnetic field external to the star, and its interaction with the 
flow of matter due to the star’s gravitational field. ‘Thus we are interested to 
know whether an initially poloidal field remains so, or whether it acquires a 
toroidal component; at what stage the stellar lines of force break off from the 
galactic field; and whether the field affects substantially the rate of accretion of 
matter by the star. 

We consider the case of a contracting proto-star-at rest in a cloud. Although 
we expect the cloud to be lightly ionized (an Ht region), we shall also discuss 
the simpler case of a fully ionized (H11) region. There is a third possibility: 
the cloud may be lightly ionized beyond a critical surface, but wholly ionized 
within it. ‘This could be due to compressional or frictional heating of inflowing 
neutral hydrogen, or, if the star is massive and has reached the main sequence, 
to stellar ultra-violet radiation. ‘This case, which will demand study of the 
hydromagnetic shock at the critical radius, is postponed for separate treatment. 

The magnetic field at large distances from the star (H,) is supposed uniform, 
it being a part of the local galactic field, and maintained by currents that flow in 
dimensions large compared with the relevant part of the cloud. For simplicity, 
the stellar contraction is taken as approximately homologous, so that within the 
star the field (H1) is uniform and parallel to Hy, except in a narrow surface region 
(cf. Section 2). The name “stellar flux tube ’’ will be given to the tube generated 
by those lines of force that cut the star. 

In the absence of a magnetic field, the velocity—density field of all the cloud 
matter (neutral and ionized) would approach the spherically-symmetric, 
temperature-limited pattern discussed by Bondi (2). For simplicity, let us 
suppose that initially the magnetic forces are weak, so that the field does not 
impede the flow, but rather the ionized matter distorts the field until it begins 
to exert forces comparable with gravity. Clearly a steady state is approached, 
in which no further dragging of the field takes place, the flow of ionized matter 
across the field being just the rate permitted by the high conductivity. The 
magnetic force is then balanced by the components across the field of the 
non-magnetic forces, including inertial forces, acting on the plasma——a familiar 
situation in steady state hydromagnetics. In particular, we shall find that in 
some regions the magnetic force compresses the plasma so as to build up a pressure 
gradient strong enough to balance it, just as in the well-known laboratory “‘ pinch 
effect ’’. 

Simultaneously, the flow pattern approaches a steady state. Except in regions 
of very low field strength, the flow of plasma across the field is negligible; but 
along field lines of infinite length, there will be a steady flow of plasma, either 
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directly into the star, or towards the equator outside the star. The plasma 
density at the equator will grow until the rate of influx down the field lines is 
balanced by the rate of recombination, so that the neutralized plasma can drift 
across the field into the star*. If the cloud were fully ionized all the flow would 
be of this type; however, in a lightly ionized cloud neutral matter can drift 
across the field, but its motion is restricted by the frictional drag of the plasma. 
Although neutral matter does flow across the field into the star, Bondi’s 
spherically symmetric velocity field will be distorted even for the neutral gas 
by the indirect effect of the magnetic field, acting through plasma—hydrogen 
collisions. 

We are assuming implicitly that rotation in the cloud is never rapid enough 
to yield a centrifugal field comparable with gravity. If this is so at any time, 
then a strong enough magnetic field will ensure that the inflow of matter does 
not lead to the subsequent build-up of a strong centrifugal field. However, 
if the star is rotating it will tend to set the cloud rotating as well; we require that 
this transfer of angular momentum from the star to the cloud is rapid compared 
with the rate of contraction of the star, so that the stellar angular momentum 
is quickly distributed over the cloud, and no important centrifugal field arises. 
This at once shows the importance of knowing in detail the structure of the stellar 
magnetic field: clearly, if many lines of force have broken off from the galactic 
field and so donot extend toinfinity, the field will be much less efficient in removing 
angular momentum. 

Although the discussion of this paper is in terms of a proto-star containing 
a frozen-in magnetic field, similar conclusions could apply on a larger scale. 
‘Thus consider a freely-falling gas cloud of say 10° solar masses, and with a strong 
frozen-in magnetic field: the rather special conditions discussed in (1) for 
uncoupling of the field and the bulk of the cloud may not obtain, so that the 
frozen-in field prevents break-up into smaller masses. Then if the scale of the 
local galactic field is large compared with the cloud’s dimensions, the problem 
is essentially that of this paper. As long as most of the cloud field lines are of 
infinite length, magnetic transport of angular momentum will probably keep 
low the rotation of the cloud, and no strong centrifugal field, either internal or 
external, will be built up. However, once the cloud field is mainly broken off 
from the galactic field, again we may expect much less efficient removal of angular 
momentum, and the probable build-up of a strong centrifugal field within the 
cloud; this will slow up or halt the cloud’s contraction, and in particular allow 
more time for the strong internal magnetic field to leak out, as discussed in (1). 
In this way excess magnetic energy may be lost, at the cost of leaving a rapidly 
rotating cloud with a weak magnetic field. Further discussion is beyond the 
scope of this paper, but we may anticipate that the further evolution of the cloud 
would be in part determined by the slow interaction of the magnetic field with 
the rotation field (3). 

Part I 
The magnetic field 

1. The equations to the field outside the star.—Although the boundary conditions 
imposed on the field refer only to the poloidal components, we allow for the 
possibility that the field has both a poloidal component H,, and a toroidal 


* The probable importance of recombinations in this context was pointed out to the author by 
Dr A. Schliiter. 
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component H,, both assumed symmetric about an axis. The magnetic force 
density can then be written 
(VaH)sH — ((VaH,)sH,+(VaH,) aH) (poloidal component) 
47 47 
(VaH,aH 
47 
the term (V,H,),H,/47 vanishing identically. 
The poloidal part of the magnetic force has a component 





(toroidal component), 


tod 
-| (VaH)aH,-H,| ai. 
47H, 


in the direction of H,, and a component 


Pp?’ 
[ vats, (ZaHt) Hayy], He 
_ 47 
across H,. Let the poloidal non-magnetic forces (including inertial forces) 
acting on the plasma have a density F,, so that they have a component 
H,(F,.H,)/H,? across H,. As remarked in the Introduction, the distortion 
of the field ceases when the magnetic forces become strong enough, so that the 
equation defining the steady-state field is 
VaH,).H FiaH 
Val, — Ce Ha 4 ee (4) 
It is convenient to introduce the following representation (4) in spherical 
polar coordinates (r,@,¢) for a field symmetric about the axis @=o: 


1 oP 1 oP T 


N= —"wein0 00’ = sin or? = * rein” 

where the scalars P and T are functions of r and 6 only. The condition V.H=o 
is automatically satisfied, and P is constant along a line of force. The current 
density j=(c/47)V,H is given by 

1 oT 47. ae: <4 
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We now discuss in general terms the different fields satisfying (6) and the 


boundary conditions. The simplest field is everywhere poloidal, so that T=o, 
and 


7s os 
oJ 


(6) 


FiaH 
UAH, aan tntthe, @ 
p 
Such a field demands that the flow of matter be strictly in meridian planes—each 


element of matter in the cloud must have zero angular momentum about the 
axis of the field. 
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Next, we allow the field to have a toroidal component, but restrict the poloidal 
currents j,, to flow everywhere parallel to H,, so that the field exerts no torque 
about the axis of symmetry. By (5) and (6) this yields the condition 

T=f(P), (8) 


so that T is constant along a line of force; also 


VaH,=f (P)H,. (9) 
In this case the toroidal field contributes to the equation (4) for H,, through the 


term , 
ms won H,- EVO, (10) 


Pp 
t being a unit vector in the ¢-direction. We can now see that such a field, with 
poloidal currents flowing strictly parallel to H,, cannot have lines of force 
extending to infinity. We are assuming that H,, reduces at infinity to the 
uniform field H,, corresponding to 

P+ — }H,y* sin* 6. (11) 
Let us consider the field line, assumed of infinite length, defined by a particular 
value of P, so that far from the star 


(rsind)=( - a) (12) 


‘The toroidal component H,= — T/rsin@; hence if we assert that the field remains 
untwisted at infinity, then 7 must vanish at infinity, and so everywhere along the 
line of force considered. Alternatively, if we assume that the cloud has uniform 
density and temperature at infinity, so that the pressure gradient vanishes, we 
need not require a priori that H, vanishes at infinity; for as r+ along the field 
line P, (r sin @) approaches (12), and the equation (4) to H,, will not hold unless 
f(P)f (P)= 0 at infinity, and so all the way along the line P. If f(P)=o, then 
H,=0, while the case {( P)=constantis ruled out as it yields a singularity on the axis. 
The strictly torque-free field must therefore consist of two parts: an inner 
region, bounded by a surface S, with a toroidal field component maintained by 
poloidal currents flowing parallel to the poloidal field, and an outer region with 
the field purely poloidal. Within S the field lines are closed loops, but they may 
or may not enter the star. There can be no continuous flow of plasma within S, 
as no lines of force cross S. Plasma initially within S will flow down the field 
lines but will not be replenished. Depending on the geometry of the field line - 
considered, it will either be swallowed up by the star, or it will flow towards the 
equator to form ‘‘dead’’ zones; some will be neutralized and flow into the 
star across the field, but there would probably remain a nucleus of ‘‘dead”’ 
plasma, drifting across the field only at the rate permitted by the high conductivity. 
The torque-free field is therefore self-consistent ; if there were a steady flow 
of plasma, it would in general require a toroidal magnetic force to keep the 
angular momentum of the plasma at the correct value for the flow to follow the 
field. Asa further consequence of the depletion of plasma, the field within S 
will approximate to the force-free field, in which all the current flow is 
to the field; for as the plasma drains away, the non-magnetic force F,>o, and 
so only in the dead zones will the steady-state magnetic field exert a force. 
(However, it is pushing the argument too far to assert that a complete vacuum 
results, with zero currents and a curl-free field. A minute electron density is 
quite capable of carrying the currents required to maintain fields of cosmical scale.) 
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We have explicitly excluded centrifugal force as an important contribution 
to the non-magnetic force F,. However, the matter in the cloud will certainly 
have some angular momentum about the axis of the field. We may expect a 
continuous inward flow of matter down infinite field lines, each element tending 
to conserve its angular momentum, and an initially poloidal field H.,, will therefore 
be twisted by differential rotation. The torque so generated opposes the increase 
of angular velocity, but a weak field H,, would still be unable to prevent the 
ultimate build-up of a centrifugal field balancing gravity. The magnetic field 
would certainly play a role in the subsequent evolution of the system, but clearly 
the problem is quite different from that considered in this paper. On the other 
hand, if H,, is strong, a very slight twist will generate a torque strong enough 
effectively to cancel out the increase in centrifugal force. This is assumed to 
be the case; it remains to verify that the slight toroidal field H, generated by 
the inflow exerts a poloidal force weak compared with that exerted by H... 

We obtain as follows a rough criterion for centrifugal force to remain 
negligible. Suppose an element of gas flows a distance dw towards the axis, 
with velocity v,, in time 5w/v,. If the cloud is fully ionized, or, although 
lightly ionized, sufficiently dense for both components to flow closely together, 
then non-uniformities in rotation between two points will be magnetically 
removed in the time that a hydromagnetic wave takes to travel from one point 
to the other. Our criterion is therefore: 


ba ow 

— > ain 

ve 4A,/(47p) 
(p being the density), or 


2 
Ae > Ipo.t; (13) 


the poloidal magnetic energy density should exceed the energy density of the 
inward motion. 

We shall return to this condition in Part II, after the magnetic and dynamical 
fields have been studied. For the moment let us assume condition (13) holds 
for field lines of infinite length, along which there is a continuous flow of matter, 
and note the effect on the equations to the magnetic field. The magnetic torque 
generated by the inflow will be of order 
(VaH,)aH, 

47 


where {2 is the angular velocity, and @ the distance from the axis. Thus the 
order of H, is 


w 


~ pv. V(Qw*) ~ pwr Q, (14) 


Hx Perse (15) 
P 

Poloidal currents perpendicular to H,, make no contribution to the correcting 
term, involving H,, in the equation (4) to H,. However, by continuity there 
must also be a component of j,, parallel to H,, or comparable order, otherwise 
there would have to be current sources on the axis of symmetry; hence the ratio 


of the correcting term in (4) to F,, is by (13), (14) and (15) of order 





pv,,.QH, ~ (PP w2)(47p27.0Q) — pOPa 6 
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We may therefore conclude that if the poloidal field is strong enough to keep the 
centrifugal force density small compared with F,—i.e. if (13) is satisfied —then 
the magnetic field has a poloidal part, hardly differing from the solution of (4), 
together with a small toroidal component. A flow of current along H,, much 
greater than that across it, would again not be consistent with the conditions at 
infinity, and so is not relevant to the case of infinite field lines with continuous 
flow of matter. 

2. The boundary conditions.—At infinity the poloidal function P must reduce 
to (11); this takes care of both components of the field. One more condition 
may be imposed on (4); the continuity of the normal component of the field 
at the stellar surface yields 


H,=FAicos@ (17) 


at r=R, where Ad is the magnitude of the frozen-in, uniform field within the 
star, given by (21) below. However, there will in general result a discontinuity 
in the tangential component, as there is no reason why H, given by the solution 
of (4) should equal —Asin@. Such a discontinuity demands a finite surface 
current and an infinite force density, and so is clearly forbidden. The field lines 
within the star must be distorted near the surface so as to link up smoothly with 
the external field. As a rough estimate of the distortion, suppose the field remains’ 
straight within most of the star, but bends in a surface layer of thickness /. The 
resulting magnetic force cannot exceed the gravitational force; let us estimate / 
by equating them: 


wai FR) Fe) (R) on 


Hence if the stellar magnetic energy is small compared with the gravitational, 
1<R, and the distortion is limited to the surface. The change in H, through the 
surface zone is found from the divergence condition to be ~/H,/R, so that the 
condition (17) remains accurate. But in a star with the magnetic energy, 
comparable with the gravitational—so that the field is at the limit for the star 
to be gravitationally bound—the distortion of the field spreads through the whole 
star. Condition (17) ceases to be accurate, and the more complicated boundary 
condition will react on the external field. Discussion of the internal field requires 
study of the equilibrium of the star and the progressive distortion of the field 
as the star contracts, and is beyond the scope of this paper. However, it is 
reasonable to assume that the star retains an internal dipole component of field, 
of strength approximately A7, so that (17) is a fair approximation to make in 
studying the main features of the external field at different stages of the star’s 
contraction. 

If the external field is twisted, there must be a toroidal component within 
the star. Consider first the force-free fields. As the lines of force emanating 
from the star are of finite length, they cannot steadily transfer angular momentum 
to the galaxy (4, 5); hence in a steady state the field must be torque-free within 
the star as well as without. If initially the internal field exerts a torque, it will 
oscillate about the torque-free state; there will be a periodic redistribution of 
angular momentum along the lines of force, but no continuous net transfer of 
angular momentum outwards. By contrast, field lines of infinite length can 
transfer angular momentum steadily out of the star into the cloud, and the internal 
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field could very well exert a permanent torque. Further discussion must await 
detailed study of the problem of a rotating magnetic system. 

3. The curl-free and force-free approximations.—We have arrived at the 
following possibilities : 

(i) A poloidal field with H,, satisfying (7). ‘The lines of force are of infinite 
length, and there is a continuous flow of matter, strictly in meridian planes. 

(ii) The field approximately as in (i), but with a small toroidal component 
superimposed, exerting a small toroidal force sufficient to keep the centrifugal 
field small, and making a small contribution to the equation (4) to H,. 

(iii) The partially force-free fields, in which the current flow is strictly parallel 
to the field within a surface S, but which is of type (i) or (ii) outside S. 

We arrived at case (iii) by showing that a torque-free field must be broken off 
from the local galactic field; no continuous flow of plasma can result, and so 
the plasma density must become very small nearly everywhere within S. 
However, as pointed out by Liist and Schliiter (6), if the magnetic field is strong, 
very small currents across the field will yield forces comparable with the 
non-magnetic forces acting; the term on the right of (4) will therefore be small 
compared with the largest terms on the left, even if the plasma density were 
maintained by flow from outside S. This argument applies as well to the case (i), 
so that as a first approximation it is tempting to drop the term in F,, in (7). The 
untwisted force-free field is the curl-free field; but this case should be clearly 
distinguished from the vacuum field, for which the curl vanishes through there 
not being any charges present to carry current. In our case the current density 
is not zero, but is given by (7). In using ¥~H, =o to compute the field, we are 
asserting that for large H,,, the term 47|F,,,H,|/H,” is small compared with 
O(rH,)/rdr and —0H,/rd@ which together make up V,H, so that these two 
terms are approximately equal. 

We now discuss the untwisted and twisted force-free fields. 

(a) The untwisted field. The curl-free field reducing to H, at infinity, and 
continuous with H,= Alcos@ at r=R, has components 


H,= | He +(A-H,) (F)"Jeose. 


R\"). 
os [ ~H,+(A- H,)(*) fine 


H,=0, 


and scalar function 


P= —}H,s%sin® al + (FZ ‘“ (=) | (20) 


The internal stellar field AT is given by 


A- HR) (21) 


where R’ is the radius at which the freezing of the field into the contracting 
proto-star is considered to have begun. If A> 3H, the field (19) has a neutral 
point on the equator at the radius 


r= RSE) ; ° (22) 
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Hence once R< R’4/3, the general shape of the field will be as in Fig. 1 so that 
many, though not all, of the lines of force in the stellar flux tube are of finite 
length. 

(b) The class of twisted fields. 'The boundary conditions on the field restrict 
P to be proportional to sin? @ at the star and at infinity ; for simplicity, we 
consider only fields with P oc sin? @ everywhere. We have seen in Section 1 that 
the torque-free condition T = f(P) implies that twisted lines of force cannot extend 
to infinity. Hence beyond a certain surface S, the toroidal function T=o, and 
again the field is approximately curl-free, with 


Pints je int (x * 4). (23) 


A being a constant. Within S, T7=f(P) and from (4) and (9) P satisfies 


or + Gce as = -f(P)f (P), (24) 





with »=cos@; hence f(P)f'(P) o P, and so T={(P)=«P, with « constant. On 
S, T must vanish; hence P also vanishes on S, which must therefore from 
(23), be a sphere of radius r,, with A=r,3. 

Thus outside the sphere r,, the field has the form 





H, = Hgcos(1— (2)’). 
H.=—Hand (: +3(2)’), (25) 
H,=0. ] 


Within r,, the twisted force-free field is of the simple type discussed by Liist’ 
and Schliiter (6) and Chandrasekhar (7), in which the scalar of proportionality 
between current density and field is constant not only along each twisted line 
of force, but over the whole region considered. P satisfies 


OP (i-w)eP_ 
pda” wags a®P, (26) 





and so is of the form (7) 


P= —19(1 —p®)[AJ yo(ar) + BI _gia(ar))], (27) 


where J , ye are the Bessel functions of order + 3/2, and A and B are constants. 
The field components are 





H, = 22088 [AJ y4(ar)+ BI as(ar)], 


H,= — S22 (ria Adyq(ar) + BY_salor)}, (28) 


¥ 


H,= SSR* [Ad ya(ar)+ BI sal ar)} 
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The continuity of H, at R and of H,, H, and H, at r, gives 
— ) + BJ _3)2(«1.) = Oo, 


A= rr [AJ yo(aR) + BJ _so(aR)], 


3H,= lA J so(ar,) + BI’ _s/o(ar,)]. 


T, J 


(29) 





Thus once Ad and H, are given, there is a single infinity of solutions, depending 
on the choice of r,. The condition of continuity of j, at r, is satisfied, as 
(VaH),=aH, vanishes as r, is approached from within. 

The general shape of the poloidal field is shown in Fig. 2. It is seen that 
between R and r, there is another sphere r’ on which H,=H,=o0. This is 
a consequence of the fields H, and H having identical senses. At r, and r’ there 
are neutral points on the axis in both hemispheres. The poloidal part of the 
field vanishes on the equator between r, and r’ at the ring of points given by 


[r?( AJ 39 + BU _g2)]' =0, 


but the toroidal field is non-zero there. Clearly Fig. 2 is the simplest of type (5) 
fields; any number of pairs of zones, each similar to the zone between r’ and r,, 
may be included within the final sphere, beyond which the field is curl-free and 
untwisted. 

An essential feature of the field (4) is the presence of linked lines of force 
within r,. The total magnetic flux across a closed loop of the field H,, does not 
vanish: as the lines of H,, are also the lines of j,,, this follows from elementary 
magnetic theory. In the present case, we have 


H,.dA=——| (VaH,).dA= ——@ H..ds¥o, 
t a =|, ' . v) % =}, sae aie (30) 


where A is the area in a meridian plane bounded by C, a closed loop of H,. This 
is the vital distinction between these fields and those which are just poloidal 
fields with a toroidal component generated by non-uniform rotation, without 
any linked flux (cf. the following note). 

A different choice of f(P) (subject to the condition that H, = —f(P)/r sin 6+o 
on the axis) yields a non-linear equation for P within the surface S (which will 
not be a sphere in general). The resulting fields will be much more difficult 
to express mathematically but will not differ qualitatively from the special field 
given by (25) and (28). 

4. The limitations of the curl-free and force-free approximations.—An immediate 
test of the validity of either approximation can be made by computing the neglected 
term 47F,,AH,/H,”, with the previously computed field H,, substituted, and 
comparing it with the other terms in (7). The non-magnetic forces F,, are given 
when the density—velocity field is known. This problem is postponed until 
Part II but we may anticipate that the plasma density along infinite field lines 
—maintained by continuous inflow—will not in general give rise to large pressure 
gradients normal to field lines, so that the force density F,, is at most of the order 
of the total gravitational force density; also, the densities do not exceed the 
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order of the densities in Bondi’s spherical flow pattern, except in the equatorial 
zone. To order of magnitude, the condition is 
4mpGM H 


et 
PH, < > (31) 


H,? _ GMp. 
— an a (32) 
i.e, the magnetic energy density should exceed the gravitational. 

Consider now the untwisted, curl-free field of Fig. 1. As previously noted, 
the plasma within XZL’ will rapidly drain away, so that in this region the 
approximation should steadily improve, except possibly in equatorial regions, 
where an abnormally high density could be built up, so affecting the field somewhat. 
Outside XLL’ the density remains comparable with Bondi’s values. © Then 
assuming a local galactic magnetic field of order 5 x 10~® gauss in regions of 
density 10 atoms/cm*, it is found that the condition (32) is satisfied except near 
the neutral point. Clearly, if the field strength becomes small, the current density 
must be correspondingly large in order that the magnetic force be sizeable; 
so far from being negligible, the curl of the field has to be enormous, so that the 
curl-free approximation is inconsistent when the field strength becomes low. 

To emphasize this conclusion, let us study conditions near X more closely. 
As the plasma has finite conductivity, there will be a steady drift (v, ), across 
the field, maintaining the current density 


j=o(“"), (33) 


(F,, a i C i *. (34) 


so that by (7), 


As long as H is not too small, then (v;), is very small compared with free-fall 
velocities; in particular, the neglect of the plasma inertia in (7) is completely 
justified. But as the neutral point is approached (v;), increases until the inertial 
force density equals (F,),; the magnetic force is then negligible and (F,), is 
balanced by mass-acceleration, yielding the velocity v; to be used in (33) for 
computing j and so VaH. Thus the objection to the curl-free field is not that 
it demands an infinite curl at the neutral point, but that before the field has 
become low enough to allow a large drift across it, the equilibrium conditions 
demand that |VaH] be enormous, though finite. For example, take F, =p g, 
where p; is the plasma density, and g the stellar gravitational field, so that 


PB, 
hed H?’ (35) 
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and the mass-acceleration 
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This becomes comparable with p,g when 


2pi(p; 

2008) aH on 
°)\ HS 

(3) 

the force pg, is balanced by an acceleration of the plasma within the narrow 
region, on either side of the neutral point, in which the left-hand side of (37) 
exceeds unity. If we take o/c?~10~, corresponding to a temperature of 
~10?°K, and p;=10~‘p(m,/m,), and the same simultaneous values of H, 
and p,, as before, then e.g. when the neutral point in solution (19) is near the 
star’s surface, r|V,H|/H~10*%. As F,, may very well be larger than pg 
(Section 8), r|VaHI|/H could be even larger. 

Thus long before the neutral point is reached, the curl-free approximation 
has completely broken down, and a field initially of the shape of Fig. 1 would 
be immediately distorted by the inward flow of plasma towards the star. As 
long as the plasma and field move together—and this condition requires careful 
study—then the progressive distortion of the field must certainly continue until 
the neutral point reaches the stellar surface, where the gravitational field is 
balanced by a pressure gradient. This motion of the neutral point reduces the 
density of field lines cutting the equator between the star and the point X; in 
order to yield a strong force the field must have a large curl, so that the field 
lines make sharp hairpin bends at the equator. The magnetic force must be 
finite just outside the star, so the X-type neutral point must be dragged well 
below the surface; it is possible it may coalesce with and so annihilate the O-type 
neutral point that field (a) must possess within the star. The general shape of 
this field (a’) is given in Fig. 3. It should be noted that well away from the 
equator, and also beyond _X, the curl-free approximation remains good, and the 
lines of force are similar to the corresponding lines in Fig. 1. 

However, in discussing the transition from field (a) to field (a’), we have 
dealt only with the requirement that the radial forces should balance. It is clear 
that the distortion of the field must introduce new transverse magnetic forces. 
The sharp curving of the field lines near the equator yields both a large current 
density, in the positive toroidal direction, and a strong, approximately radial field 
on either side of the equator. Hence there is a strong magnetic force acting 
towards the equator, which must compress the plasma, and build up a density 
gradient of the right magnitude to withstand the magnetic force. The situation 
is very similar to the well-known ‘‘pinch-effect’’. As indicated in the 
Introduction, although in equilibrium the magnetic and non-magnetic forces 
must balance, it does not follow that the magnetic field plays a passive role: as 
in the present problem, the field may build up, through compression, a thermal 
field able to balance the magnetic forces. 

The motion of the field of type (a) into one of type (a’) is accompanied by the 
progressive destruction of the flux of the closed loops, such as C’CD in Fig. 1. 
The neutral point moves onto loops of smaller and smaller area until the star’s 
surface is reached; at the same time more lines of the type F’E’EF appear from 
the neutral point. This destruction of flux at an X-type neutral point is the 
reverse of the process discussed by Dungey (8) in another context. Near an 
O-type neutral point the drift of plasma through the field, due to the finiteness 
of the conductivity, is always such as to destroy magnetic flux: this is the gist 
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of Cowling’s argument against the dynamo maintenance of fields containing 
O-type neutral points (9g). Near an X-type neutral point, however, the ohmic 
field may either create or destroy flux, according to circumstances. For example, 
if another strong non-magnetic force were brought into play, acting away from 
the star, and stronger than gravity, the field would evolve in the reverse direction, 
new closed lines of force being produced as the neutral point moved outwards 
with the plasma. 

In fact, the field in our problem would never take on the structure of Fig. 1. 
As the star contracts from the initial state with A =lH,, the neutral point would 
appear near the surface of the star when R=R’/,/3, but in the absence of a 
strong force overcoming gravity there is nothing to drive plasma and field outwards. 
A state of equilibrium, with the field having a topology similar to that in Fig. 1, 
could be maintained with the help of a plasma pressure gradient near the neutral 
point, and an equatorial dent in the field lines of finite length so that they can 
hold up the plasma*. Quite apart from the question as to how it could arise, 
such a field would be unstable against lateral perturbations, which would enable 
the plasma to flow down the finite field lines into the star, so that again the neutral 
point would be dragged into the star. But this possibility serves to emphasize 
that the field adopts the structure (a’) rather than (a) because of the absence, 
all the way to infinity, of a non-magnetic force able to balance gravity; we may 
anticipate, for example, that if there is a strong centrifugal field in the cloud, one 
of the effects will be to break off the stellar field from the galactic field, through 
the appearance of a neutral point. 

Further discussion of field (a’) is postponed until Section 6. We now turn 
to fields of type (6), with twisted field lines with a sphere S of radius r,, and an 
untwisted curl-free field beyond S (Fig. 2). There are four neutral points, all 
on the axis, at r=r,,r’, and 6=0,7. Along the axis the plasma can flow 
straight into the star, but near the neutral points at N, N’ the flow beyond S 
must be deviated so as to follow the lines of force such as O’TO, which swerve 
away from the axis. Hence again the field beyond S and near N and N’ must 
have a large curl in order that the magnetic force may overbalance gravity and 
deviate the flow as required, the drift across the field due to the finiteness of the 
conductivity being again too slow for inertia to be important. The field (6) is 
therefore inconsistent with the equilibrium conditions, and the star’s gravitational 
field drags inwards the plasma and field; the neutral points N and N’ coalesce 
with M and M’ respectively and a steady state is reached when the field has 
acquired a structure as in Fig. 4, with the two neutral points at the stellar poles. 
The distorted field lines will exert a force towards the axis, so compressing the 
plasma flowing down towards the star. 

Plasma initially within the sphere r’ (MV M’ in Fig. 2) will flow directly 
into the star and will not be replenished, but plasma initially between r’ and r, 
flows down the field to the equator. Some will be neutralized there and flow 
into the star, but there will probably be formed an equatorial zone of dead plasma, 
unable to flow into the star except at the very slow rate allowed by the finiteness 
of the conductivity. A negative outward pressure ‘gradient in this zone will 
supply a force opposing gravity and friction, and could get over a difficulty at 
the O-type neutral point of the poloidal field (cf. Fig. 4). At O, not only 
(VaH,)sH,=0, but also (VAH,)AH,=f'(P)H,AH,=0, and the magnetic _ 

* This suggestion was made by Mr D. Lynden-Bell. 
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force would seem unable to balance the non-magnetic forces. However, in fact 
the difficulty near O is not very serious and can be got over even without the . 
help of a pressure gradient; whereas in the two previously discussed cases the 
neutral point problem led to striking changes in the fields, a very slight change 
in the toroidal field will yield a finite magnetic force. We write, instead of T=aP, 


T = P{a + B(P—P*)"*], (38) 
where P* is the value of P at O, and f is a function of P, finite at O, and dropping 
rapidly to zero as P deviates from P*. Then we find 


| VaH,|~ P*BP,,*, (39) 


which is finite near O, and so yields with H, a finite magnetic force. In an 
accurate treatment of the field this extra twist would need to be included in the 
equation for P; for our purposes it suffices to remark that Fig. 4 remains an 
adequate picture of the poloidal field. 

The correction to field (6) is less drastic than to field (a). Figs. 2 and 4 are 
alike in the topology of their lines of force; whereas in Fig. 3 all finite lines of 
force have disappeared, and the whole of the stellar flux-tube extends to infinity. 

5. The origin and stability of the twisted and untwisted fields.—The fact that 
we can construct, besides the unique purely poloidal field, a whole class of twisted 
fields that satisfy the conditions of equilibrium and the boundary conditions. 
shows that the steady state condition is not in itself adequate to determine the 
field uniquely. This is because it is only the component of current across the 
field that enters explicitly into the magnetic force; currents parallel to the field 
enter only through the condition 47j=cVaH that holds for all components. 

In the time-dependent hydromagnetic problem there is no indeterminacy. 
The rate of change of the magnetic field is given by dH/dt= Va (v,aH), where 
v; is the velocity of the plasma (1); together with the equation of motion of the 
plasma this determines the magnetic and velocity fields at subsequent times in 
terms of the fields at the initial time. This is not altered if at some points the 
drift of plasma across the field is not negligible, due for example to the curl of 
the field becoming very large. Hence when steady-state solutions have been 
found the attempt must be made to discriminate between them either through 
a stability criterion, or through a qualitative argument showing that one or other 
steady state could not have been reached from the initial state through a physically 
allowed motion of the lines of force. 

Suppose, therefore, that the local galactic field from which the stellar field 
has arisen is unbroken initially, so that its lines of force extend to infinity. (This 
assumption has been implicit in the discussion so far.) Clearly such a field can 
approach an untwisted equilibrium state as in Fig. 3. However, the transition 
from Fig. 3 to Fig. 4 is forbidden, because it demands creation of flux from an 
O-type neutral point. Creation or destruction of flux is dependent on the relative 
drift of plasma and field, which arises through departures from the relation 
E+v,aH/c=o. The correct analogue of Ohm’s law for a gas is 


AH VaH ; 
=: ~s Shed (40) 


470 


E+ 





where o is the conductivity in e.s.u., and E‘ represents the ‘‘ impressed ’’ e.m.f.s 
due to the electron partial pressure, the Hall field and other effects tending to 
cause drift of electrons relative to the ions. These extra effects are small, and 
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in any case take far too long to build up to their maximum, and so the only 
important cause of flux creation or destruction is the Ohmic field cVaH/470; 
and as previously remarked, at an O-type neutral point this always acts so as to 
destroy flux, in contrast to its behaviour at an X-type point. In order to reach 
Fig. 4 from Fig. 3 a circle of neutral points would have to be formed in each 
meridional plane, which would then have to separate into the two bounding 
curves NVN’ and NWN’. This requirement is enough to dispose of the 
possibility. 

However, it is possible that the local galactic field may have initially a poloidal 
component as in Fig. 3, with a broken-off part separating the lines of the stellar 
field and those stretching to infinity. Even so, the field cannot take up the modified 
force-free structure (b’) unless there is also a non-vanishing toroidal flux across 
a region bounded by a closed loop of H,,. The distortion of H, by non-uniform 
rotation will not generate such flux. Closed toroidal loops as are required 
must be generated on the axis of the field, where H,=o. Again, Cowling’s 
argument shows that dynamo action near a neutral point is unable to overcome 
the destruction of energy and flux by Joule dissipation. ‘To generate flux from 
the axis we require a poloidal ‘‘battery’’, represented by E' in (40); but 
again the time-scale is far too long for the process to be relevant, even if E' 
exists. 

Hence if the poloidal field has initially the structure of Fig. 3, but there is 
no toroidal flux, it cannot approach the equilibrium state (6’); any toroidal 
field component within the surface S must exert a torque, and along each loop 
of force toroidal hydromagnetic oscillations will be set up about a mean state 
with H,=o. Consider then the equilibrium of the mean field. Field lines 
within the star are gripped by the pressure and gravitation fields; their external 
parts adjust themselves to the external forces, if any. But the part of the field 
between NV N’ and NWN‘ hasan intrinsic tendency to contract, for (V AH) aH/47 
is directed inwards in the general direction of the neutral point O. If there is a 
thermal field present, of energy comparable with the magnetic energy, then this 
contraction will be prevented by thermal pressure, and the rate of disappearance 
of the poloidal field loops into O will be controlled simply by the very slow Joule 
decay. ‘This is a possibility if the region considered is at a distance not small 
compared with the Bondi sphere, at which the scale height is comparable with 
radial distance (cf. Part II, Section 7). But if the region considered is well 
within the Bondi sphere, any plasma originally within NWN’'VN will drain away 
either into the star or towards an equatorial zone of thickness comparable with 
the scale height, and moreover it cannot be replenished from outside. Hence the 
contracting flux-tubes cannot build up a pressure field able to withstand them, and 
there is nothing to prevent the continuous contraction of lines of force into loops 
of smaller and smaller area around O. Simultaneously, the field lines within NV N’ 
and outside NWN’ will move towards each other. Clearly if there were no 
plasma between NV N’ and NWN’, even in an equatorial zone, the flux-tubes 
would just disappear into O, NV N’ and NWN’ would coalesce, forming a neutral 
line, field lines on either side of this line would link up, and the system would 
ultimately settle down into an untwisted field of type (a’). However, if there 
were enough equatorial plasma, the collapse of the flux-tubes would be halted 
once a sufficiently strong pressure gradient had been built up near O. The 
equilibrium state would then be similar to type (a’), but containing as well a 
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tube of broken-off flux surrounding O and pinching the plasma within it. Such 
a tube would steadily decay through resistance; but in any case, it is unlikely to 
persist for long, as the highly compressed plasma will recombine rapidly, and it 
cannot be supplemented from outside. As the plasma is neutralized, its partial 
pressure decreases, the pinching magnetic force becomes unbalanced and 
compresses the plasma still further, so accelerating the rate of recombination, 
and the disappearance of the flux-tubes into O. Hence the field would approach 
configuration (a’) in a comparatively short time-scale. 

The last possibility is that the poloidal galactic field has initially a structure 
as in Fig. 3, but with toroidal flux linking the field H,, so that the modified 
force-free structure (b’) is the relevant equilibrium state. This equilibrium field 
is likely to be stable against perturbations that are not too violent, even if there 
is no thermal pressure field; for there is no hydromagnetic motion that can 
remove the toroidal field, so leaving an unstable poloidal field. Contraction of 
the poloidal loops towards O is resisted by the toroidal field. If a loop of H,, 
contracts towards O, then H, increases roughly as 1//, where / is a mean distance 
from O to the loop considered ; whereas the toroidal field dragged with it increases 
like 1/7. Similarly, contraction of a toroidal tube into the axis is prevented by 
the poloidal field. 

Given enough time, Joule dissipation will cause a simultaneous contraction 
of the loops of H,, into O and the toroidal loops into the axis, so that again the 
field would be approached. (Quite generally, it is clear that no dynamo action 
can prevent the slow decay of a strictly force-free field, in which there is no 
current flow across the field.) This drift is very slow, but violent enough motions 
could so distort the field as to accelerate the decay. 

In all the discussion it has been assumed that the stellar field is approximately 
uniform and given by equation (21)—it is just an isotropic compression of the 
field H,, assumed originally the uniform field of the whole cloud. We have seen 
that the equilibrium fields arising during the contraction of the proto-star are 
of type (a’); fields such as (6’) are not relevant unless the initial cloud field 
contained toroidal flux. However, if more complicated cloud fields are allowed, 
then equally there is no compulsion to adopt (21) for the stellar field, so that 
one boundary condition can be relaxed. For example, the twisted force-free 
field (b’) may extend from the surface S to the centre of the star. Or the stellar 
field could be uniform, but with a sense opposite to that of the field at infinity Hy. 
If the external field is purely poloidal, then it has a structure similar to Fig. 3, 
except that the region between NVN’ and NWN’ is not present, and the field 
within VV N’ has the opposite sense. The corresponding field (6’) would have 
a similar poloidal part, and toroidal flux within NVN’. 

6. The untwisted field—The computation of field (a’) is a formidable problem, 
and is postponed for separate treatment; however, anticipating some results of 
Part II on the velocity-density field, we may supplement the conclusions of 
Section 4 with some semi-quantitative arguments. 

Plasma flows down the field lines, at a rate dependent on the structure of the 
field, and concentrates in the ‘‘equatorial zone’’ in approximate mechanical 
equilibrium. There is a strong tendency for the neutral matter to be dragged 
by the plasma into the zone. Neutral matter in the zone (including recombined 
plasma) drifts inwards across the field under gravity, but resisted by the 
plasma-hydrogen friction; this friction is thus approximately equal to the 
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gravitational force acting on the meutral gas. Hence in the equatorial zone the 
equation to the equilibrium of the plasma is 


v 
awe — Vpi+pVb=0, (41) 


where ¢ is the gravitational potential GM/r. The plasma pressure p;=2n,kT, 
n; being the ionic number density, k Boltzmann’s constant, and T the absolute 
temperature, the factor 2 arising from assuming singly ionized matter. We 
write n;= nn, where n,, is the number density of neutral atoms (mainly hydrogen) ; 
the factor 7 is given by a balance between the rate of inflow of matter into the 
zone, the rate of drift of neutral matter inwards, and the rate of plasma 
recombination. Thus » will be a complicated function of radial distance from 
the star, but as a first approximation we shall assume it varies slowly enough for 
it to be taken as constant. 
Taking the curl of (41), we have 


7 reed BR Vd=0. (42) 


Also, (41) can be rewritten as 


(VaH),aH 


- —2nkT Vn, +pVb=0, 


so that 
[((VvH)aH]a ve 
47 
Equations (42) and (44) combine to give 


[(VaH) aH], ¥6-+ ZT 94 [(VaH)aH]=0, (45) 
H 

which is the equation to the field in the equatorial zone. Besides the same 
boundary conditions as before, we impose the restriction that H, should be 
negative all the way from the star to infinity, so that there is no neutral point 
outside the star. The equation (45) requires one more boundary condition ; 
this may conveniently be taken as the value of H, at the star’s surface. Assuming 
a plausible value for y, (45) can in principle be solved by substituting for H 
in terms of the scalar P, given by (5). Equation (43) then yields n,, and so also 
the plasma density n;. The densities in the zone should strictly be consistent 
with those given by the rate of flow down the field lines into the zone (cf. Part II), 
the drift of neutral matter inwards, and the rate of recombination. In practice, 
one can only hope, by a trial and error method, to choose value of » and (H,), 
which yield density fields consistent as to order of magnitude. 

Even so the solution of the problem remains formidable. Although the field 
outside the zone is approximately curl-free, its analytic form will not be simple 
because of its having to fit on to the rather eccentric zone field. 

However, in the discussion of Section 4 it was pointed out that the field and 
plasma do not move strictly together, the drift of plasma across the field being 
given by (34). With low field strength and high current densities this drift 
can become important: as discussed in Section 4, once the inertia of plasma 





= 2nkT( Vay, A Vv). 
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drift becomes equal to the inward-pulling forces on the plasma, then the objection 
to the existence of a neutral point disappears. At this stage the structure of 
the field becomes (a”), a combination of fields (a) and (a’). Within the 
field line passing through the external neutral point the field will be nearly curl-free, 
for plasma drifting through the neutral point is soon forced to follow the field 
lines, so that most of this region has a low plasma density as well as a high field 
strength. Outside this limiting field line, the field will be similar to field (a’), 
with hairpin-shaped field lines at the equator, but curl-free field lines away from 
the equator and far enough away along the equator. 

We may anticipate that the appearance of a neutral point will be assisted by 
the pinching effect of the magnetic field. The more distorted the field becomes, 
the greater the compressive force of the field, and the higher the plasma density 
required to balance it. Again at high densities the plasma recombination rate 
goes up. Unlike the case of the pinch effect discussed in the last section, there 
is now an external supply of plasma, as the field lines are of infinite length; but, 
as shown in Part II, the more distorted the field, the greater is the inhibiting 
effect it has on the flow of gas from infinity. Hence at some stage in the 
contraction of the star, further distortion of the field demands for equilibrium 
a plasma density greater than can be maintained. The plasma density decays, 
and the field lines contract laterally ; this accelerates still further the recombination 
of the plasma, and at the same time inhibits still more the replenishment from 
outside. Thus the plasma density must steadily decrease, so that the frictional 
coupling with the neutral gas diminishes, until the dominant inward force on the 
plasma is just the gravitational force pg. Again, the curl of the field, at the point 
where inertial force takes over from magnetic force as the principal opponent of 
gravity, is given by (33) and (35), with H satisfying (37); hence at this point 
|V,H| is approximately proportional to p;}*. Thus as p; decreases, with the 
consequent lateral motion of the field lines and increase in | V4 Hi], simultaneously 
the value of |VaH| required to allow substantial plasma drift across the field 
decreases. Once, therefore, this lateral instability sets in, there is nothing to 
prevent the appearance of an X-type neutral point and broken-off field lines. 
It is hoped that a subsequent paper will determine at what stage in the star’s 
contraction the neutral point appears. 

It is of interest to note the evolution of the field, assuming that the present 
theory were relevant to a magnetic star that has reached the main sequence, and 
so stopped contracting. The field within the star is not curl-free, and so would 
gradually decay; lines of force would slip out of the star, and the external field 
would steadily adjust itself to the changing boundary condition. Ultimately, the 
field would be given by the solution of (7), reducing to Hy at infinity, and 
continuous in both components at r=R with the ultimately uniform, curl-free 
field within the star: the solution of (7) would determine the amount of flux of 
the original galactic field retained by the star. 


Part II 
The dynamical problem 
7. The flow of gas in a fully ionized cloud.—Although we are mainly interested 
in a star embedded in a lightly ionized H1 cloud, we begin by considering a 
fully ionized cloud. The problem could be relevant to the case of a bright star 
surrounded by an H1 region, but the main reasons for studying it are its 
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comparative simplicity, and because even in a lightly ionized gas the plasma and 
neutral components flow close together if the density is high enough. We shall 
assume that the flow is isothermal. In an Hu zone compressional energy 
generated is rapidly absorbed by a slight increase in the already high degree of 
ionization (10); while in a dense H1 cloud a high concentration of dust grains 
acts as an efficient cooling mechanism. ‘The magnetic field assumed is of type (a’) 
or (a’). 

Although the plasma flow is along the field, it should be noted that there will 
not be a continuous flow along all lines of force. This can be seen by comparison 
with Bondi’s spherically symmetric non-magnetic problem (2). There exists 
a critical sphere, of radius r,=GM/c,*, where M is the stellar mass, G the 
constant of gravitation, and c,, the constant velocity of sound. The stellar radius 
must be less than r,, otherwise the star could not hold itself together against 
thermal pressure. Within r, the gravitational energy per gram exceeds the 
thermal energy per gram, so that in a hypothetical state of hydrostatic equilibrium 
the scale height would be small compared with the radial distance; this results 
in the impossibly large densities near a fully contracted star found by Bondi for 
the case A=o (zero accretion), and also in his solutions with subsonic velocities 
near the star. Hence once the stellar radius is well below r,, matter within r, 
must flow in radially at approximately free-fall (supersonic) velocities. Beyond 
r, continuity demands a radial, subsonic flow; departure from hydrostatic 
equilibrium is very slight, and the density is nearly constant. 

Now consider the present problem. If a line of force is part of the star’s 
field, the gravitational pull of the star must again lead to flow down the line 
of plasma within r,, and by continuity the flow extends itself outwards beyond r,. 
Similarly, if the line does not cut the star, but crosses the equator well within r,, 
there will be a continuous flow down the line into the equatorial zone. Hydrostatic 
equilibrium along such a line of force would not be reached until the density at 
the equator is enormous. Quite apart from the loss of plasma from the zone 
through recombinations and inward drift across the field, there will not be 
enough matter in the cloud to supply the density field necessary for hydrostatic 
equilibrium. Hence the picture is as outlined in the Introduction: steady flow 
down the field lines into an equatorial zone containing both plasma and neutral 
matter, their densities being determined by the inflow of plasma, the rate of 
recombination, and the dynamics of the drift of neutralized gas across the field 
into the star. 

By contrast, if a line of force cuts the equator beyond r,, hydrostatic equilibrium 
along such a line is soon reached by a slight drift of plasma towards the equator, 
and no continuous flow will result (provided that the source of radiation keeping 
the cloud ionized is maintained). Hence there is a critical tube of flux that crosses 
the equator near the Bondi radius; outside it the gas is at rest, but within it there 
is a continuous flow down the field lines. 

We now discuss the dynamics of the flow. Once a steady state has been reached, 
the plasma velocity v is parallel to the field: 


v=«H, (46) 


with « a scalar function of position. If p is the plasma density, the equation of 
continuity yields 


O= V.(pv)= V. (px) =H. V(px), (47) 
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so that 
px =constant along a line of force =p,.x... (48) 


Hence in a steady state there is a non-vanishing velocity field v,, = «..H, far from 
the star: even if initially the gas is at rest at infinity, the flow will asymptotically 
approach a state with v,, finite. This is because the magnetic field prevents 
flow of gas from one flux-tube to another; if the field does not vanish at infinity, 
each flux-tube remains of finite width, with continuity demanding a finite 
velocity everywhere along the flux-tube. 

The equation to the steady flow along a field line is just Bernoulli’s equation, 
with the gravitational field as a body force: 

catlog F _ on + $v? = 40,,’. (49) 

In a steady state the magnetic field does no work, as the flow is parallel to the 
field. Introducing non-dimensional variables x, z, y and h, defined by 


GM 
pa ([3)* P=Pa's V=Coy, H= Hh, (50) 
we have from (48) and (49) 
P na 
and a2" (St) 
- - +logz+ $y?= }y,,”. (52) 


$(y)—Y(x)=y, (53) 


Eliminating z, we get 


where 
d(y)=4y*—logy, (x)= =—logh, y=4ya*—logys, (54) 


Equations (51) and (53), applying to each line of force, replace Bondi’s equations 


xyz = | (dylogy)-(Z +2log) = —logA (55) 
for spherically symmetric flow. 
The discussion of (53) follows closely Bondi’s discussion of (55). The 


minimum value of ¢(y) is 0-5, and it is attained for y=1; hence the parameter y 
must be chosen so that for all x along the line of force considered 


y+Y(x)>0'5. (56) 

Now ¢(x)>o0 as x>0; also, for all physically relevant fields, (x) 00 as x0. 
Hence along each line of force must have a finite minimum value. Suppose 
first that %¢>0o everywhere along the line except at infinity. Then (56) is 
satisfied for y>o-5, and so ¢(y)20°5. By inspection of the graph of ¢(y) 
against y, we have 

if Yo>t, then y>i 

if Yo<I, then y<!l for all finite x. 
either y>1 
or y<! 
If y is non-negative for all x, but has a zero for x =x’, then the solution with 


Y~ = 1 will by (53) attain the sonic velocity again at x’; hence a solution 
beyond the zero may change over to a subsonic solution at x’, and vice versa. 


and if Yo=l, then 
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Now suppose ¢ becomes negative for some range of x along the line of force 
considered, and let x,, be the point at which it attains its minimum. Depending 
on how far the star has contracted and on the line of force considered, x,, may 
be a true minimum and off the equator; or just at the equator; or not a true 
minimum, but the point where the line enters the star’s surface. The inequality 
(56) will hold everywhere along the line provided it holds at x,,. Thus if 
y +u(x,,)=0'5, then y(x,,)= 1, and the flow on either side of x,, may be supersonic 
or subsonic. The velocity at infinity is then given by 


$(¥x) “er [o-5 —p(x,,)]; (57) 


approximately, 


Vx = [1 — 2yb(x,,)}** if the flow is supersonic at infinity, (58) 


Va = exp [y(x,,)—0°5] if the flow is subsonic at infinity. (59) 


All values of y,, yielding $(y.,) larger than (58) are allowed, for then ¢=y + y >0°5 
at x,,, and so for all other values of x; but values of y,, between (58) and (59) 
will not yield steady state solutions. If y,, is less than (59), the flow is subsonic 
everywhere ; if greater than (58), it is supersonic everywhere. For either of the 
two values (58) and (59) the flow is just sonic at x,,, and so flow supersonic for 
x >x,, may become subsonic for x<x,, or, by a discontinuity in the velocity 
gradient at x,,, may remain supersonic for x<x,,. Similarly, flow subsonic 
for x >.x,, may become supersonic or remain subsonic for x <x,,. 

However, not all steady state solutions are physically plausible. As % gets 
large for small x, by (53) the velocity in the solutions subsonic near the star 
becomes very small for small x, and leads by (51) to very large densities. This 
is again due to the smallness of the scale height well within the Bondi sphere : 
for low velocity solutions, the equation of motion approximates to the equation 
of hydrostatic support. Such solutions may be relevant when the star’s radius 
is near r,, but once R<r, they can be ignored, along with the zero flow solution. 

Assuming, then, that the flow near the star is supersonic, the problem arises 
as to which of the allowed steady-state solutions will arise in practice. Let us 
suppose that initially the cloud is at rest at infinity. Within the Bondi sphere 
the gravitational field causes inward flow, and the region of flow is gradually 
extended outwards by continuity. Once the system has reached that velocity 
field, supersonic near the star, consistent with the equations of motion and 
continuity, and with the least energy per gram of matter, there is no reason why 
the velocities should increase above this level. Hence we expect the system to 
approach the flow with the minimum allowed velocities: if ~20 everywhere 
along a field line, we take y,.=1, while if % has a negative minimum at x,,, we 
take y(x,,)=1, and adopt the solution that is subsonic for x«>x,,. The point x,, 
on each line of force plays the role of the point x=} in Bondi’s spherically 
symmetric solution (2): beyond it the dynamical pressure is small, and 
approximate hydrostatic equilibrium holds. 

However, it should be emphasized that the equations of motion and continuity 
do not forbid all higher velocities at infinity, with correspondingly greater flow 
along the field lines. If >, all supersonic velocity fields are allowed; while 
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if {<o for some range of x, all values of y,. greater than (58) are allowed. 
Solutions of this type, with y, greater than the minimum allowed, will be 
relevant to the case of a cloud not initially at rest at infinity, but having matter 
being driven in rapidly enough by an external pressure field (e.g. by an external 
H11 region). This result is in striking contrast to Bondi’s spherically symmetric 
flow, and is again due to the confining effect of the magnetic field. If Bondi’s 
spherically symmetric cloud had a non-vanishing inward velocity at infinity, 
the radial flow would rapidly build up a high pressure gradient and dam back 
the flow, which would asymptotically approach the steady pattern with zero 
velocity at infinity. In the present problem the magnetic field prevents the 
build-up of such a pressure gradient, by restricting the flow to be along the 
lines of a field that is finite at infinity. The flow along each infinitesimal 
flux-tube is analogous to flow through a pipe of variable but finite cross-section : 
not only are there low velocity solutions, in which the pressure controls the 
flow until it becomes sonic, but there are also high velocity solutions in which 
the pressure gradient is negligible everywhere. These results are analogues 
of well-known results in laboratory gas dynamics, the reciprocal of the 
field-strength playing the part of area of cross-section of the pipe. 

In Bondi’s spherically symmetric problem there is a definite upper limit 
to the total rate of accretion by a star: 


GM\? 
A~ (SS) CoP = 4120 Pao. (60) 
Again, a higher rate of inflow would lead to a large pressure gradient beyond the 
sphere r,, which would dam back the flow. In the present problem, the rate 
of inflow across a sphere S surrounding the star, with outward unit normal n, 


and of radius less than r,, is 


- | pv.ndS=af \x|pH .ndS= I, Ivo .m)dS, (61) 


where S’ is the northern hemisphere. We have seen that there exist steady 
state solutions with v,, supersonic; these yield a rate of flow across S greater 
than 


Pala i) H.ndS =20L,%0¢., (62) 
H. Je 


where L,, is the radius at infinity of the flux-tube containing S. The total 
accretion rate is given by taking for S the Bondi sphere; hence if L,, for this 
sphere exceeds 4/2r,, the accretion rate exceeds Bondi’s maximum (60). Even 
if L,. were less than 4/2r,, high enough velocities at infinity would still lead to 
enhanced accretion. 

However, supersonic velocities at infinity will arise only if there is an external 
pressure field driving the matter in. We have seen that a cloud initially at rest 
at infinity will spontaneously approach a state with at most the sonic velocity 
at infinity. If the freezing of the field into the proto-star starts at a radius 
R,>+/2r,, then L, > Ry>+/2r,, and so the “‘cross-section’’ for capture of 
material will be greater than the area of Bondi’s sphere; but this will not lead 
to an enhanced accretion rate unless the velocities at infinity are of sonic order. 

Thus the crucial question as to whether a magnetic field can assist spontaneous 
accretion of matter by a star (still assuming that centrifugal force is negligible) 


19 
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depends critically on the structure of the field, as determined by equation (7); 
for it is the field-strength which determines the variation of along a line of 
force. If y% remains positive along the line considered, then y,, = 1 for that line; 
but if the field-strength increases too rapidly inwards along the line (implying 
a sharp decrease in the area of the flux-tube along which the gas flows), ys becomes 
negative, and y,, is by (59) strongly subsonic. Although an accurate estimate 
of the spontaneous accretion rate must await detailed computation of the 
fields (a) and (a’), our qualitative discussion of the field suggests strongly that 
the accretion rate will be well below Bondi’s value. Those lines of force that 
approach the equator are crowded together, as illustrated in Fig. 3; the large 
radial component of H will yield a negative value for 4(x,,), and a consequent 
large reduction in the accretion rate. 

It should be noted that field (@) would also cut down the accretion rate. 
Only a few infinite field lines reach the star; the rest cutting the equator beyond 
the neutral point X. Now the position of X in the curl-free field (a) is near the 
radius at which the freezing-in of the field commenced, and this is likely to be 
at least as great as. Hence most of the infinite field lines cut the equator beyond 
r,, and so there will be no continuous flow of plasma down them. 

If the steady state theory is to be relevant, the time for the system to approach 
its asymptotic state should not be too long. Once the flow has extended outwards 
to a distance x, at which the field lines are nearly straight, we can expect the flow 
to approximate to the steady state pattern: the gas confined by the straight 
field lines will just expand into the partial vacuum created by the inflowing matter. 
For a rough estimate, let us suppose that before the asymptotic state is reached 
the matter flows in with Bondi velocities: the time for the flow to reach x, is 
then ~ }x,3(GM/c,.*)~6 x 10° Mx, where M=M/©. Now x, is the distance 
at which the freezing of the field into the star may be considered to have started. 
If x,~ 1, the time-scale is very short; and even for x)~ 10 (near the upper limit 
if the star is not to contain too much magnetic flux for stability), the time-scale 
is still short compared with the Kelvin-Helmholtz time-scale, and of the same 
order as the time of free-fall. 

8. The flow of gas in a lightly ionized (H 1) cloud.—As in the fully ionized case, 
we must distinguish between lines of force that cut the star or cross the equator 
within Bondi’s sphere r,, and those that cross the equator beyond r,. Again 
there is no continuous flow of plasma down the distant field lines, but instead 
hydrostatic equilibrium is reached after a slight drift of plasma towards the 
equator. By contrast, there is a continuous flow of plasma along the field lines 
within r,, towards the star or the equator, and the plasma velocity field extends 
outwards, leading asymptotically to a finite velocity field at infinity. 

In addition there is the problem of the flow of the neutral gas. In the absence 
of the magnetic field, all the matter would flow in radially with Bondi velocities, 
but in the present problem, with the plasma tied to the field, the flow of neutral 
gas across the field is resisted by the frictional drag. The density of the mutual 
friction between the two gases is (1): 


+nynom,V(v,, —V,), (63) 


where n,,, m, are the number densities of hydrogen and ions, v,,, v, are the velocity 
vectors of the two gases, o is the mutual collision cross-section, taken as 10~#* cm?, 
m,, is the mass of the hydrogen atom, and V =|v,, —v,| or c,, whichever is the 
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larger. If the hydrogen is to flow across the field essentially as in Bondi’s theory, 
the component across the field of the friction, with Bondi’s radial velocity field 
substituted, must be small compared with the gravitational force on the hydrogen : 


Nyt (Vq »(*) (su) = #(Z) (S22) ng(r6 x 104M) <1. (64) 


Ca ~ 
We have assumed n/n, = 10~+, a standard value for the ratio of the densities of 
ionized metals to neutral hydrogen in an H1 cloud. 

Beyond r, (x=1), the Bondi velocity is = c,,/x* so that V=c,, and ny = (mq) .0- 
Condition (64) is then independent of x, and yields 

(mao POE. (65) 
Within r,, however, Bondi’s velocity ~c,./x"? is supersonic; V is then |v, —v;, 
which is of the same order as vy, except for regions where the field is nearly radial 
(where in any case both gases flow parallel to the field). The ratio (64) is then 
oc1/x"%, and so friction increases more rapidly than gravity as the star is 
approached. Whenever the frictional drag on Bondi flow exceeds gravity, the 
transverse flow must in fact be given by a balance between friction and gravity, 
and is therefore cut down to well below Bondi’s velocity. 

When friction cancels the transverse component of gravity, the neutral matter 
is effectively forced to follow the field lines. Hence for densities above the 
limit (65), the two gases flow together down the field, and the theory of the last 
section applies to the gas as a whole. In particular, both gases have asymptotically 
a finite velocity parallel to the field at great distances from the star. As Bondi’s 
spherical velocity field vanishes at infinity, and the transverse drift of the neutral 
gas is reduced below Bondi’s field by friction, the approximation of taking the 
neutral gas flow as parallel to the field improves at great distances. 

However, although for (nq)..>6 x 10°/M the flow of neutral matter within 
the critical flux-tube is everywhere closely parallel to the field, there is a slight 
transverse drift, such that the frictional drag on the neutral gas balances gravity. 
This flow will be maintained by drift of neutral matter across the field from the 
‘*dead ’’ plasma zone outside the critical flux-tube. If the flow within the tube 
were strictly parallel to the field, the neutral matter in the ‘‘dead’’ zone would 
reach hydrostatic equilibrium as well as the plasma; as it is, continuity demands 
some inward flow. Thus the neutral matter that reaches the equatorial zone 
consists mainly of matter that has flowed down the field lines within the flux-tube, 
but includes some that has drifted. in from the ‘‘dead’’ plasma region. 

If (n,).. satisfies (65), frictional resistance to flow across the field is negligible 
beyond a radius r, <r,; in particular, neutral gas will flow into the critical flux-tube 
from the ‘‘dead’’ region, and will continue radially with Bondi velocities until 
r,, after which it will flow closely parallel to the field lines. Within the critical 
flux-tube the plasma builds up its finite velocity field far from the star, and this 
flow is gradually transferred to the neutral gas by mutual friction. Hence 
ultimately both gases flow nearly parallel to the field far from the star; but the 
time-scale for this state to be built up is probably longer in this case. 

Assuming 1,/n,=10-*, we find the friction everywhere exceeds the 
gravitational force on the plasma if (nq)..>0°6/M; for lower densities friction 
dominates near enough to the star. Depending on the density at infinity, and 
the distance from the star, the correct non-magnetic force F,, to be used in (7) 


19° 
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for an accurate computation of the field is the gravitational force on the plasma 
or the frictional drag on the plasma, whichever is the greater; but the friction 
cannot exceed the gravitational force on the neutral gas. 

g. The centrifugal field —Consider finally the condition (13) for the magnetic 
field to be able to keep the centrifugal force small compared with gravity. Beyond 
the Bondi sphere r,, the spontaneous velocity field is at most sonic, so that (13) 
certainly holds if the magnetic energy density is large compared with the thermal, 
as is the case for the galactic magnetic field usually assumed in H1 clouds. 
Within the Bondi sphere the inflow is approximately free-fall; v2? ~GM/za, 
and (13) is just the condition (32) for the magnetic energy density to exceed the 
gravitational. Again, the magnetic field is strong enough for this to hold in 
general: in fact it is the condition which would imply the curl-free approximation, 
if this did not lead to a neutral point. If the cloud field were somewhat weaker 
than is usually assumed, then the two energy densities would be roughly equal 
—the initially straight field lines would be distorted by the inflowing matter 
until a rough equipartition of energy is built up. 
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A NOTE ON THE MAGNETIC BRAKING OF A ROTATING STAR ~ 


L. Mestel 
(Received 1959 March 7)* 


Summary 
A poloidal magnetic field connecting a rapidly rotating star and a 
slowly rotating envelope will be twisted so that it transports angular 
momentum across the intermediate region of low density. The field exerts 
equal and opposite torques on the star and envelope. In the intermediate 
region the field will probably not exert strong compressive forces, but will 
be approximately force-free. By contrast with the completely force-free fields 
of the previous paper, which cannot be generated by distortion of a poloidal 
field, there is no toroidal flux linking a poloidal loop, and current does not flow 

in the same sense at every point of a loop. 





In a recent paper (1) List and Schliiter have discussed the magnetic braking 
of a rotating star, in the context of the von Weizsicker cosmogony. It is assumed 
that a star has been formed in a rotating cloud; a magnetic field emanating from 
the star ensures that a sphere of radius rg surrounding the star rotates with it 
(the Ferraro—Alfvén effect). If Q) is the star’s angular velocity, rg is taken 
as of the order of the ‘‘ Kepler’’ radius, at which gravity and centrifugal force 
balance : 

Q,°7 3 = GM, (1) 
M being the stellar mass. Beyond rz each sphere has the appropriate Keplerian 
angular velocity Q, where 
273 = GM =Q,27,3. (2) 
The non-uniform rotation is supposed to lead to turbulence, and the consequent 
rapid dissipation of energy to a simultaneous outward transport of mass and 
angular momentum. 

The form of the magnetic field beyond 7, is not discussed by the authors, 
as they assume that turbulent friction is the dominant toroidal force responsible 
indirectly for the outward flow of mass and angular momentum. They concentrate 
attention on the torque-free field between the star and rz; the field is supposed 
to have both poloidal and toroidal components, the twist being sufficient to 
transport from the star to rg just that amount of angular momentum that the 
turbulence beyond r, is transporting to infinity. The rate of loss of angular 
momentum by the star is estimated as a function of the decaying turbulence, 
but not of the magnetic field, which plays a subsidiary role: the turbulence 
beyond rx disperses the envelope, and generates just the right twist in the field. 

The authors also discuss the possibility that the edge of the uniformly rotating 
zone may be given not by (1), but by r,,, defined by the condition of equality of 
magnetic and turbulent energy densities. In this case the magnetic field does 
appear in the formula for the loss of angular momentum, but only through the 


* Received in original form 1958 May 16. 
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quantity 7,,; the picture is still of a turbulent transport beyond r,,, and a 
torque-free field between the star and r,,, its twist being determined by conditions 
beyond r,,. 

The authors give an elegant discussion of the conditions under which a 
torque-free field may transport angular momentum. If 7;; is the magnetic stress 
tensor, and k,; the magnetic force density, then 


: I 
T= - = 3H%5,;), (3) 


and 


Pa 
= Ty (4) 


t 
T,; may be interpreted as the rate of flow in the 7-direction of the j-component 


of linear momentum, and the integral | 7;,;n; dS as the total flow of j-momentum 


across a Closed surface S with outward normal m,;. Similarly the authors introduce 
a tensor density describing the torque-density of the magnetic field: 


C) r) . . 
d;;= (x;k;— x,k;) = = Ox, Die = — Fy, Cin Tie)» (5) 


where D,;, is the density of flow of angular momentum. The index-pair (i,/), 
(with Dj;,= — Dj.) gives the direction of the angular momentum component, 


and k the direction of flow. The integral | D,;,2, dS over a closed surface S 


s 
gives a unique value for the net rate of transport of the (7,7) component of angular 
momentum out of the volume enclosed. For an axially symmetric field defined 
by scalars P and 7 as in (2), the total flow out of a sphere of the component 
parallel to the axis is 


I moe 


ef TS d= MRO. 1))- FPG, (6) 


where .=cos 6, and F(P) is defined by the torque-free condition 


a _ aF(P) ig 
P=f(P)= 5. (7) 
It follows from the condition that the field be source-free that this flow is 
independent of r, as must be the case in a steady state. 

The relevance of Liist and Schliiter’s work does not necessarily depend 
on the correctness of the cosmogony with which it is linked. Certainly, if the 
Keplerian velocity field does not give rise to turbulence, then an important 
feature of their process will be missing. However, the magnetic field beyond rz 
can play a role analogous to turbulence, in that constant twisting by the Keplerian 
field will generate toroidal forces and lead equally to a radial transport of matter 
and angular momentum. Even if one does not picture a star as being born in 
a medium which already possesses a strong rotation field, the magnetic field will 
transfer angular momentum from the proto-star to the cloud and in time will 
probably build up a rotation field similar to that assumed by Liist and Schliiter. 
If the loss of angular momentum from the star is not rapid enough, compared with 
the rate of contraction, the star will take up a disk-like structure, symmetrical 
about the rotation axis. The subsequent evolution of the disk may very well 
be determined by the magnetic redistribution of mass, the inner parts contracting 
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and the outer parts dispersing with excess angular momentum (3). The relevant 
feature is that there will again arise a zone in which the magnetic field exerts no 
torque and so does not accelerate or decelerate the matter in the zone (which 
will be of low density), but acts as a means of transport of angular momentum 
from an inner to an outer region—-just the problem discussed by Liist and 
Schliiter. 

The purpose of the present note is to criticize the particular field chosen 
in (x) and to suggest more plausible alternatives. The authors choose 


nn HEEB) 


where P,= P(r, + 1) (assumed non-zero), and 7, is a constant to be determined 
by the conditions at rg. The toroidal function T is then 


T=F(P)= +47. (z) (5 S ‘). (9) 


T can have no discontinuities, and F(P) must have opposite signs at the poles 
if the field is to transport angular momentum; hence there must be a line of force 
P=o so that F(P) can change sign. The simplest function satisfying this 
condition and the other conditions for a real non-singular field is 


P= P, cos? 6. (10) 


Hence the field chosen, and the currents maintaining it, have components 


2P, e T, 7 
H, = —* cos 8, ie= — Fe se (3008? 8-1), 


8a 

H,=0, Je= (11) 

H,=— 22sin 6 cosé, j= £ SBsine, 

This field transports angular momentum 4 7,P, per second from the star to the 
radius r,. The meridional lines of force are straight lines through the origin. 
It should be noted that the field (11) is not force-free; the authors point out 
that in some cases considered the magnetic forces are far larger than the 
gravitational, and if the conditions of equilibrium (equation (4) of (2)) is to be 
satisfied, it must be through a large plasma pressure gradient, built up by the 
magnetic forces themselves. It was pointed out in (2), Section 5, that there is 
a limit to the amount of compression a plasma can tolerate, especially if 
neutralized plasma cannot be replenished from outside, as will be the case with 
a Keplerian rotation field beyond ry. Further, field lines that are to assist in 
decelerating a star must be part of the stellar field, and there is nothing to prevent 
plasma initially within r, from flowing into the star. Hence a field such as (11), 
with lines of force having opposite senses on either side of the equator, is likely 
to approach a state with an X-type neutral point and broken-off field lines. The 
high conductivity of the plasma is no argument against such a transformation, 
any more than in the work of (2). The value of the conductivity merely determines 
the scale of the distortion required for relative drift of plasma and field to be 
noticeable; if either the magnetic or non-magnetic forces produce a field with 
a large enough curl, the theorem of the freezing of the field is no longer valid. 
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In fact, the reasons for adopting the rather eccentric field (11), with its 
singularity, are not very compelling. The authors justify their choice by appeal 
to an ‘‘economy-principle’’. They suggest that the field is determined by the 
transport of angular momentum, in the sense that the field adjusts itself to transport 
angular momentum “‘ with least expense’’: the field lines link the origin and r, in 
the simplest way, without being anywhere retrograde. But this principle can 
hardly be more fundamental than the condition of equilibrium. If the equilibrium 
condition and the boundary conditions at rg do not suffice to fix the field, then the 
remaining arbitrariness may be removed by appeal to the economy-principle ; 
but the choice of the field (11), in spite of the forces it exerts, rather puts the cart 
before the horse. 

Without denying that cases may arise in which the magnetic forces are 
balanced by gas pressure, it is certainly plausible to look for a force-free field 
that can transport angular momentum across the region of probable low density 
between the star and rg. The equation to the field is (4) of (2), with Fp~o. 
The field is also torque-free, but it should be noted that relation (7) holds only 
between the star and r,; within the star (7) cannot hold, as the field must exert 
a decelerating torque. In particular, a line of force (defined by a constant P) 
which enters the star at 6 (> 7/2), travels through the star and leaves it at (7— 0) 
is not forbidden in its external parts to have associated toroidal components of 
opposite sign in northern and southern hemispheres, as is required if the line 
is to assist in the transport of. angular momentum. Outside the star the field is 
torque-free up to rg, and so 7 must be constant along each segment of the line 
of force within ry and outside the star. Within the star, however, T will gradually 
change its sign as the line traverses the star. Similarly, beyond rg, T will vary 
along a field line in such a way that the field exerts on the matter a torque equal 
and opposite to the torque exerted on the star. Liist and Schliiter’s difficulty 
in finding a force-free field capable of transporting angular momentum seems to 
arise from their demanding that 7 be the same function of P along the whole 
length of a line of force. 

In the preceding paper we discussed force-free fields, with 7 constant along 
each field line, and hence with toroidal flux across loops of the poloidal field. 
Such fields cannot be maintained by dynamo action against Joule decay, nor 
can they be generated by distortion of a poloidal field by non-uniform rotation : 
any linked flux must be primeval. However, for the present problem this type 
of force-free field is irrelevant, as we do not demand that 7 be constant along the 
whole length of a field line. As there must be currents across the field within 
the star and beyond r,, in order that the field may exert a torque, the equation 
of continuity of charge can be satisfied without T having the same sign everywhere 


along a field line. Hence even if there are closed field lines C of Hp, | jp-ds 
c 


can vanish, so that there need not be any linked toroidal flux: nor is there any 
contradiction in picturing the field as arising by toroidal distortion of an initially 
poloidal field. Further, dynamo maintenance of the field against Joule decay 
cannot be ruled out a priori. In the regions where the field exerts a torque, 
dynamo action can drive the transverse currents, while as the currents parallel 
to the field need not flow in the same sense at every point of a line of force, they 
can be driven by an electrostatic field. (Of course general arguments such as 
Cowling’s theorem apply to forbid maintenance of certain types of field; “but 
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the breakdown in dynamo action is due to the induced currents being too weak 
near an O-type neutral point, not to their having the wrong direction). 

It is of course possible that the initial field is not poloidal, but for example is 
one of the twisted fields of (2), with toroidal flux. In that case non-uniform 
rotation will superpose a further toroidal component, and the problem is 
essentially unchanged. However, we shall consider only the case of an initially 
poloidal field, with axis identical with the axis of rotation, and with a toroidal 
part imposed by the condition at rz, where either turbulent friction or a magnetic 
torque begins to act on the matter. Suppose, for example, that the poloidal field 
in the star and in the torque-free region were as in Fig. 3 of (2), up to the sphere 
rg, and let P* be the value of P on the line of force that just grazes the star at 
the equator. Consider the following function T defining the toroidal field : 


(a -*) 
To (—}i = -1 o>P>P* 
(zs r and @<z/2 


re) P*>P 


re) P*>P 


and @2>7/2. 
-To(Fs)(ps -*) o> P>P* | 





T vanishes at @=o and on the line P=P*. The function F(P) (yielding 


+aTeP*(5-1) (H+) o> P>P*, 


T=F'(P)) is 


as r 


aaa. (13) 
° P*>P, 


and the angular momentum transported per second—by those lines of force for 
which P > P*—is by (6) }7,P*. The poloidal current maintaining the toroidal 
field is 


1) Hp o> P>P* 


P*>P. 
The poloidal force due to H, is 


T,? P P P 
- <p (pe -*)( gs) (Fe 1) He — Hed o>P>P* 
(15) 
° P*>P 


and the term in H, in the equation to the field ((4) of (2)) is 


(me -1)(F-1)/rsino o> P>P* 


° P*>P. 
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The equation to P within ry is then 


ep sindd(_1 aP\__T#(P\(P _,\(2P_, (17) 
or "7? 96\sind 00)  P*\P*/\P* oll : 


for o> P>P*. For P*>P the equation is (17) with the right-hand side zero, 
provided the density is low in this region also. The boundary conditions are 
continuity of H, at the stellar surface and at rx, and also continuity of H on the 
field line P=P* (so that the magnetic pressures may balance). Numerical 
integration of the equations will decide whether the field can in fact have the 
postulated structure, with all the field lines infinite. 

The conditions will be relaxed somewhat if there is a high plasma density 
for P< P*, so that the field in this region is not curl-free, but exerts a com- 
pressive force balanced by the plasma pressure gradient. As the lines with P< P* 
are not linked with the star, this plasma will not rapidly drain away. However, 
too great a compression would again lead to recombinations and decay of the 
plasma density unless it were replenished from beyond ry. This in turn requires 
that the Keplerian velocity field beyond rz is restricted to the region P > P*. 

If 7, is small enough (corresponding to a slow dispersal of the envelope 
beyond r,), the term on the right of (17) is small, so that for P> P* the field is 
effectively curl-free. If also the pressure gradient is negligible for P< P*, the 
field as a whole is nearly curl-free within rg; its structure will then certainly 
not be as in Fig. 3, but rather as in Fig. 1. Along loops such as CDC’ 
which lie within rg, the function 7 must be zero, as there is no toroidal flux; 
hence no net angular momentum transport can occur, as is physically obvious. 
But lines that cross rg, whether they form closed loops or go off to infinity, will 
again not be torque-free beyond rg, but will be twisted by the non-uniform 
rotation so that they exert a torque of sign opposite to that exerted on the star ; 
hence 7 will have opposite signs at points such as G and G’. The toroidal 
function (12) can now be applied to this field, provided P* is defined as the 
parameter of the limiting loop of force that just grazes the sphere rg. In contrast 
to the previous case, only a fraction of the stellar field lines assist in the transport 
of angular momentum. 

Clearly there are many generalizations of (12). If the field has no toroidal 
flux (i.e. it is just a twisted poloidal field), the toroidal component is determined 
if the conditions beyond rz are completely known. It may be that in no case 
of interest does the field have the structure of Fig. 3, and that the transport of 
angular momentum is primarily from the parts of the star near to the rotation 
axis, a result of possible importance during the contraction of the star to the 
main sequence. This note is merely concerned to show that the general theory 
of Liist and Schliiter can be applied to the type of field that naturally arises from 
equilibrium conditions, and no special fields such as (11) need be postulated. 


St John’s College, 
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STELLAR GROUPS, IV. THE GROOMBRIDGE 1830 GROUP OF 
HIGH VELOCITY STARS AND ITS RELATION TO THE 
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Summary 

The available proper motions and radial velocity data have been used 
to establish the existence of a moving group of subdwarfs (Groombridge 
1830 group) which includes RR Lyrae. On the basis of the relationship 
between the observed ultra-violet excess and displacement below the normal 
main sequence, the subdwarfs in the Groombridge 1830 group are identified 
with main sequence stars in the globular clusters. This identification gives 
a modulus of m—M=14™-2 for the globular cluster M13 with the result 
that My~+o™-5 for the 2 RR Lyrae variables and My=—2™-3 for 
the brightest stars in the cluster. RR Lyrae itself, for which we derive 
My~+0™-8 from the moving cluster parallax, is shown to obey the 
period—amplitude relation for the variables in M3 and to be reddened 
by o™-o5 with respect to those variables. By equating the luminosity of 
RR Lyrae to the mean of the variables in M3 we obtain a modulus of 
m—M=15™-o for the cluster. We have not derived this modulus in the 
logical way of fitting the M3 main sequence to the main sequence of the 
Groombridge 1830 group because the colour observations of the M3 main 
sequence stars may contain a systematic error. 

Because the presence of RR Lyrae variables in stellar groups may 
provide the only accurate calibration of the luminosities of these stars, 
it is important to make a systematic search for such groups. We have 
tabulated the available information for (1) all stars classified as subdwarfs 
on the basis of their spectral features, (2) stars that are probably extreme 
subdwarfs, on the basis of their observed ultra-violet excesses, although 
classified as normal main sequence objects, and (3) stars that may be 
globular cluster giants, subgiants or horizontal branch stars. Space motions 
of objects in these three categories have been computed from parallaxes 
derived by fitting them to the (My, B—V)-diagram of M13, which was 
calibrated by the Groombridge 1830 group. It is shown that there is a 
high probability that several of these stars are travelling together in space 
and that, when additional material becomes available, several groups 
similar to the Groombridge 1830 group will be found. Also, the space 
motions of 4 RR Lyrae variables with relatively well-determined proper 
motions are identified with those for some nearby subdwarfs with a resulting 
My ~ +o™-6 for the four variables. 





In previous papers of this series (Eggen 1958a, b,c; Papers I, II and III) 
several groups of stars have been discussed on the basis of the apices and total 
amounts of the stellar motions. Aside from the interesting dynamical questions 
raised by the existence of these groups (Woolley 1959), the resulting 
(M,, B—V)-diagrams are of importance in supplementing the 
material from galactic clusters. This is especially true if stellar groups can be 
found that contain stars for which a luminosity calibration is not available from 
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other sources. Of particular interest in this connection are the field RR Lyrae 
variables and subdwarfs that might be identified with globular cluster stars. 
In this paper we will discuss (1) the evidence for the existence of a group of very 
high velocity subdwarfs moving together in space and containing RR Lyrae 
(Groombridge 1830 group), (2) the use of this group in obtaining distance 
modulii for globular clusters by the identification of the main sequences of the 
clusters with those of the group stars, and (3) the possible existence of other 
similar groups. 

1. The Groombridge 1830 group.—Among the stars known to be within 
20 parsecs of the Sun, and for which radial velocities and proper motions are 
available, the object with by far the largest space motion is the subdwarf 
Groombridge 1830 (HR 4550). ‘The relevant data are as follows: 

p= 7"'042/yr, U = — 263 km/sec, 
p= —98 km/sec, V = —151 km/sec, 
m=0"-116, W=— 22 km/sec, 
where the velocity vectors, relative to the Sun, are directed as follows: 
U, away from the galactic centre (/=148°, b=0°); V, in the direction of 
galactic rotation (J=58°, b=0°); W, toward the north galactic pole (6=go°). 
Taking the velocity of the Sun, relative to the stars within 20 parsecs (Woolley 
1958) as U(O)= —11, V(O)= +17 and W(O)= +7 km/sec, the space motion 
of Groombridge 1830 relative to the mean of these nearby stars becomes 
U + U(Q)= —274 km/sec, 
V +V(O)=— 134 km/sec, 
W+W(Q)=— 15 km/sec. 

A search through a catalogue of 700 high velocity stars yielded the following 
stars which appear to share the large space motion of Groombridge 1830; the 
stars were chosen by the convergent point (i.e. moving cluster) method described 
in Papers | and II: 

V,=285 


a ADe p 
Star 6, A@sinA mn donlios ies 


gy 
+ 72°94 309 + 5 0°25 — 266 —6 o-o10 
Gmb. 1830 145 ° 7042 — 98 —6 0126 o-116 
HD 163810 209 © 0°82 +192 -—8 0-020 0-018 
RR Lyr 210 4 0°24 — 72 —3 0°004 0°0035 
HD 219617 = 203 I 1°30 + 15 —3 00225 — 4 00205 


The values of A@ were formed from the observed position angles of the proper 
motions, 4, and the position angles computed from the assumption that the stars 
are moving with Groombridge 1830 toward A=16"19™ and D= —54°-4; 
A is the angular distance between the star and the point (A, D). The observed 
radial velocities, p, and their quality, O, were taken from R. E. Wilson’s 
catalogue and the computed values were obtained from the expression V,cosA, 
with two values for V,, the total velocity referred to the Sun, determined from 
two values for the parallax of Groombridge 1830; 285 km/sec for a parallax 
of 0”:126 and 300 km/sec for a parallax of 0-116. The group parallaxes, 7,. 
were computed from the expression 7, = 4-74 4/V,sinA. 

The proper motion and radial velocity determinations for the group members 
are in Table I where the uncertainties in some of the data are apparent; the 
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proper motions of + 72°94, HD 163810 and RR Lyrae, and, to a lesser extent 
the radial velocity of +72°94, HD 163810 and HD 219617, could be greatly 
improved by additional observations. 


Tas.e | 


Proper motion and radial velocity determinations 
for members of the Groombridge 1830 group 


Star Me PE He PE s* S* 


wo a” 


+ 72°94 —0'240 o°'O10 +0°218 G 
—0°195 10 +0°157 G 

Gmb. 1830 3°995 I — 5800 GC W (11) 
L (5) 


Md (2) 


HD 163810 —0°473 W (3) 
—0°321 Md (2) 


GC Md (300) 
Ss W (20) 


HD 219617 —0°512 8 — 1°186 7 GC Ww (8) 
¥ 


Md (2) 


—O°513 O'014 —1'200 o’o10 
—0'512 —1°190 


* GC=General Catalogue, Y=Yale Zone Catalogues, C=Cape Observatory, 
G =Greenwich Observatory, and S=Sproul Observatory. The sources of the radial 
velocities are Md=McDonald, W=Mount Wilson, and L=Lick Observatories; 
the observatory abbreviations are followed, in parentheses, by the number of plates. 


The space motions of the group members have been computed from the 
parallaxes derived above on the basis of the group motion and are listed in 
Table II. In those cases where there is a large difference between two values 
of the observed proper motion or radial velocity, space motions have been 
computed for both determinations to demonstrate the sensitivity of the three 
components of the space velocity to variations in the observed data. For example, 
because of its apparent place in the sky, the large observed difference in the 
distribution of the total proper motion of HD 163810 between the two coordinates 
is reflected almost entirely in the value of the W-velocity. From the spread of 
values of the space velocities in Table II it can be seen that, within the uncertainties 
in the data, the members of the Groombridge 1830 group can be considered to 
have isoperiodic galactic orbits (Woolley 1958) and, therefore, to form a 
kinematical group. 

In Paper II an attempt was made to estimate, by the use of a synthetic 
convergent point, the number of spurious group members that might be present 
in groups of various velocities as a result of the method of selection. It was 
estimated that two or three spurious members would be included for groups 
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with a total space motion of 100 km/sec. The same material (FK3) indicates 
that no spurious members will be included in groups with total space motions 
in excess of 150 km/sec, but because the stars in the FK3, which were tested 
by the synthetic convergent point (6= —go°) have generally very low velocities 
compared to the subdwarfs, this result may not be applicable to the present case. 
On the other hand, the fact that relatively few stars have space motions near 
300 km/sec makes it highly unlikely that spurious motions of this size will 
simultaneously (1) result from placing the stars on a definite sequence in the 
(M,, B—V)-diagram, (2) have vectors U, V and W similar to those for 
Groombridge 1830, and (3) represent stars that have the same ultra-violet 
excess as Groombridge 1830 (Section 3). For example, there are 33 F-, G- and 
K-type subdwarfs in Tables VI and VII that have ultra-violet excesses comparable 
with those of Groombridge 1830 and the globular cluster stars. Only seven 
of these—the first three stars in Table VI, the second in Table VII and the three 
members of the Groombridge 1830 group for which U—B values are available 
—are inbound toward the galactic centre with velocities greater than 200 km/sec. 
Since the values of V+V(Q) and W+W(Q©) for these seven stars range, 
respectively, from —120 to —450 km/sec and from —140 to +60 km/sec, it 
seems unlikely that nearly fifty per cent of the stars would accidentally satisfy 
the requirements for membership in the Groombridge 1830 group. More 
quantitative statistical tests for spurious members of high velocity groups will 
be possible when extensive astrometric, photometric and spectroscopic 
observational programmes for subdwarfs, which are now in operation, are 
completed. 

The magnitudes and colours, which were observed in either the (P, V),-system 
or, by Roman (1955), in the (B, V)-system, are listed in Table II together with 
the spectral types, also by Roman, and the range in luminosity resulting from the 
range of parallax shown in the table. ‘The colours on the (B, V)-system are used 
here since our ultimate aim is to compare the group stars with globular cluster 
stars and the latter have mostly been observed on that system. The colours and 
luminosities of the group members are shown in Fig. 1 where the Pleiades~Hyades 
main sequence is also plotted. The estimated uncertainties in the luminosities 
of the group members, based on the extreme range of the parallaxes shown in 
Table II, are indicated by the length of the lines in Fig. 1. The trigonometric 
parallaxes of the group members, other than Groombridge 1830, are small and 
offer little check on the derived luminosities. The individual values are as 
follows; the observatory abbreviations and the weights, in parentheses, are 
those given by Jenkins (1952): 


Star a, (unit = 0” -oo1) 
+ 72°94 o M(6) 
Gmb. 1830 111 A(28), 107 M(8), 164 St(3)— mean = 116 
HD 163810 26 M(7), 11 Y(12), 31 C(7) —mean=20 
RR Lyr 11 A(28), —12 Yk(6), 5 W(7), 25 S(7)—mean=9 
HD 219617 14 M(8), 39 ¥(12), 43 D(4), o C(7)—mean=24. 


2. The double stars.—Fig. 1 shows that members of the Groombridge 1830 
group are undoubtedly subdwarfs lying 1™-5 below the Pleiades-Hyades main 
sequence in the (M,, B—V)-diagram. The question of the existence of these 
and other subdwarfs in the (Mj,,,, log T,)-diagram is discussed elsewhere in 
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Fic. 1.—The (My, B—V)-diagram for members of the Groombridge 1830 group. The 
length of the vertical lines indicates the range of luminosities for the various values of the parallax 
given in Table II. The continuous curve represents the Pleiades—Hyades main sequence. 


Monthly Notices (Sandage and Eggen 1959). From the equations developed in 
that discussion it is apparent that if subdwarfs do exist in the (M,, ,, log 7,)-diagram, 


they will violate the mass—luminosity relation for normal dwarfs and it is, 
therefore, very important that masses of subdwarfs be independently determined. 
As pointed out previously (Eggen 1956) there are very few known or suspeeted 
subdwarfs in binary systems. The best mass determination available is that for 
85 Pegasi A, but since the derived mass depends upon the trigonometric 
parallax, there is necessarily considerable uncertainty in the value. On the 
other hand, membership of visual binaries in stellar groups permits the use of 
the group motion parallax in deriving the mass. Representative observations 
of the two binaries in the Groombridge 1830 group, ADS 10938 and 16644, 
are shown in Table III. Very few observations of ADS 10938 have been made 
since its discovery in 1900; the three listed in Table III show an angular motion 
of 1°/year and a fixed separation near 0”-5 but orbital elements obviously cannot 
yet be derived. The observations of ADS 16644 extend back to 1876 and more 
than one-half an orbital period may have elapsed since then, but unfortunately 
this system is seen nearly edge-on (inclination near go”) and any orbital elements 
must depend mainly on the observed separations. The two components of the 
system are of nearly equal magnitude and are indistinguishable, but a double 
reversal of quadrant (180° of motion) has been assumed to have occurred between 
1g02 and 1929. ‘The apparent orbit is nearly rectilinear with one end as yet 
unobserved. Using the observations to 1930, Parenago (1946) derived the 
value of P=145 years and a=o"-52 from the observed rate of change of the 
observed separations and the assumption that van den Bos’s observations in 1928 
define one end of the rectilinear apparent orbit. These elements do not satisfy 
later observations and it doubtful whether even preliminary elements can be 
derived until it can be shown from observations made now that apastron has 
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Tasie III 
Observations of ADS 10938 and 16644 
ADS 10938 ADS 16644 
PA Sep. . t PA Sep. 


218 1876°28 42 0°83 
198 1891°59 49 «0°63 
1898 -66 48 058 
1901°16 42 058 
1901°73 47. «(0°63 
1927°00 ‘Too close 
1927°80 180 o715 
1929°60 252 oO'15 
1932°33 Too close 
1937°53 224 0°26 
1942°56 226 0°34 
1943°55 229 0°47 
1944°84 223 0°40 
1948°64 222 0°44 
1949°85 224 0°42 
1955°62 226 0°72 
19§7°90 223 0°68 
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been reached. An attempt to find the extreme range of periods and semi-major 
axes that will fit the available observations indicated that the total mass of the 
system is less than two solar masses, that is, the two (assumed) equal masses are 
both less than the mass of the Sun, but this result is very uncertain. 

3. Identification of the subdwarf group members with globular cluster stars.—In 
addition to the dynamical questions raised by the possible existence of a moving 
group of subdwarfs near the Sun, the existence of RR Lyrae stars in this and 
possibly other stellar groups (Section 5) provides a' powerful method of calibrating 
the luminosities of the cluster type variables themselves. Furthermore, since 
the Groombridge 1830 group shows that subdwarfs and cluster variables are 
related, it seems probable that the globular cluster main sequences should be 
compared to the subdwarfs rather than to the ordinary dwarfs near the Sun. 
In the past this has been a question of some importance because, as Fig. 1 indicates, 
an error of more than a magnitude can result in the luminosities of the globular 
cluster stars if the cluster main sequences are fitted to the Pleiades-Hyades main 
sequence rather than to the subdwarfs. In this section we use the main sequence 
of the Groombridge 1830 group (subdwarfs) to compare with the main sequence 
of M13 and discuss the resulting luminosities of the M13 stars. A similar 
discussion of the available M3 data is given in Section 4 but in this case greater 
reliance is placed on the absolute magnitude of RR Lyrae itself as determined 
from the assumption that it is moving with the Groombridge 1830 group than on 
available observations of the M3 main sequence stars. Justification for comparing 
RR Lyrae with the M3 variables comes from the similarity of the characteristics 
of these stars (Section 4). The divergence of the M3 main sequence from that 
of the Groombridge 1830 group, which results from this comparison (Fig. 8), 
is interpreted as an error of about o™-1 in the colours of the faint stars in M3. 

The most reliable data for faint globular cluster stars are probably those 
published by Baum, Hiltner, Johnson and Sandage (1959) for M13._— Earlier 
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results for M13 by Arp and Johnson (1955) provided accurate three-colour 
observations to visual magnitude 17-5 and the later results carried this to 
magnitude 20°5 in three colours and magnitude 22 in two colours. The two 
colour indices, U—B and B—V, obtained from more than one photoelectric 
observation by Arp and Johnson for giant and subgiant members of M13 are 
shown in Fig. 2 as open circles. The filled circles in the figure represent the 
photoelectric colours obtained for main sequence stars in the cluster by Baum, 
Hiltner, Johnson and Sandage. The continuous curve represents the relationship 
between U—B and B—V for the Hyades stars (Sandage and Eggen 1959). 
The following values of U — B are also available for three subdwarf members of 
the Groombridge 1830 group: 


Star B-V U-B 
m m 
Gmb. 1830 +0°74 —O'17 
HD 163810 +0°625 — 0°06 
HD 219617 +0°47 —0'20 
these are shown in Fig. 2 as crossed circles. The faint companion to HD 163810 
probably has little effect on the colour of the system: the equal components of 
HD 219617 probably have identical colours. 
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Fic. 2.—Three colour diagram for the giant and subgiant members of M13 (filled carcies), 
main sequence members of M13 (open circles) and Groombridge 1830 group members (crossed 
circles). The continuous curve represents the (U—B, B—V)-relation for the Hyades main sequence 
Stars. 

It is evident from Fig. 2 that the observed ultra-violet excess of the 
Groombridge 1830 group members and the M13 stars is very similar. 
Since we have shown elsewhere (Sandage and Eggen 1959) that, at least 
to a first approximation, displacement below a standard main sequence in 
the (M,, B—V)-diagram is directly correlated with the observed ultra-violet 
excess, we will fit the main sequence of M13 to that defined by the Groombridge 
1830 stars. The result of this procedure is shown in Fig. 3. It has been 
assumed that M13 is free from reddening since the bluest stars are observed to 
be bluer than any stars known in the Galaxy. Arp (1955) and Arp and Johnson 
(1955) have derived a reddening of o™-16 and o™-12, respectively, but we believe 
that the procedures used in these derivations are probably unsatisfactory. 
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The filled circles in Fig. 3 represent mean points from the photographic 
observations by Baum, Hiltner, Johnson and Sandage (1959) for luminosities 
fainter than M, = +2™ and from those of Arp and Johnson to +3™; the open 
circles represent the individual photoelectric results for cluster members observed 
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Fic. 3.—The (My, B—V)-diagram for M13 (circles) and the Groombridge 1830 group 
(vertical lines). The filled circles represent photographic mean points and the open circles 
individual photoelectric results. The two crosses represent the RR Lyrae variables in M13. 
The modulus of M13 was determined as m—M-=14:2 by fitting the main sequence cf the cluster 
to the subdwarf members of the Groombridge 1830 group. M13 is assumed to be unreddened. 


by Arp and Johnson. A mean, colour—luminosity relation for the stars fainter 
than M,= +5™ in Fig. 3 is tabulated in Table IV based on a distance modulus 
of 14™-2 obtained from the above procedure. 


Taste IV 
Mean colour-luminosity relation for the main sequence of M13 
(B-V) (P-V),z My 


m m m 
+0°45 +0°565 +5°1 
+0°50 +0°605 +5°5 
+0°55 +0°655 +58 
+0°60 +0°705 + 6°05 
+0°65 +0°75 +6°3 
+0°70 +080 +6°55 
+0°75 +0°85 +6°8 
+080 +0°90 +7°° 
+0°85 +0°945 +7°15 
+0°90 +0°99 +7°4 
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One consequence of the above procedure is that the brightest stars in M13 
are at M,= —2™-3. This value should be compared with the usually adopted 
one of — 3™ which is based on (1) the observed fact that in most globular clusters 
the brightest giants are three magnitudes brighter than the horizontal branch 
(RR Lyrae variables) and (2) the assumption that this branch is at M,=o"™-o. 
The value of M, = — 3™ for the bright end of the globular cluster giant sequence 
has been used in many estimates of galactic and extra-galactic distances and 
therefore, if the procedure in fitting the main sequence of M13 to that of the 
Groombridge 1830 group is correct, these distances will need a substantial 
revision. 

A second consequence of the identification of the main sequence of M13 
with that of the Groombridge 1830 group is that the RR Lyrae variables have 
M,~ +o™-5 instead of M,~o™-o. Additional evidence on this point will be 
discussed in the following sections. 

4. RR Lyrae and its relation to the M3 variables.—The only cluster for which 
there are available three colour data for a large number of variables is M3 
(Sandage 1959). ‘The relation between the periods and light amplitudes, B, 
V, and U, for 33 variables in this cluster is shown in Fig. 4 where the separation 
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Fic. 4.—The period—colour relationship for the variables in M3 more distant than 180" from 
the cluster centre. The position of RR Lyrae and T Sextantis and of two variables in M13 
(crosses) are also shown. 


of the Bailey type a and type 4 stars (large amplitudes) and the type c stars 
(small amplitudes) is apparent. The three colour results for RR Lyrae (Hardie 
1955) and for 'T’ Sextantis (Tifft and Smith 1958) are also plotted in the figure 
and indicate that, in light amplitude at least, these two field variables are very 
similar to the M3 stars. The results for RR Lyrae and T Sextantis are as 
follows : 
U-B B-V Ay Ap Ay 
RR Lyr +0°05 +038 0°95 I°IS 1°08 
T Sex + 0°06 +0°24 0°82 0°55 0°44 
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On the other hand, the available two-colour data for the M13 variables (Arp 
1955), which are represented in Fig. 4 by crosses, show marked deviations from 
the mean period—amplitude relations for the M3 variables. Both of the M13 
variables stand off, in opposite directions, from the mean period-amplitude 
relations for the ab and the c type M3 variables more than any individual variable 
in M3. This result is a particular example of the fact that the period—amplitude 
relation can vary markedly from cluster to cluster (Belserene 1954). This 
variation is probably related to the fact that small variations in the mean luminosity 
of the variables may occur from cluster to cluster (Sandage 1958) and to the 
great variation in the spectral features observed in field RR Lyrae stars 
(Iwanowski 1953, Preston 1958). Therefore the field variables probably 
represent an inhomogeneous group and it is only through the period-amplitude 
relations discussed above that we can identify RR Lyrae with the variables 
in M3. However, before the final comparison can be made between RR Lyrae 
and the M3 variables, the possible reddening of the former must be investigated. 

We have assumed that the stars in both M3 and M13 are unaffected by 
absorption on the basis of their high galactic latitudes and the presence of stars 
observed to be bluer than any known galactic stars. On the other hand, 
RR Lyrae lies nearly in the galactic plane (6= +11°) and, at a distance of 
250 parsecs, may be reddened. The values of B— V and U—B for the variables 
in M3 are plotted in Fig. 5 as crosses; non-variable stars on both sides of the 
variable star ‘‘gap’’ in M3 are also shown in the figure as dots. Also, the mean 
values of U—B and B—V for T Sextantis and RR Lyrae, both of which were 
shown above to be similar to the variables in M3, are plotted in Fig. 5 as large 
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Fic. 5.—Three colour diagram for the M3 B-v. oO2 4 334s 
variables in Fig. 4. T Sextantis falls on the mean 


relation for the cluster stars but RR Lyrae needs Fic. SS ae es ee 
to be moved 0™-04 in U—B and o™-o5 in B—V tion for the M3 variables in Fig. 
along the reddening line indicated by the arrow. T Sextantis (unreddened) and RR jes 
Non-variable horizontal branch stars are indicated (with a reddening of o™-05 indicated by an 
by filled circles. arrow) are plotted as open circles. 


filled circles. ‘There appears to be little or no reddening affecting T Sextantis, 
because it falls on the mean (U—B, B— V)-line defined by the Mg stars, but 
RR Lyrae needs to be moved o™-05 in B—V and o™-04 in U—B along the 
reddening line (U—B)=0-74(B—V) (Lindholm 1957). The value of 
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E.y — p= +0™-05 derived in this way differs considerably from +o0™-12 derived 
previously (Sandage 1956) but, within the framework of the assumptions that 
RR Lyrae is similar to the M3 variables and that the latter stars are unreddened, 
the present value is undoubtedly more reliable. The previous determination 
rests on the questionable assumption that the intrinsic colour of RR Lyrae is 
the same as that of the supposed unreddened, non-variable stars with which 
it was compared. The consistency of the result that RR Lyrae is reddened by 
+o™-o5 and the assumption that it is similar to the M3 variables is demonstrated 
in Fig. 6 which contains the colour (blue)—amplitude relation for the M3 variables. 
A correction for +0™-o05 reddening in this diagram puts RR Lyrae on the mean 
colour—~amplitude relation of the cluster variables. 
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Fic. 7.—Two colour diagram for M3 giants and subgiants (filled circles) and main sequence 
stars (crosses). The broken curve represents the mean relationship for M13 giants and main 
sequence stars and the continuous curve that for the Hyades main sequence stars. 


The above discussion shows that the period-light amplitude relation is the 
same for RR Lyrae and the M3 variables, and it would therefore appear logical 
that their absolute magnitudes might also be the same. This could be verified 
by fitting the main sequence of M3 to that of the Groombridge 1830 group and 
comparing the resulting mean absolute magnitude of the cluster variables with 
that of RR Lyrae derived on the assumption that it is a member of the group. 
However, as mentioned above, with the exception of M13, the photometric 
observations of faint stars in globular clusters are not of high precision because 
the photographic interpolation usually depends on very few photoelectrically 
calibrated standard stars and it is difficult, in the photographic interpolation 
method, to control the colour indices at faint light levels with an accuracy greater 
than o™-1. This is particularly true of M3 (Johnson and Sandage 1956) where 
only 8 photoelectric standards were observed fainter than V =19™ whereas 
23 such standards were obtained for the investigation of M13. Because an error 
of only o™-1 in the colour indices of the faint M3 stars will cause an error of o™-4 
in the position of the main sequence, the M3 data are considered to be of 
insufficient accuracy to compare main sequences. A further indication of possible 
difficulties in the M3 photometry is the anomalous U — B values derived for M3 
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main sequence stars shown in Fig. 7. The accurate M13 data show that the 
main sequence stars have an excess of about 6(U—B)=o™-2 whereas the 
observed M3 values rise to o™-4 to o™-5. 

Therefore, instead of fitting the main sequence of M3 to that of the 
Groombridge 1830 group, we have fitted the mean magnitude of the M3 
variables to that of RR Lyrae in Table II, after correction for the absorption 
discussed above. The resulting (M,, B—V)-diagram for M3 is shown in 
Fig. 8 which also contains the colours and luminosities of the Groombridge 
1830 group members; the colour and magnitude of RR Lyrae have been 
corrected for absorption before plotting in Fig. 8. 
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Fic. 8.—(My, B—V)-diagram for M3. The modulus of 15™-o was derived from a comparison 
of the variables with RR Lyrae. The vertical lines represent members of the Groombridge 1830 
group. The mean colour and luminosities of representative variables in M3 are plotted as crosses. 





The (M,, B—V)-diagrams for M13 (shaded) and M3 (unshaded) are 
compared in Fig. 9 where the two variables in M13 are shown as crosses and 
the variable star gap in M3 by broken lines; the blue extension of the horizontal 
branch in M3, which is populated by only a few stars, is shown as a thin line. 
Figs. 8 and g show that the stars in M3 fainter than V=19™ do indeed 
systematically deviate from the main sequence of the Groombridge 1830 group 
and of M13. This deviation is not considered to be real however, because of 
a possible error in the measured colour indices of the faint M3 stars. These 
same stars, fainter than V =19™, also show an anomalous ultra-violet excess 
(Fig. 7) which may be caused by the same error. 

















Fic. 9.—A comparison of the (My, B—V)-diagrams for M13 (shaded) and M3 (unshaded). 
The two crosses represent the variables in M13. The systematic deviation of the M3 main sequence 
from that of M13 in this figure and from that of the Groombridge 1830 group in Fig. 8, is believed 
to be due to an error in the observed colours. 


5. Other possible groups.—The above discussion of the fitting of the main 
sequences of globular clusters to that of the Groombridge 1830 group is 
intended as an illustration of a procedure that, if it can be extended to other stellar 
groups, might well afford the only definitive method of calibrating the luminosity 
of the globular cluster stars. It is therefore important that a search be made for 
other groups and in this section the data available for such a search are discussed. 


” 


Stars have been labelled ‘‘ subdwarfs’’ on the basis of a variety of criteria; 
for example, velocity, position in the HR diagram, the appearance of the spectrum, 
and the spectral type-colour relationship. The interrelations of some of these 
criteria are discussed elsewhere in Monthly Notices (Sandage and Eggen 1959) 
and we will deal here only with subdwarfs so classified by the appearance of their 
spectra, with some additional stars included on the basis of an observed 
ultra-violet excess in their colours. Stars for which radial velocity observations 
are available and which, as far as can be determined, have been classified as 
subdwarfs on the basis of their spectral features, are listed in Table V. The 
table contains the following information : 

Wilson No.: The number in R. E. Wilson’s (1953) catalogue of radial 
velocities. A few stars, not in that catalogue, are identified by their HD or BD 
numbers. 

m,: The visual magnitude taken from Wilson’s catalogue except for stars 


© 


observed photoelectrically, which are indicated by two decimal magnitudes. 
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Sp.: The spectral type assigned by the observer mentioned in the sixth 
column. Other classifications are noted at the end of the table. 

p,Q: The radial velocity and its quality, a, 6, c, d, taken from Wilson’s 
catalogue. A few values not in that catalogue, or improved by more recent 
observations, are indicated by the observer; R=Roman (1955), E=Evans 
(1958), G=Greenstein (1954). 

Ob.: The authority for the classification listed in the third column; 
P = Popper (1942, 1943), J = Joy (1947), R= Roman (1955), M = Miinch (1944), 
E = Evans (1958) and G=Greenstein (1954). 

u, 9: The total annual proper motion and its position angle, taken from 
various sources. 

B—V, U-—B: Colours observed by Roman (1955) on the (U, B, V)-system 
A few values, given in parentheses, have been (1) converted from the 
(P, V)g-system when no U — B is listed, (2) observed at Herstmonceux with the 
20-inch reflector or at the Cape Observatory with the 18-inch reflector, or 
(3) observed by the Cape observers (Stoy 1958) with the McClean refractor— 
in the latter case the (U — B), on the refractor system is given in the notes to the 
table. , 

This list of subdwarfs is undoubtedly very incomplete, even among stars 
brighter than visual magnitude 10 and radial velocity, astrometric and photometric 
observations are now in progress to increase the number of known subdwarfs 
as well as to strengthen the data in Table V. 

On the assumption, discussed elsewhere (Sandage and Eggen 1959), that 
displacement below the main sequence is dependent upon the ultra-violet excess, 
we have computed space motions for those stars in Table V whose (U—B) 
and (B—V) colours place them on or near the M13 (and Groombridge 1830 
group) line in Fig. 2. The parallaxes used in these computations were derived 
by fitting stars to the (M,, B—V)-diagram in Fig. 3. In addition to the number 
in Wilson’s radial velocity catalogue, or BD/HD number, Table VI includes 
(1) the values of the two components of the proper motion, (2) the radial velocity, 
(3) the adopted parallax and (4) the three components of the space velocity freed 
of the solar motion with respect to the nearby stars (cf. Section 1). 

There are several objects in Roman’s (1955) catalogue of high velocity stars 
that, although not noted there as subdwarfs, have an ultra-violet excess similar 
to the globular cluster and Groombridge 1830 group members. A few of these 
stars are globular cluster giants. Others, although spectroscopically classified 
as luminosity class V, may be subdwarfs. Some of these stars are listed in 
Table VII where, in addition to the information supplied for the stars in Table VI, 
we have included the magnitudes, colours and spectral types; the parallaxes 
were obtained by fitting the giant or main sequence of the (M,, B—V)-diagram 
in Fig. 3. Two A-type stars, which show an ultra-violet deficiency, are included 
in Table VII. The parallaxes of these stars have been computed on the 
assumption that they belong to the horizontal branch of the (M,, B— V)-diagram 
for globular clusters since a similar deficiency has been observed for A-type stars 
in a globular cluster (Sandage and Walker 1955). Although there may be a 
range of o™-3 in luminosity of the horizontal branch from cluster to cluster 
(Sandage 1958), we have fitted these A-type stars to that of M13 in Fig. 3. 

In Table VIII we have gathered together the stars in Tables VI and VII 
that appear to have similar space motions. It should be remembered that in 
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TaBLe V 
Stars classified as subdwarfs on the basis of their spectral features 
Wilson No. My Sp. p pe tf] B-V U-B 
g1* 12°0 sdF4 — 100 0°25 79 ies ee 
390 10°27 sdAg —114 0°35 72 +0°36 —o'16 
646* 5°20 Gs VI — 97 3°76 115 +0°70 +O'1! 
987 10°08 sdF2 — 266 0°32... += +0°44) 

HD 16031 9°75 sdF4 + 63 O17 153 +0°45 —o'22 
1735 8-06 sdF7 — 139 0°82 194 =+0°46 —0'24 
1800 10°5 sdK5 — 167 1°360)=—-156 see ae 
1828 10°18 +154 0°3: 140: (+0°46) (—o-25: 
1892 10°8 — 33 0°45 # I4!1 is oe 
1904 10°6 —173 1°58 135 wa ea 

HD 21543 8-22 3 + 59 o'41 120 +0°62 +0'03 
1919 9°7 — 160 1°63 «131 San beat 
2420 9°22 +339 0°54 122 +0°44 —O'lg 
2851 6°97 ; 066 146 +060 +0°23 
3314 II's 0°45 #8126 ae se 
3562 II'9 060 86170 
3668 10°8 0°43 165 cee nee 
4005 8-08 0708 «225 4=%+0°86 +0°22 
4400 10'0 056) 175 
5058 11°6 0°62 «141 es oes 
5240* 8°34 1°97 158 +0°60 —o'12 
5494 9°69 066 182 +048  —o'19 
573° 9°59 0°62 142 (+0°42) (—0°25) 
5757 11's 0°43 ~—«180 lal ee 

25°1981 9°35 0°35 197 +0°30 —o'12 
5784* 8°74 009 «149 «= +0°907)S_ +:0°53 
5793 10°10 0°37. 247 (+0°38) (—0°24) 
5815 10°8 0°64 155 sae sae 
5908 II's bas 
5909 9°8 
5951 11°8 

+9 2190 10°5 
6236* Ir°o 
6237 10°2 
6294* 8°24 
62906 10°6 
6704 10°3 
728 12°6 
6734 8°23 
6869 10°8 
6880 98 
7017 113 
7117 10°70 
7i3I 9°5 
7170 6°49 
7266 
7296 
7308 
7404 
7547 
7507 
7913 
7923 
7932 
7951 


2) 
a 


MO OS OS TT I TU AD TUM ADVE VO | 


°°79 138 
0°47 ~=—«180 


0°54 ~=—« 160 


0°88 +0°4I —0°25 
o18 180 oie 


(+0°49) 


+0°47 —o'1g 


(++i 41 


(+0°53) 
+0°74 o'l7 


+0°48 —o'18 
+0°44 —0'25 


l+t++ei ++ i 
aoanneene nen gThonn nn gTanngTenQHanasrerwananvaageoeaTonDPeaorTrwmeoe oO0O 
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TaBLe V (cont.) 
Stars classified as subdwarfs on the basis of their spectral features 
Wilson No. Me Sp. p Ob...» 6 . B-¥ 
8236 10°04 sdF3 — 162 0°54 161 +0742 
11°6 sdG3 + 16 088 164 oe 
6°20 sdFsp — 22 021 46290 «=©6—6+ 0°90 
9°7 sdF7 +126 O50 «215 ins 
10°46 sdG2 —138 087 270 +0°70 
9°74 sdF4 + 34 035 185 +042 
8-5 sdAg +167 0-75 229 ei 
9°06 Ko VI +292 3°68 196 +080 
7°20 sdF'5 —171 I'IQ 254 sees 
10°0 sdF5 —157 Osi 213 
11‘4 sdG2 +167 0-44~—sa183 
9°13 sdG1 — 51 0°29) 0.229 +067 
+2°3375 9°95 sdFs5 — 389 0°37 +0°44 
10287 6°98 sdA2 — 363 0°08 +0°16 
10289 10°8 sdF8 + 9g! 044 «148 igen 
10367 II‘o sdF8 —216 0°60 te 
10404* 9°62 sdF8 +191 0°83 0°62 
11060 11*4 sdAs — 68 o'26060=— 184 
HD 181743* 9°64 sdF + 18 082. 180 (+0. -46) 
11949 10°20 sdF7 —172 0°53 +0°48 
+ 26°3578 9°36 sdF4 —128 o'17 +0°37 
12730 8°5 sdGo —172 1°21 sted 
12754 10°6 sdFo —319 
13327 738. FoVI - @ 
13543* 10°52 sdF +103 
13782 11‘7 sdF6 — 67 
+ 17°4708 9°47 sdF6 — 296 
14122* 14°4 sdKs5 — 157 
14539 10°00 sdG2 —113 
14568 10°62 sdF6 — 32 
14598* 8-17 sdF8 + 10 
14704 10°5 sdFo — 100 
14884* 7°OV r — 1§ 
15042 8-92 + 50 


Notes to Table 


0-90 +0°51 
0°48 (+0°40) 


+0°43 
1°68 sis 
nl (+0°48) 
(+0°42) 
I 6 +0'47 


o"90 obo 
0°34 (+0°33) 


BaeoaracaRoarnanwDomhonerRMonnocganonnnaecd 
Se UD VU DV SRPSSKPOMOZ—SARvvvPD vy Peeu vy 


S 


91 LDS 1A. 
646 p Cas. 
4005 Roman (1955) classifies this star as G2 III and states: ‘ All lines are 
very weak’’. 
5240 Roman (1955): “* very weak lines ”’. 
5784 Roman (1955) assigns the classification Gs IV and states: “ very weak 
lines and strong CH ”’. 
6236 Spectroscopic binary with a range of 83 km/sec. 
6294 Roman (1955): sdF2. 
8554 Roman (1955): G6 V. 
8691 Astronometric binary, P~13004, 
8771 A companion, 300” distant, has common proper motion and velocity. 
9840 Roman (1955): G2 V. 
10404 ADS 10938. 

HD 181743 (U—B)co= +1™-49, which is roughly equivalent to U- B= —o™-17. 
13543 (U—B)cg= +1™-45, which is roughly equivalent to U- B= —o™-20. 
14122 Spectroscopic binary with a range of 37 km/sec. 

14598 ADS 16644. 
14884 SX Phe. 
15042 ADS 17173. 
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Tasie VI 


Space motions of stars in Table V with observed ultra-violet excesses 
comparable to that of the M13 and Groombridge 1830 group stars 


Wilson No. Me Be p Tp U+U(0) V+V(O) W+W(o) 


8771 —1°OI2 —3°544 +300 0°037 — 282 — 450 
5240 +0°725 —1°836 -—240 0°035 —273 —211 
8236 +0'104 —0°438 —162 0-008 — 204 —177 
1735 —0'203 -—o'805 —139 0°026 — 163 — 103 
6294 +0°370 —0o°'790 — 18 ovolg — 160 — 135 
8603 —0°03I —0°347 + 34 O°'OIO0 — 100 
74604 —O'OI2 —0°473 +159 O°'OI2 — 195 
5784 +0049 —0°082 —168 00017 —272 
5730 +0°373 -—0°488 +203 o-o12 —270 
13327 —O'IIZ —0'900 — 45 0°034 — 86 
HD 16031 +0°078 —O'155 63 o-o12 
7308 —0o'184 —0°430 95 Ov’oll 
6734 —0'268 —0°465 62 0°025 
HD 181743 -—0°'077 —0°818 18 0013 
HD 21543 +0°354 —0°204 59 070395 
4005 +0°051 —0°056 167 0°003 
8249 —0°203 —0°'062 — 22 0°0055 — 132 
5494 —o'o10 —0°648 + 59 0°0135 — 187 
10287 —o'050 —0°064 —363 0°006 — 285 
+17°4708 +0°496 +0°071 —296 0:'0135 +198 —250 
+ 2°3375 —0°368 +0'069 —389 o-or! + 338 —213 
2420 +0°455 —0'283 +339 o'o14 + 344 — 109 


l++++ 


ie +1 


each case the space motion was computed on the assumption that the stars could 
be fitted to the dwarf (or giant) sequence of the (M,, B— V)-diagram for M13 
and the Groombridge 1830 group. This assumption is based on the observed 
ultra-violet excess in the colours of the stars. When more photometric and radial 
velocity material becomes available for subdwarfs, we believe that additional 
stars will undoubtedly be found to be members of these possible groups. In 
support of this view we include in brackets in Table VIII a few stars from Table V 
for which photoelectric, photometric data are not available. For these stars we 
have tried several parallaxes, guided by a rough spectral type—colour relationship, 
until one was found that gave space motion vectors similar to those for some 
star in Tables VI and VII. It should be emphasized that this method is exactly 
equivalent to the convergent point method described in Papers I and II since, 
for an individual star, a necessary and sufficient condition for group membership 
is that the components of the space motion be identical with those of the group. 
An illustration of the accuracy of the parallax which can be obtained by this 
method, on the assumption of group membership, is given below for Wilson 
Nos. 6294 (HD 84937) and 5909 (—3°2525). This illustration, as well as the 
discussion of the Groombridge 1830 stars given in Table II, also shows the 
sensitivity of the computed space motions to the accuracy of the fundamental 
data. 

It should be noted that, although we have stated that group membership 
implies that the space vectors of the group stars are identical, this may not be 
precisely true. Deviations from this requirement can ocur, in U and V, of 1 to 
2 km/sec/100 pe as a result of small differences in the eccentricities of isoperiodic 
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Wilson No. 6294 No. 5909 


” 


Ha= +0°380 = j4g= —0°790 Ha=+0°33 we=—0'58 Cape 
= +0°359 = —0'791 


+0°370 —0°790 p=+20 W (6) 
= +40 Md (2) 

— 18:1 W (5) 

—15°3 Md (1) +25 


p 


—17°6 


7 U+U(0O) V+V(O) W+W(O) 7 U+U(0) V+V(0O) W+W(O) 
0018 — 168 — 147 +10 o°o10 — 241 — 202 —5§ 
[o-o19 — 160 —135 o] O°Ol4 —172 — 145 +1 
0°020 —153 —128 —10 [o-o15 —159 —135 +2] 

0°016 — 148 —128 +4 


orbits of widely separated stars (Woolley 1958). Also, if W is not tightly coupled 
with U and V, we would expect its value for the group members to be correlated 
with the height of the stars above or below the galactic plane because of the 
differences in the gravitational field at different Z-distances (Woolley 1957). 
However, an evaluation of this effect for the Z-distances involved here shows 
that extreme differences of only 10 to 15 km/sec in W will be encountered and 
in this preliminary investigation they will be ignored. 

One of the most important uses of the possible groups discussed above is 
in calibrating the luminosities of field RR Lyrae variables. As is evident in 
Table I, even the available data for members of the Groombridge 1830 group 
are uncertain and this is true of subdwarfs and RR Lyrae variables in general. 
There are not more than a dozen RR Lyrae variables with well determined annual 
proper motions in excess of 0”-o5 and, perhaps, another dozen with relatively 
large but insufficiently observed motions. The same situation exists for the 
radial velocities which in many cases depend upon only one or two observations 
and often have a probable error as large as 25 km/sec. There is also an additional 
uncertainty in the radial velocities of the RR Lyrae variables caused by the 
uncertainty in the phase of the light variation. 

We have tentatively identified the motions of a few of the field RR Lyrae 
variables, that have large and relatively well determined motions, with those 
of subdwarfs. ‘These stars are listed in Table [X where various determinations 
of the proper motion and radial velocity and the space motions of the variables, 
for two values of the parallax are given. The mean absolute visual magnitude 
of the variables, determined by the group parallax method, together with RR Lyrae 
from the Groombridge 1830 group, is M,=+o0™-6. This value must be 
considered as preliminary since accurate photometry is not yet available for 
most of the variables. This preliminary value of +0™-6 should be compared 
with the usually adopted value of o™-o determined solely from a statistical 
discussion of the motions of these stars. It should also be noted that the most 
recent determination based on these statistics is +0™-5 (Paviovskaia 1953). 
But the agreement of +0™-6 and +0™-5 obtained by the two methods should 
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not be stressed since the inherent inaccuracies of the traditional statistical parallax 
method make it inferior to the method demonstrated here. 


Taste VIII 


Possible groups 


Wilson No. 7 U+U(O) V+V(O) W+W(o) V B-V_ Sp. My, 


“4 m m bi) 
6294 o'oI9 — 160 — 135 ° 8-24 o41 sdF2 +4°6 
7708 0°029 — 157 — 125 10 8-38 o53 FSV +5°7 
5909 O’OI5 —159 — 135 2 OOS Nas +5°7: 
6728 O°OT05 — 160 —126 12°2: +7°3: 


2302 0°050 + 38 — 190 8-52 0°88 
4005 0003 + 38 — 184 8-08 


+7°0 
+0°5 


HD 5916 0:0030 + 92 — 138 6°84 
8249 070055 + gI — 132 6°20 
6120 0°0055 + 99 — 128 


—o8 
o'r 
+03 


5730 O’o12 — 279 
5784  0°0017 —272 


+5°° 
—o'2 


[5058 O°ols — 198 +7°5 3] 


$236 0°0075 — 204 —177 


+4°4 
[6880 O°OI0s — 203 —174 


+6°0:} 


7O17 0°020 — 71 — 133 


+7°7: 
8296 o’ors — 74 — 145 


+5°6: 


10287 0°006 +160 — 285 


+0'9 
[10367 0009 +143 — 303 


+5°8:] 


3 
sdGo 
Ki V 
G2 III* 
t 
+ 
Tt 
sdF6 
t 
..  sdG2 
7464 o°’or2 — 195 S ° sdF4 -+5°1 
sdF3 
sdF7 
sdGo 
sdF7 
sdAz§ 
sdF8 


* 4005. Roman (1955): “ All lines are very weak’’. 

+ HD 5916 and Wilson Nos. 8249 and 6120. HD 5916 has been classified as G8 I1I-IV 
by Roman (1955) and as dG2 at Mt Wilson (Wilson 1953); the observed colour is much 
too red for either classification. No. 8249 (HD 122563) has been classified Go VI by 
Roman and sdFsp by Wilson and Joy (1950) who state that the hydrogen and metallic 
lines are weak; the colour is much too red for either classification. No. 6120 (HD 81192) 
has been classified as G8 II by Roman, who states, “‘ Strong CH and very weak lines ’’, 
and as G7 III+ by Keenan and Keller (1953), who state, “‘ CN and metallic lines weak 
and 4172 slightly strong”’. ‘The parallaxes used in computing the space motions of all 
three stars were obtained by fitting the stars to the (My, B— V)-diagram of Fig. 3 near 
the bifurcation of the giant sequence. 

t 5784. This star (HD 74462) has been classified as Gs IV by Roman with the 
statement, ‘‘ very weak lines and strong CH”’, and as sdGo by Wilson and Joy. The 
colour is too red for either classification and the space motion has been computed from 
a parallax derived on the assumption that the star is a globular cluster giant. 

§ 10287. This high velocity A-type star (HD 161817) is probably a horizontal branch 
star. The spectra are described in detail by Slettebak (1952). 





F 
t 
& 
§ 
x 
“S 
5 
4 
A 
£ 


gt 4 


7.94 


ozi4 $6 4 910-0 


2000.9 
6000.0 


Lzi 4 gGzi-4 
ori 4 $6 + 


£10.90 


6100.0 


0700.0 


$19.0 


£190.0 
4190.0 


(OVA+A (oun+n 


(9 


(¥Z) Mw 


(Z) PIN 


(%) M 


eS 


‘998/UIy OT — JO aNnTeA 2anBUsaITE US Buruunese fq 
Parwseni]s 6 AI°90[94 [wIpEs 944) UO 8101994 UOTjOUL azeds ay; Jo aoUapUsdap ay) “wuEsBOI}909d8 YOST AA 1)4 2UO UOdn spusdap Ayrojaa jeipes ayy 4 
PIECOCIY = PIX PUR LOST II = WM Sade > ‘Bunsdezis9y{ ~ ff ae Am A S861 “6 ‘sami agra =I, 


ztv1 4 4 Gh0.9 4 


£.qz + 910.0 - 


zZEOL Mm 


9 ttt0.0~ 





+ 400.0~— 
+ 020.0— 
I V1I9.0— 


d x 


sfapmpaqns sof nou, yjim pasodmuay saqiisne ashy py maj v sof suonou ands 


Xi iv] 


£1 J o70.0- 
6 TLto.0 


£1 199.0 ADM 


LlLo.-o- tvLigi aH 


+ ¥R0.0 
t$0.0 


ORO.0 
160.0 


"4 1S 


goly_Li+ 


FA LL 


IBC 





No. 3, 1959 Stellar groups; IV 277 


\ithough the accuracy of the group parallax value of M, < +06 should also 
not be emphasized at the present time, the potential accuracy of this method 
is So high that its exploitation promises to provide a definitive calibration of the 
luminosity of RR Lyrae variables when good astrometric and photometric data 


are available. 


Reval Greenwich Odservatery, Mt Wilson and Palomar Observatories, 
Herstmoncenx Castie, Carnegie Institution of Washington, 
Sussex : Caltiforma Institute of Technology, 
1959 Febreary 27. Pasadena, Cal., U.S.A, 
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ON THE EXISTENCE OF SUBDWARFS IN THE 
(M,,» LOG T,)-DIAGRAM 


Allan R. Sandage and Olin 7. Eggen 


(Communicated by the Astronomer Royal) 
(Received 1959 March 4) 


Summary 


The effect of weakening the Fraunhofer lines on. the observed colour 
indices of stars is examined. The blocking of the continuum radiation by the 
lines and the change in the temperature distribution of the solar atmosphere 
due to the lines is taken into account. If all of the lines were removed from 
the Sun, the change in the observed colour indices is predicted to be 
A(U — B)=o0™-32 and A(B—V)=o0"-17. These changes define a “ blan- 
keting line ’’ in the (U— B, B—V)-diagram from which the correction to the 
observed B—V colour due to line weakening can be found from the observed 
ultra-violet excess §(U — B). 

The magnitude and colour data for subdwarfs lying more than o™-25 below 
the normal main sequence are given. Colour corrections derived from the 
observed ultra-violet excesses are applied to the observed B—V colours. 
The corrections are sufficient to move the mild subdwarfs on to the Hyades 
main sequence. The conclusion is that most subdwarfs are below the main 
sequence in the (M,, B—V)-diagram only because of the effect of the weak 
Fraunhofer lines on their observed colours. Members of the Groombridge 
1830 moving group are not completely corrected to the main sequence and it 
is possible that these stars do in fact lie below the main sequence in the (M z,,, 
log T,)-diagram. ‘Tentative conclusions from the data are (1) the internal 
opacity of subdwarfs is probably due to free-free transitions of H and He 
rather than the photoelectric opacity of the heavy elements, and (2) the H/He 
ratio in the Groombridge 1830 stars may differ slightly from that of the 
Hyades. 

The blanketing model gives a relation between the magnitude difference of 
any star from the Hyades main sequence and the observed ultra-violet excess. 
The relation is applied to stars in different galactic clusters and to stars in the 
general field with the result that (1) the method of fitting main sequences to 
find distance moduli of galactic clusters is in error if corrections for the 
different ultra-violet excesses are not applied, and (2) the observed range of 
0°10 mag. in the ultra-violet excesses of field stars which are usually called 
‘““main sequence dwarfs ’’ shows that the total intrinsic width of the main 
sequence defined by these stars is at least o™-60. Because the strong-line field 
stars, the subdwarfs, and stars in different galactic clusters form a continuum 
in the (B—V), (U—B)-diagram, it seems probable that (a) the subdwarfs 
form a continuum in the (M,, B— V)-diagram, merging into the more normal, 
strong-line stars and (4) differences in the heavy element abundance exist 
between individual field stars and between stars in different galactic clusters. 


Introduction.—The presence of stars between the normal main sequence and 
the white dwarfs in the HR diagram has been known since the work of Adams and 
Joy (1922) and of Adams, Joy, Humason, and Brayton (1935). Modern 
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trigonometric parallaxes for the nearby stars leave no doubt that there are stars 
which lie as much as 1™-5 below the normal main sequence. However, because 
subdwarfs exist in the (My, spectral type)- and the (/,, colour index)-diagrams 
does not guarantee that they exist in the (M,,,, log T,)-diagram. Different bolo- 
metric corrections and a different colour-effective temperature relation may apply 
to subdwarfs and main-sequence stars. The present paper is concerned with 
the observed ultra-violet excess for known subdwarfs with well determined 
absolute magnitudes. Following Strémgren (unpublished) and Schwarzschild, 
Searle and Howard (1955), this excess is interpreted as an effect of the Fraunhofer 
lines on the continuous spectrum. A “‘blanketing model’’* is given using the 
blanketing coefficients for the Sun by Michard (1950). This model predicts 
that the change in B— V for a given observed ultra-violet excess, 5(U— B), for F 
and G subdwarfs is just sufficient to move most of the subdwarfs with known M, 
on to the Hyades main sequence. 

The observational data on subdwarfs with well determined absolute magni- 
tudes are given in Section 1, where the ultra-violet excess is discussed. The semi- 
empirical blanketing model giving the A(B— V) =f[8(U — B)] relation is derived 
in Section 2. Here the predictions of the model are compared with the obser- 
vations. The model is applied to the observed B—V values and the corrected 
positions of the subdwarfs in the (M,, B— V)-diagram are shown. To interpret 
the observational results in terms of stellar structure, the well-known theoretical 
equations are then reviewed. Other applications of the theory are given in 
Section 3. 

1. Observational data for subdwarfs.—The term subdwarf was first used by 
Kuiper in 1939 to describe stars lying between the main sequence and the white 
dwarfs. Subdwarfs had been found several years earlier in the Mt Wilson 
spectroscopic parallax programme when a list of 6 so-called ‘‘ intermediate white 
dwarfs’’ was published (Adams, Joy, Humason and Brayton 1935). It is 
important to point out that most ‘‘ subdwarfs ’’ have been so termed not by definite 
knowledge of their position in the HR-diagram but rather by the peculiar spectral 
features or by large proper motions. For example, in the searches carried out by 
Kuiper from 1939 to 1942, the parallaxes of proper motion stars were inferred 
from the assumption that the tangential velocities were less than 474 km/sec, that 
is, that 7/u 0-01. This procedure, and its justification, is described by Kuiper 
in the publication of the Paris Colloquium on white dwarfs (1941). Clearly the 
method can give only very approximate estimates of the real absolute magnitudes 
of these stars. Kuiper’s early estimates placed the subdwarfs in a band extending 
to 3 magnitudes fainter than the main sequence in the HR-diagram (Kuiper 1941, 
Fig. 1). 

P. P. Parenago later discussed the data available for subluminous stars in 1946. 
He used absolute magnitudes derived entirely from trigonometric parallaxes 
(Parenago 1946, Table 2) but only eight stars, other than white dwarfs, had 
a2>0"-050. Parallaxes smaller than this are usually so uncertain as to be of little 
individual value. The HR-diagram determined by Parenago shows a band of 
subdwarfs lying 2™-5 below the main sequence. 

In recent years it has been pointed out that the subdwarfs are usually assigned 
an earlier spectral class than that compatible with their effective temperature. 


* The nomenclature will be described in Section 2. 
21* 
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This misclassification, which is partly due to the weakness of the spectral lines, 
has the effect of placing the subdwarfs too far to the left in the HR-diagram and 
therefore introducing a false displacement from the normal main sequence. It 
has been pointed out by Greenstein (1956) that the colour-index is a better indi- 
cation of the true temperature of subdwarfs than the spectral type. It has been 
shown, however, that the colour is also affected by the weakness of the spectral 
lines (Schwarzschild, Searle and Howard 1955). A fundamental question is 
whether a subdwarf sequence remains in the (M,, B—V)-diagram when proper 
correction is made for the effect of the weak Fraunhofer lines on the observed 
colour. 

To investigate this question we have selected all the stars, later than Fo, which 
fall more than o™-25 below a standard main sequence (Eggen 1955 6, Table III). 
Only stars with photoelectrically determined magnitudes and colours and 


TABLE I 
Subdwarfs more than o™-25 below the main sequence 


HD 7 P-—-V B-V U-B Sp m+P.E. Wt. My 
m Unit 0-001 m 
55 0°94 wii ae 72+ 780+: 
6582 o-59 «(«O'71 o11 GoVI 136+ 581+: 
10700 062 0°72 020 G8Vp 275+ 5°7O+ ° 
14412 062 o'72 o17 GsV 88+ 6-03+ ° 
17925 ‘ 078 = 8600°87 0°56 dKo 122+ 6-46+ - 
20766 =: o'52 0°64 006 ~ G2V 106+ $°3t+ ° 
20807 5° 0°46 058 —o'04 G1IV 97+ 5‘o2+- 
25329 0-78 «60°88 0°36 KiV 54+ 718+ - 
ae 080 «089 065 K2v 56+ 7°06+ - 
38230 . 0°73 0°83 dK2 77+ 669+ - 
64379 0°34 0°44 —0°06 F5V 67+ 4°11t° SDS, 3° 
65583 o°59 G8V 56+ ae: 
74377. (8: 084... .. dKs 59+ 7°38t ° 
89125 5° 0°39 “ —o-o5 dF3 60+ 4°70+ ° 39 Leo. ADS 
7712A 
6-72+ ° Gmb 1830 
5°04+ °23. 10CVn 
6°47+ ° 
7°32+ ° 
519+ ° 
3°48+ ° 
3°87+ ° 
6-19+ - 
729+ ° 
655+ ° 
691+ ° 
os a : KoV 79+ 6-95+ ° 
0°52 062 —0'06 sdF8 group... 60 +° Member 1830 Gr. 
o-76—_(i.. ves, 171+ 8 649+ ° SDS, HR 7703 
0748 058 0106 Gs5V 66+ 5 504+ - 
a6: @%  .... ee) 111+ § 4°49+° y Pav 
17 0°36 047 —0o-20 sdF8 group... 48 +: Member 1830 Gr. 
76 055 067 4«®0'05 G2V 86+ 4 69 5°43+°- 85 Peg. ADS 
17175A 


UuUuuUn usu NN AOU DO 


103095 
110897 
112758 
114606 
126053 
128167 
134083 
144579 
144872 
145417 
151877 
154577 
163810 
191408 
193664 
203608 
219617 


0°64 o17 sdG8 116+ 
0-42 054 —o'or GoV 65+ 
a, ee osu a 61+ 
o'52 O°! o'o5s G3V 52+ 
0°52 0°63 008 dG3 61+ 
0°24 O36 —o-08 F2V 63+ 
O31 o42 —0°03 F5V 61+ 
0°63 0°73 o:21 dG8 80+ 
082 096 0-71 K3V 56+ 
0°70 : KoV 63+ 


c+ 
OO 


NUUMswus ADD Sw 


0-72 o82 044 K7V 50+ 


MW Op UMN Os Ds DAH HOM DH 
j CH AnreuunQaset Sslw 
-+#NNWNS DAW ADS 


* ADS 10938. Am=1™. 

+ ADS 16644. The magnitude of one component is assumed to be My= 5-55 because Am=o for the 
binary. 

t85 Peg. Not plotted in Figs. 1, 2, 3 or 4 because the photometric values refer to both components. 
Am= 2°93. 
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zt 2 0"-05 based on observations from more than one observatory have been used. 
The available data for these stars, and for the members of the Groombridge 1830 
group discussed elsewhere (Eggen and Sandage ig59) are listedin Table I. The 
HD and the parallax catalogue (Jenkins 1952) numbers identify the stars. The 
values of V and P— V are taken from a previous publication (Eggen 1955 a); the 
values of U— B and B— V are taken from the photometric catalogues by Roman 
(1955), Johnson and Morgan (1953) and Johnson and Knuckles (1957). We 
have also determined U — B and B-— V for three southern stars with the 18-inch 
reflector at the Cape Observatory. 














O-2 
B-V. 


Fic. 1.—The (My, B—V)-diagram for stars with n>o"-05 determined at more than one 
observatory which lie more than 0™-25 fainter than the normal main sequence. The three crossed 
circles are members of the Groombridge 1830 group. The gap between the points and the main 
sequence is artificial and results from the requirement that AMy >o™-25. The open circle indicates 
Yale No. 3300. The uncertainty in My due to the probable error of the parallaxes is indicated by 


vertical lines. 


It might be argued that the sample of subdwarfs in Table I is too small to be 
representative, especially since only five stars are common with a previous tabu- 
lation (Eggen and Sandage 1959, Table V) of subdwarfs discovered by their 
spectral characteristics. However, the previous tabulation includes the extreme 
subdwarf members of the Groombridge 1830 group, which are also included in 
Table I, and since the departure of the group members from the main sequence is 
nearly matched by other stars in Table I, it is thought that, except for one doubtful 
case discussed below, the data in Table I represent the total spread of subdwarfs 
in the (M,, B—V)-diagram. As usual, the individual trigonometric 
must be considered with some caution, as can be seen from the separate deter- 
minations listed in Table II; note particularly the disagreement for Yale 
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Nos. 1002, 2974, 3300. Yale No. 3300 is indicated in Figs. 2 and 6 as an open circle. 
The stars in Table I provide the only absolute magnitudes of sufficient precision 
to define the subdwarf region in the (M,, B— V)-diagram. 


Taste II 
Individual parallaxes for stars in Table I 
(Unit 0”:oo01) 


Te 
3 89(Y), 53(C) 
143(A), 138(M), 105(Yk) 
365 312(M), 277(Y), 257(S), 288(C) 
479 115(Y), 76(C) 
599 133(M), 113(Y), 131(C) 
7O1 107(C), 104(C) 
705 98(C), 63(C) 
887  45(M), 39(Yk), 72(S), 54(V) 
1002 106(M), 30(V) 
1314  86(A), 75(M), 54(W) 
1864 71(Y), 68(C) 
1889  42(A), 73(M), 55(Yk), 61(W), 66(V) 
2082 51(A), 45(M), 49(W), 104(S) 
2413 56(A), 56(M) 
2745 108(A), 106(M), 164(St) 
2935 58(A), 65(M) 
2974 78(M), 37(C) 
3010 46(A), 53(M) 
3264  61(A), 65(M), 44(W) 
3300 ©=_- 70(A), 28(M) 
3416 54(A), 67(M) 
3650 = 77(A), 78(M) 
3654  50(A), 63(Yk) 
3669 «= 71(Y), 53(W), 75(C) 
3827 42(A), 67(M) 
3884 56(Y), r19(C) 
4782 177(Y), 167(C) 
4849 = 50fA), 85(M), 71(G), 66(S) 
5152  124(Y), 98(C) 
5807  84(A), 57(M), 77(Yk), 127(W), 114(S), 63(St) 


Adopted = 91 


Fig. 1 shows the (/,, B — V)-diagrams defined by the stars in Table I that are 
bluer than B — V = o"'-80 and for which values of U— B are available. The main 
sequence of the Hyades (Heckmann and Johnson 1956) is shown for comparison, 
and members of the Groombridge 1830 group are shown as crossed circles. One 
star, Yale No. 3010, falls considerably outside the region of the (M,, B—V)- 
diagram populated by the other starsin Table I. The spectroscopic observations 
of this star, discussed previously (Eggen 1955 5), indicate that the spectrum is no 
more abnormal than the intermediate subdwarf 4 Cassiopeiae and certainly not as 
abnormal as the extreme subdwarf Groombridge 1830. Although there are other 
stars, whose absolute magnitudes depend upon only one determination of the 
trigonometric parallax, that also fall in the same region of the (M,, B— V)-diagram 
as Yale No. 3010, we will base our arguments on the remaining stars in, Table I 
with well determined parallaxes. It is important, however, that the parallax of 
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Yale 3010 be checked and that a second determination be made for such stars as 
Yale Nos. 291, 595, 1282 and 4148. 

Many of the stars in Table I show an ultra-violet excess with respect to a 
standard (U—B, B—V)-relation. The (U—B, B—V)-diagram is shown in 
Fig. 2 where the crossed circles are stars in the Groor ’ idge 1830 group and the 
broken curve, parallel to the adopted, normal relation, 13 the position of the (U —B, 
B—V)-relation for giants and subgiants in M3 (Johnson and Sandage 1956). 
Note that the Groombridge 1830 stars follow the globular cluster relation fairly 
closely, a point discussed in a previous paper (Eggen and Sandage 1959) when 
fitting the globular cluster main sequence to that of the Groombridge 1830 group. 














 * 
Fic. 2.—The (U—B, B—V)-diagram for the stars in Table 1. The lower line is defined by 


Hyades stars. The upper line is defined by the giant and subgiant stars in M3. The same symbols 
as those in Fig. 1 are used. 


The standard (U—B, B—V)-relation in Fig. 2 is that defined by Hyades 
main sequence stars (Johnson and Knuckles 1955, Heckmann and Johnson 1956) 
which are probably unaffected by companions. Because the (U-B, B—V)- 
relation is known to vary from cluster to cluster, we will use both this relation and 
the (M,, B—V)-relation for Hyades members alone in what follows. This is 
important because the theory in Section 2 gives AM, from a main sequence 
defined by certain stars as a function of the ultra-violet excess with respect to the 
same stars. : 

We have not used the intrinsic (U — B, B— V)-relation given by Johnson and 
Morgan (1953, Table 14). Their mean values are nearly o™-04 brighter in U—B 
than a sample of nearby stars from B— V = 0™-47 to o™-60 (Johnson and Knuckles 
1957). _Wesuspect that some of the stars used by these authors in forming means 
in this range of B—V are binary stars, evolved subgiants, or subdwarfs which 
should not have been used if a pure main sequenc: relation is desired. In any 
case, because the intrinsic relations for the Hyades, Praesepe, and Coma Berenices 
are known to differ from each other (Johnson and Knuckles 1955) the gross mean 
of all field stars is expected to have a scatter of at least + o™-04 at a given B—V, 
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‘The differences can be interpreted on our present model if the metal abundance is 
lower in Coma Berenices stars than in Hyades members. ‘This point is discussed 
in Section 3. 

2. The blanketing model.—The essential point of the theory presented here is 
that the Fraunhofer lines affect the U, B and V regions of the spectrum in different 
ways so that a weakening of the lines produces changes in the observed colour 
indices, U—- Band B—V. If the relationship between the effect on U—B and 
on B— V is known, then the correction to the observed B—V can be computed 
from the observed ultra-violet excess. Because the observed B—V for weak 
line stars will be bluer than that for strong line stars of the same temperature, the 
weak line stars will fall below the standard main sequence. Therefore, because 
of the relationship between A(U— B) and A(B—V) we should expect that the 
displacement of a weak line star below the standard main sequence will be 
correlated with the observed ultra-violet excess. 

The Fraunhofer lines affect the observed energy distribution of a star in at 
least three ways. (1) The lines block the radiation of the continuum and at a 
given wavelength the effect is to reduce the observed emergent flux. (2) The 
blocked radiation affects the temperature gradient of the atmosphere itself and 
therefore changes the continuum flux between the lines. (3) A change in the 
metal abundance could change the opacity with a resulting change in the emergent 
flux and the size of the Balmer jump. 

We now consider each of these effects quantitatively. 

(1) Let (A) be the fraction of the energy (in terms of the continuum) absorbed 
by lines at A. Let S(A) be the sensitivity function of the particular photometric 
system (telescope, filters and photocell) and let F(A) be the continuum flux of the 
source between the lines. When all lines are removed, the change in the magni- 
tude of the source in a spectral band from A, to A, due to the first effect is 


Amag=2'5log| | * rsdn] |" FS(~e)aa], (1) 


Unfortunately, adequate values of ¢(A) do not exist for stars other than the Sun 
in the wavelength interval of the (U, B, V)-system. For a first approximation 
of this effect we are forced to use the Sun, for which values of ¢(A) from A = 3100 to 
A= 8700 are available (Michard 1950). Extensive blanketing measurements in 
stars of different spectral types are now in progress by the Burbidges, Sandage 
and Wildey, but the results are not yet available. A major difficulty in these 
measurements is in assigning a precise value of the continuum level in the shorter 
wavelengths of the spectra. 

Michard’s values for «(A), determined from a study of the continuum level in 
the Utrecht Atlas, were used. They refer to the centre of the disk. The S(A) 
values are taken from Johnson’s definition of the (U, B, V)-system (1955). The 
values of F(A) were taken from Minnaert (1953, Table 2) and refer to the con- 
tinuum between the lines at the centre of disk. ‘The values of «, S and F substi- 


tuted in equation (1) show that if all lines are removed from the Sun the changes 
in U, B and V due to the first effect would be 


AU= 0-653 

AB=0:282 

AV =0°075 (2) 
A(B-—V)=o0-21 


A(U — B) =0°37 
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(2) The most recent study of the change in the temperature gradient of the 
solar atmosphere due to the Fraunhofer lines is that by A. Przybylski (1957) 
where a historical summary with references should be consulted. The effect was 
first studied by Milne (1928) and was labelled ‘‘ blanketing ’’ because the lines 
act as a blanket on the solar atmosphere causing the continuum temperature to 
rise over what it would have been in the absence of lines. (In the present paper 
we use the term “‘ blanketing ’’ to denote the total effect of the lines on the emergent 
flux observed with wide filter bands, that is, the sum of our three effects.) Using 
Pecker’s computations (1951) of the line absorptions for various optical depths, 
Przybylski shows that in the presence of lines, the temperature at the mean depth 
for continuum formation (7 0-6) is higher by about 200°K than it would be if 
the lines were removed. Hence, on removing all the lines, the continuum colours 
should be redder than those predicted in (2). Approximating the 200° change 
by a change in the emergent flux of a black body gives the correction due to this 
second effect as 


A(B- V)= ros6CA(1/T)| 5 - x|- aes a4 , 


where A, and A, are the effective wavelengths of the B and V kaisdiaiinan and C, 
is the second radiation constant. Substitution of T’=5700°K, AT=200°K, 
A, = 4430 and A, = 5540 gives A(B—V)=0™-04. The corresponding A(U — B) is 
estimated from the slope of the U- B=f[(B-—V)] relation (Table III) near 
B—V=0"™-6 to be A(U— B)=0™-05. 

(3) Subdwarfs are known to have low abundances of the heavy elements 
(Chamberlain and Aller 1951, Aller 1958). The lowering of the metal abundance 


is expected to change the electron pressure at a given temperature in Go stars with 
a resulting change in the atmospheric opacity. However, computations by 
Kandell (unpublished) show that the change has no effect on the emergent conti- 
nuum flux. A computation of the expected change in the colour temperature— 
eifective temperature relation as a function of electron pressure gives the same 
result as is shown by comparison of Figs. 1 and 3 of Sandage (1956). Fig. 3 
shows that a decrease in the metal abundance will decrease P, at given 0, and log g 
values if @,>0-8. But Fig. 1 of the same reference shows that the (@,, 0.) relation 
is insensitive to P, for 8, >0-8, which proves the result. Consequently no effect 
on the colour is expected from the change in P, by decreasing the metal abundance. 

In summary then, the total change in the colour of the centre of the solar disk 
from these three effects is 


A(U— B) =0™-37 —0™-05 =0™-32 
A(B— V)=0™-21 —0™-04 =0™17. (4) 
The observed colours for the Sun are B—V=0"-63 and U—B=o0-08 
(Stebbins and Kron 1957). These values refer to the integrated light of the entire 
solar disk. The B — V value agrees well with that expected foradwarf G2. Because 
the centre of the solar disk has a spectral type of dGo, we expect its colour to be 
B—V=o0™-s5gand U— B=0™-04; values obtained by interpolation in the standard 
spectral type-colour relation. These values indicate that the Sun has an ultra- 
violet excess of about o™-o8 with respect to Hyades stars and, on the present model, 
this suggests that the Sun has weaker Fraunhofer lines than the Hyades. With 
the ratio A(U — B)/A(B — V) =0-32/0-17 from equations (4) and with A(U— B) = 
+o™-08 the predicted A(B— V) for the Sun is o™-04. Consequently, if all lines 
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were removed from the stronger-line Hyades stars, the predicted colourchanges are 
A(B— V)=0™-17 + 0™-04 =0™-21 and A(U— B)=0™-32 +0™-08 =0™-40. These 
changes define what might be called the ‘‘ blanketing line’’ for Go Hyades stars 
inthe (U—B, B—V)-diagram. ‘The line is shown in Fig. 3 and is the path which 
a Go star in the Hyades would follow if we progressively weakened the Fraunhofer 
lines and if the weakening proceeds linearly for all wavelengths. When all lines 
are removed, the observed colours should be B— V =0™-60 —0™-21 = 0™-3g and 
U— B=0™-13 —0™-40= — o™-27. 

For the present investigation we need the blanketing line for a range of spectral 
types from about Fo to Ko. These can be derived in the same manner once 
(A) values are available for a variety of spectral types. However, until work 
now in progress is completed, we must use the Go line and assume that the error 
will not be too large for small differences in spectral type. 





INTRINSIC LINE. 


os 











° 


Fic. 3.—The blanketing line for dGo stars is shown starting at B—-V=+0™-60, U-B= 
+oM-13 and ending at B—V = +0-39, U-B=—o™-27. The end point is the predicted position 
of a dGo Hyades star if all the Fraunhofer lines were removed. The two stars analysed by 
Chamberlain and Aller are shown as open circles. 


Let 6(U — B) be the observed ultra-violet excess [Fig. 4 illustrates the defi- 
nition of A(U—B) and 6(U—B).] The values of A(B— V) corresponding to a 
given 5(U — B) can now be determined from Fig. 3. As an example, suppose a 
dGo star of normal colour, B— V =o™-60, is made bluer, A(B— V)=0™-10, by 
slightly weakening the lines. ‘The blanketing line in Fig. 3 predicts that the 
observed U— B colour will be — o™-06 at B— V=o0™-50. According to Table III, 
this observed colour corresponds to an excess of 6(U—B)=o0™-o0g. Con- 
sequently, on this model any star with an observed excess of 5(U — B)=0™-o9 
should have the correction A(B—V)=0™-10 applied to its observed B—V.,to 
reduce it to the colour system of strong line stars (Hyades). The predicted 
values of A(B— V)=/[8(U — B)], derived in this manner, are listed in Table IV. 
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Table IV also lists the values of AM,, the displacement below the Hyades main 
sequence, due to the effect of Fraunhofer lines on the colour alone. These values 
of AM, are derived from the blanketing line which a dGo star will follow in the 
(M,, B— V)-diagram as the lines are weakened. This line startsa B— V =o™-60, 
M,= +4™-78 on the Hyades main sequence and ends at B—-V =o™-39, M,= 
+4™-70, when all the lines are removed. Combining the A(B— V)=/[8(U— B)] 
relation previously found with this AM, =f[A(B — V)] relation, gives the required 
theoretical relation between distance below the main sequence and the observed 
ultra-violet excess. 














Fic. 4.—The relation between observed ultra-violet excess, 8(U—B) and the blanketing corrections 
4(U—B) and 4(B—V). 


Taste III 


Adopted main sequence and two colour relation 


B-V My U-B B-V My 
m m 


m m m 

0-20 re 0°13 0"70 5°32 
0°25 oe oll o"75 5°60 
0°30 Pa 0°07 0°80 5°85 
0°35 3°10 0°03 0°85 6°10 
0°40 3°42 ool "90 6°34 
0°45 3°78 0°00 0°95 6°38 
0°50 4°10 0°03 1°00 6°80 
O55 4°45 0-08 

o'60 4°78 o"13 1°05 7°02 
0°65 5°02 o'19 1‘ro 725 099 


A check on the procedure to this point is highly desirable. The first test 
concerns the extreme value which 6(U—B) can take when all the lines are 
removed from a dGo star. Table IV shows that this value is predicted to be 
5(U — B) =o™-28 for the observed colours of B— V =o™-39 and U— B= —o™-27 
(see Fig. 3). A plot of all stars of large ultra-violet excess in Roman’s (1955) list 
of high velocity stars shows that stars at B—~V=0™-39 indeed are not bluer 
than U— B= —o™-27 and this upper limit is well defined observationally. 

The second test is concerned with the two subdwarfs HD 19445 and 
HD 140283 which Chamberlain and Aller (1951) have analysed. These stars, 
previously classed as A4p and Asp, are plotted in Fig. 3 with their observed colour 
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indices. If the blanketing model is correct, when these stars are moved back along 
the blanketing line they will fall on the (U— B, B— V)-diagram at a value of B— V 
which corresponds to the effective temperature predicted by Chamberlain and 
Aller. The results are shown below where the tabulated 5(U— B) values were 
read from Fig. 2 and the corresponding A(B— V) values taken from Table IV. 


HD (B-V),, (U-B),, 8(U-B) A(B-V) (B-V)pg,, 
19445 +0™-46 —o™-24 o™-21 o™-18 +0™-64 
140283 +o™-48 —oM-10 o™-08 o™-08 +o™-56 
Fom6o 
‘The mean of the colours that the stars would have if the lines were not weakened 
is (B—V).,,,= +0™-60 which corresponds to that of a dGo star for which 
T,=6000°. This is in good agreement with Chamberlain and Aller’s pre- 


dicted value of 6300° and differs considerably from the temperature of A-type 
stars. 


TABLE IV 


Theoretical \(B — V) =f [(8(U — B)| and 3(U — B)=f(AMy) relations 
for dGo Hyades stars 


3(U — B) A(B-V) AMy 
m m m 
0°00 0°00 0°00 
0°02 0°02 O°14 
0°035 0°04 0°26 
0°05 0°06 0°37 
0°075 0°08 o's 
0°09 o"10 0°63 
o'r! or12 o"75 
0°13 O'14 og! 
O17 016 1°02 
o'21 o18 I°I4 
0°26 0°20 1°26 


A similar conclusion is available from Bahng’s (1958) six-colour results and 
interpretation of Groombridge 1830 itself. _Bahng shows that the red and infra- 
red colours of Groombridge 1830 correspond to a temperature of a dK3 star but 
the observed B—V of o™-75 corresponds to an earlier spectral type; a difference 
which is attributed to the effect of the lines. The ultra-violet excess of Groom- 
bridge 1830 is 6(U—B)=o™15. For this 6(U-—B) Table IV shows that 
A(B-—V)=o™15, giving a corrected colour of (B—V),,,,=+0™go. This 
value corresponds to a dK2 star, in good agreement with Bahng’s estimated 
““temperature type’’ of dK2-8. Agreement between the observed and pre- 
dicted values in both of these tests of the blanketing model gives some confidence 
in the AM, =f[5(U — B)] relation tabulated in Table IV. 

Using Table IV and the observed 5(U — B) values read from Fig. 2, we now 
correct the observed (B— V’) for the stars in Table I earlier than Ko. The results 
are given in Table V where successive columns contain (1) the Yale No., (2) the 
observed M, copied from Table I, (3) the observed displacement below the 
Hyades main sequence, read from Fig. 1 at constant (B—V), (4) the observed 
ultra-violet excess, (5) corrections to the observed B—V taken from Table IV, 








No. 3, 1959 Subdwarfs in the (M,.,, log T,)-diagram 289 


(6) the observed B—V corrected for the line weakening and (7) added identi- 
fications for the stars. The correlation between AM, and 5(U—B) in the third 
and fourth columns of Table V is shown in Fig. 5 together with the predicted 
correlation given in Table IV. The agreement between the observed and pre- 
dicted correlation in Fig. 5 is fairly good, especially when the uncertainties in 


TABLE V 
Corrections to B—V for stars in Table I earlier than Ko 


(1) (2) (3) (4) (5) (6) 
Yale No. My AMyt &U-B) A(B-V) (B-Vcon 
™m m 

219 «= § 81 + “10 0°42 +°10 o'16 0°16 0°87 p Cass 
365 5°70+°04 0°28 + °04 0°08 "09 o°81 7 Cet 
479 6:03+°17 0°60 + °17 o'l2 o°13 0°85 
7O1 5°31 4°25 o'41+°25 ol oz 0-76 ¢1 Ret 
705 5°02 +°27 0°40 +°27 ? o"14 o'72 ¢? Ret 
1864] 4:11+°23 0°40 +23 . 0°07 o°51 SDS, 3” binary 
1889 §= 55°75 £"19 0°44 4°19 ; 0-09 o"79 

2413* 4°70+'18 0°60 +18 ° 0-09 o'59 39 Leo ADS 

7712A 

2745 6°72+°09 1°17 +09 ols 0°89 Gmb 1830 

2935 5°04+°23 0°65 +°23 0°07 0°08 0°62 10 CVn 

3010 7°322°25 2°51 +°25 o"09 oro o’7I 

3264 5"19+°21 0°27 +°21 0°09 o"10 o'7”73 

3300 «=. 3°48 “21 o'30+°21 oll or12 0°48 a Boo 

3416 83°87 +-25 0°28 +°25 0°03 0°03 0°45 45 Boo 

3650 «= 6 19+°14 0°69 +°14 0°09 o'10 083 

HD 163810 6:0 + ‘1 rizt ‘1 o-21 018 080 Member 1830 Gr 
4849 5°04+°'16 0°40 +°16 0°05 0°06 0°64 
HD 219617) 5555+ *1 165+ ‘1 o21 o18 0°65 Member 1830 Gr 


(7) 


t Observed AMy taken from Fig. 1. 
} Am=4 mag. Secondary can be neglected. 
* Am=5°5 mag. Secondary can be neglected. 


the absolute magnitudes of some of the stars is considered. This agreement is, 
perhaps, better seen in Fig. 6 where the (M,, B-- V)-diagram for all the stars is 
plotted after applying the theoretical corrections in Table IV. Omitting one 
star (Yale 3010), which we have already discussed, and the members of the 
Groombridge 1830 group which are indicated by crossed circles, no star deviates 
from the Hyades main sequence in Fig. 6 by more than three times the uncertainty 
in AM, due to the probable error of the parallax. Furthermore, if we allow for 
the possibility of errors in the observed values of U—B up to o™-03, all of the 
subdwarfs in Fig. 1, with the exceptions already noted, will fall within their pro- 
bable errors on the Hyades main sequence in Fig. 6. We conclude, therefore, 
that these subdwarfs do not in fact form a separate region in the (M,,, log T,)- 
diagram, even though it is probable that their heavy element content is lower 
than that of the Hyades. Tointerpret these results we now consider the equations 
from the theory of stellar interiors. 

It is well-known (Strémgren 1952, Massewitch 1954, Swarzschild 1958) 
that subdwarfs can, theoretically, be produced in the (M,,,, log T,)-diagram by 
at least two processes. (1) If the internal opacity in normal main sequence 
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Fic. 5.—The relations between observed ultra-violet excess, 8(U—B) and the distance below the 

Hyades main sequence due to the weakening of the Fraunhofer lines alone. The observational 

data are from Table V. The theoretical line is tabulated in Table IV. The Groombridge 1830 

members are indicated by crossed circles. 














Fic. 6.—Same as Fig. 1 but with the corrections to B—V applied to the individual stars according 
to the size of the observed ultra-violet excess. 


stars is primarily due to bound-free photoelectric transitions of the heavy elements, 
a decrease in the heavy element content will increase the luminosity, L, but will 
keep the radius, R, nearly constant, thereby placing the star below the main 
sequence (cf. Eqns. 9-12). If, on the other hand, the opacity of normal main 
sequence stars is due to free—free transitions of hydrogen and helium, no such 
effect is expected from changes in the heavy element abundance. (2) A decrease 
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of the hydrogen abundance with a corresponding increase in helium will also 
move a star below the main sequence. This is brought about by an increase in 
the mean molecular weight which causes L to increase and, because as above, R 
is nearly insensitive to the change, the star moves below the main sequence. 
Only properly computed models can show these effects precisely but the general 
outline can be demonstrated by the following simplified arguments. 

Consider first the case where the opacity is primarily due to bound—free tran- 
sitions of the heavy elements. Then 


«x =const Z(1 + X)pT,-* (g/t) (5) 
where g is the Gaunt and ¢ the guillotine factor. For any star in hydrostatic 
equilibrium with energy transport by radiation, and with the opacity given by (5), 

const wp *°M55R-%5 / ) 
L= -}. 6 
Z(1+X) \g (6) 
A second equation for L comes from the requirement that the energy generated 


by nuclear processes be equal to the observed luminosity. If the C—N cycle 
dominates, with the energy generation given by 





€oy =const ZXpT" ergs/gm sec (7) 
then 


const ZXpu"M"t2 
bs wr (8) 





L 


Equations (6) and (8) give two relations in 4 variables (L, M, R, composition). 
Equating (6) with (8), with n= 20, gives 
L=const Z-*%X-4*4>8 (9) 
R=const Z°X-0:16 40-78 (10) 
where the chemical parameters have been approximated by power series. The 
(g/t) factor is absorbed in the constant because we are interested only in differential 


effects at constant mass and nearly constant radius. From 7,=const (LR-*L)"* 
it follows that 
T,=const Z-0-*X-0- egos (11) 
Elimination of the mass from (g) and (11) gives 
L=const Z°@X1-87) 56 (12) 
which is the equation for the main sequence in the (L,,,, T,)-diagram. 

We must emphasize the approximate nature of these results. For one thing, 
convective energy transport and its effect on the surface conditions has been 
neglected (cf. Schwarzschild 1958). For another, if the structure of the star 
changes with changing X and Z (as it must to some degree) then the constants in 
equations (9g), (10), (11) and (12) will themselves be functions of X and Z and 
the argument fails. However, our approximation is probably adequate because 
there is a check on the equations from the model computations by Blackler (1958) 
for main sequence stars with different hydrogen contents. Miss Blackler’s 
models were computed for X values of 0-76, 0°85, 0°93, and o-99. Equations 9 
and 11 predict the change in L and T, for constant mass, and constant Z, with X 
changing from 0-99 to 0°76, to be 


Log (L,/L,) = 4 log(X,/X,) =0°46 (a) 
Log (74;/T.2)= 0-92 log (X,/X,)=0-105 _ (0) 
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Miss Blackler’s computed values are 0-64 and 0-47 for (a) with masses of 1 and 
2Q, respectively and 0-203 and 0-126 for (6) with the same mass values. The 
agreement with our equations is satisfactory. 

Equation (12) illustrates the two methods described above for producing sub- 
dwarfs. Atconstant 7’,, either a decrease in Z or in X will decrease L, moving the 
star below the main sequence. 

It should be pointed out that a decrease of Z will not produce subdwarfs if 
equation (5) does not correctly represent the opacity because of the predominance 
of free-free transitions of hydrogen and helium. If this is true, then equation 
{5) should be replaced by 

«=const (X + Y)(1+X)pT,-**(g/t) (13) 
and the resulting equations for L, R and T, become 
L=const Z-°8X-4+)4*"18 (14) 
R=const Z2°%X-0-16 40-78 (15) 
T,=const Z-°8 X—092 40 (16) 
L=const Z°4X1-47' 56, (17) 
where L, R, and T, are nearly independent of Z. 

The above equations are derived on the assumption that the energy is gener- 
ated by the CN cycle because we wished to compare our results with Blackler’s 
model computations already mentioned. However, it is probable that the 
energy generation in most of the subdwarfs discussed here is by the proton- 
proton chain, for which n=4 in equation (7). Although the exponents in the 
equations are changed if n= 4, the sign of the dependency of L on chemical com- 
position in equations (12) and (17) is the same and therefore our general con- 
clusions are independent of whether the energy is generated by the CN cycle or 
the p-p reaction. 

The question is whether the opacity in real stars of solar mass is due primarily 
to bound-—free or to free—free transitions. Nauer’s models (1954) show that near 
the centre of the Sun most of the opacity comes from the free—free transitions 
while near the surface it is due to bound-free transitions. It is therefore not 
known, from these considerations, whether equation (12) or (17) will best relate 
T, and L for main sequence stars of solar mass. 

Because a change in the heavy element abundance must form subdwarfs if 
the opacity is due to bound-free photoelectric transitions (equations (5) and (12)), 
and because Fig. 6 shows that the difference in the heavy element content does 
not change the structure of the star as equation (12) predicts, it is suggested that 
equation (5) does not in fact represent the opacity correctly but rather equation 
(13) applies. This conclusion is not completely certain because equation (12) 
shows that some combination of increased Z with decreased X would cancel the 
effect on L of either difference working separately. If the combination was 
properly chosen (Z°*X1"~ const), no change in L at a given T, would result 
and Fig. 6 could be explained even if equation (12) were valid. Such a com- 
bination might be justified from theories of element synthesis together with a 
changing chemical composition of the interstellar medium with time; ideas 
which are now in vogue. In particular, if heavy elements are made in stellar 
interiors and scattered into the interstellar medium by stellar explosions, an 
increased Z abundance might be expected to be accompanied by increased He 
abundance. Thus, an increase in Z might logically go in step with decreased X. 
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On the other hand, if it could be independently shown that equation (5) does 
not control the opacity, but that equation (17) correctly describes the situation, 
then the conclusion from Fig. 6 would be definite that no difference in the H/He 
abundance exists between Hyades stars and stars in Table I. 

The case of the three Groombridge 1830 group stars is more interesting 
because the colour correction of the theory does not appear sufficient to move the 
stars to the main sequence. Before correction they average 1™3 below the 
Hyades main sequence. After correction they are still below the main sequence 
but only by o™-35 on the average. If this is significant, then equations (17) shows 
that the hydrogen content in Hyades stars is higher than that of stars in the Groom- 
bridge 1830 group. If, on the other hand, the opacity is controlled by equation 
(5), then the Z(Hyades)/(Z(Gr. 1830) ratio is approximately 2 by equation (12) if 
X is kept constant. However, Chamberlain and Aller (1951) and Kondo (1957) 
show that the heavy elements are down by at least a factor of 10 in stars like Groom- 
bridge 1830, a factor which is considerably greater than 2. Therefore, the 
present data suggest that a change in the heavy element content does not affect the 
stellar structure of the Groombridge 1830 Group members and that equations (13) 
and (17) probably represent the case better than (5) and (12). If this is true, then 
the data show that the H/He ratio must be different between the Hyades and the 
Groombridge 1830 group members. 

Clearly these conclusions are quite tentative because the accuracy of the 
blanketing corrections is not high enough to decide if the deviation of the Groom- 
bridge 1830 stars from the main sequence of Fig. 6 is significant. Although the 
two tests of the theory previously mentioned seem definite, it should be remem- 
bered that the assumption that the blanketing corrections obtained from the Sun 
are also valid for early F and late G-type stars is an approximation which can be 
eliminated only when ¢(A) values are available for stars of different spectral types. 
The approximation is probably incorrect in detail and may account for the slight 
over-correction of a few stars in Fig. 6. It will be interesting to see if the blanketing 
line for Ko dwarfs in the (U— B, B— V)-diagram lies along the intrinsic (U—B, 
B-—V)-line because such a possibility is suggested by the five subdwarfs later than 
Ko in Table I which have (U—B)~o. Clearly, the present method of cor- 
recting for line blanketing will fail in this case because the A(B— V) correction 
cannot be inferred from a non-existent 5(U — B). 

As Schwarzschild has indicated (1958, Chapter IV), a check on some of these 
results is possible if the masses of subdwarfs were known. Equations (9) and (14) 
are the theoretical mass-luminosity relations on the two opacity assumptions. If 
the previous conclusion is correct that the colour correction is enough to make the 
mild subdwarfs disappear in the (M,.,, log T,)-diagram and therefore that the 
established low abundance of heavy elements does not affect the structure, then 
the subdwarfs should not deviate from the empirical mass-luminosity relation for 
strong line stars. The binary 85 Pegasi is the subdwarf with the best determined 
mass but uncertainties in the parallax, the semi-major axis, and the mass ratio 
(Eggen 1956, Wyller 1956) are great enough to make the decision of deviation 
from the mass-luminosity relation difficult. At present, all that can be said is 
that the observed mass neither confirms nor denies the theoretical expectation. 
As is pointed out in the previous paper (Eggen and Sandage 1959) the mass of 
ADS 16644 in the Groombridge 1830 group is also too uncertain for a decision 
on this critical point but if observations are continued, a check may be possible. 

22 
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3. Application to the nearby stars and galactic clusters —The (U, B, V)-data 
available for nearby stars and for the Hyades, Praesepe and Coma Berenices 
clusters show considerable differences in the U-B=f[(B-—V)] relation for 
F and G-type stars from B—V=+0™-40 to B—-V=+0™ 70. Within each of 
the three clusters, stars assumed to be free from companions show an almost 
negligible scatter in the (U—B, B—V)-diagram but the systematic differences 
between the clusters is appreciable. For a given B—V, the U—B for the Hyades 
is larger than for Praesepe and Coma Berenices. The Praesepe stars show 
§(U — B) z~ o™-020 with respect to the Hyades in the range considered, while Coma 
Berenices stars show 6(U— B)=o0™-045. Although the Hyades-Coma difference 
is large compared to the quoted probable errors in U— B of o™-008 for the photo- 
metric transfers to each cluster (Johnson 1952, Johnson and Knuckles 1955), the 
Praesepe-Hyades difference is only three times the probable errors of these 
transfers and may not be real. Table VI gives the mean values for the (U-—B, 
B—V)-relation in each of the clusters and illustrates the effect. 





Taste VI 
U-—B, B-—V relations for Hyades, Praesepe, and Coma Berenices 

Hyades Praesepe Coma 
B-V U-B U-B U-B 
0°44 —0°005 os —0°035 
0°46 +0008 — 0°005 —0°027 
0°48 0°020 +0°005 —0O°OI7 
050 0°035 0016 —0°006 
o°52 o"050 0°030 +0°007 
0°54 0°067 0°044 0°022 
0°56 0086 0060 0°040 
0°58 0°107 0'079 0°060 
o'60 0°127 o'Io! 0080 
0°62 o"150 0°126 "102 
0°64 Or175 O'lss O'125 
0°66 0°200 0°184 O'155 
0°68 0°225 ‘a 0°185 
0"70 0°248 ‘ 


The (U—B, B—V)-diagram for the nearby stars shows a much larger scatter 
than is present for the individual clusters. The entire spread places some of these 
stars beyond the extremes of the Coma Berenices and Hyades lines respectively. 
The field stars are apparently much less homogeneous in chemical composition 
than the stars in a given cluster. 

Very similar results were presented by Strémgren (1958) and Chalonge 
(1958) at the Rome Conference on Stellar Populations. Strémgren, using narrow- 
band, photoelectric photometry, found a difference between Coma Berenices and 
Hyades stars of the same order as that derived above but he found no difference 
between Praesepe and Hyades stars. (Since Strémgren’s result is independent 
of inter-stellar absorption, the small difference in (U—B) between the Hyades 
and Praesepe stars may indicate a reddening of o™-o15 in B— V for Praesepe.) 

There are two conclusions from these data. (1) Differences in heavy element 
abundance probably exist between stars in different clusters and in the general 
field, and these differences are not related by a simple time sequence. Ill evi- 
dence from composite (M,, B—V)-diagrams (Eggen 1955 c, Sandage 1957) put 
Coma Berenices younger than Hyades and Praesepe by the order of 2 x 10° years, 
yet the heavy element abundance is smaller in Coma Berenices. The conclusion 
is that the heavy element abundance of prestellar matter does not increase 
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uniformly with the time from the formation of the galaxy because, if it did, 
Coma Berenices stars would have higher metal abundances than either Hyades 
or Praesepe. Although an average increase of heavy elements with time may 
indeed take place in the galaxy, evidently there are local irregularities of chemical 
composition in the galactic plane. 

(2) The second conclusion from Table VI is that differences are expected in 
the position of the main sequences in the (M,, B—V)-diagram for the Hyades 
and Coma Berenices due to differences in the A(B—V) values. The method of 
fitting main sequences to find photometric distances for these clusters is therefore 
more complicated than was once believed. Fortunately, the distance to the 
Hyades is fixed from geometrical arguments by the moving cluster method (van 
Bueren 1952, Eggen 1958). Therefore, if the Hyades is used as the standard 
main sequence to which other clusters are fitted, either the B — V of these clusters 
must be corrected to the (U— B, B— V)-system of the Hyades by Table IV, or 
the AM, corrections of Fig. 5 should be applied at constant B—V. Thus, the 
distance moduli derived by both Johnson and by Sandage for several clusters 
probably require small corrections. Eggen’s procedure of using the slope of 
the Praesepe main sequence and the zero-point of field stars with II > 0”-200 to 
define the standard main sequence also has the same objection because the 8( U — B) 
values of the field stars are not zero. The corrections to the separate values of 
Eggen, Johnson and Sandage are not negligible. For example, Table IV shows 
that the correction to the distance modulus of Coma Berenices is AM, =0™-32 if 
8(U — B) =0™-045. 

The problem of fitting galactic cluster main sequences is a milder case of that 
encountered in globular clusters where the extreme excesses of 6(U — B)=0™-20 
are found. The preceding paper (Eggen and Sandage 1959) shows that if we 
neglect this excess for the globular clusters, an error of about 1™-2 results in the 
position of the main sequence. The entire problem of photometric distances to 
clusters where a fit to the main sequence is made in the Fo to Ko range should be 
re-examined from this point of view. Fortunately, the corrections are 
to be smaller than o™-4 for most galactic clusters because the 5(U — B) values are 
known to be smaller than o™-o5 from data already in the literature. As an illus- 
tration, the tentative correction just mentioned changes the distance modulus 
of Coma Berenices from m—M=4™-80 (Eggen 1955 c) to m—-M=4™-48. It 
should be mentioned that no effect is expected for clusters where the fit to the main 
sequence is made to the late B and A-type stars because a changing metal content 
should not affect the colours and no large change in the position of the observed 
main sequence is expected. 

Finally, we enquire into the expected intrinsic width of the main sequence. 
On the present ideas this sequence can stretch from the Hyades stars, with no 
ultra-violet excess, through the Coma Berenices stars with (U— B)=0™-04, to 
the Sun and mild subdwarfs like those in Table I with (U — B)~o™-10 and finally 
to the Groombridge 1830 and globular cluster stars with 6(U —B)=0™2—4 
range of 1™-5 ata given B—V in the (M,, B—V)-diagram. Since 
1830 group members probably are at one edge of this continuum in the (M,, 
B—V)-diagram and have very weak lines, they have usually been excluded from 
consideration as main sequence stars. The stars which are usually called ‘‘main 
sequence dwarfs ’’ are not distinguished by their line weakening and lie ina band ~ 
width of 5(U-B)=o™-10 in the (U-B, B— V)-diagram. If these stars are 
taken to define the main sequence in the (M,, B— V)-diagram, Table IV shows 
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that this sequence should have a total width of AM, ~ o™-60 due to the effect of 
weak Fraunhofer lines on the B—V colour alone. The (M,, B—V)-diagram 
for stars with 7 20"-100 (Eggen 1955 5, Sandage 1958) does show some spread 
in the Fo G5 range but the ever present question of the accuracies in the trigo- 
nometric parallaxes makes the comparison uncertain. 


Mount Wilson and Palomar Observatories, Royal Greenwich Observatory, 
Carnegie Institution of Washington, Herstmonceux Castle, 
California Institute of Technology, Sussex 
Pasadena, California: 
1959 February 27. 
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THE RADIO EMISSION FROM NORMAL GALAXIES 


J. OBSERVATIONS OF M31 AND M33 at 158 Mc/s AND 237 Mc/s 
R. Hanbury Brown and C. Hazard 


(Received 1959 March 9) 


Summary 


A survey of the radio emission from bright normal galaxies has been 
carried out at 158 Mc/s. This paper reports the results obtained on M31 
and M33. The integrated radio magnitude (m,) of these two galaxies has 
been measured and compared with their photographic magnitude (m,); 
it is found that the difference (m,—my,) is almost the same for both. 

The distribution of intensity across M31 shows two major components: 
a disk and a corona. The corona is responsible for 90 per cent of the total 
emission, has a true axial ratio of about 0-6, and has a major axis which is 
about three times greater than that of the visible nebula. The results show 
that M33 also has a corona; the ellipticity, emissivity and dimensions of the 
corona relative to the visible object are closely similar to the values found 
for M31. 

A comparison between M31, M33 and the Galaxy suggests that the 
emissivity of the corona of the Galaxy is significantly higher than that of 
the other two galaxies. 





1. Introduction.—In previous papers (1, 2, 3) we reported some observations 
of galaxies which were carried out at 158 Mc/s with the 218 ft fixed paraboloid 
at Jodrell Bank. The recent completion of the 250 ft steerable paraboloid (4), 
with its wider field of view and greater gain, has made it possible to observe a 
greater number of galaxies and also to improve the previous measurements. 
A survey has therefore been made at 158 Mc/s of the brightest normal* galaxies 
in the northern sky in an attempt to establish the relationship between their 
photographic and radio magnitudes. Subsequently a few galaxies have also 
been observed at 237 Mc/s; this was only done in those cases where the improved 
resolving power was considered to be useful in interpreting the results obtained 
at 158 Mc/s. 

The account of this survey has been divided into two parts. The present 
paper (1) has been restricted to the two brightest spirals M31 and M33; these 
two objects have been studied in greater detail than the other galaxies and require 
individual discussion. A second paper (II) contains the results on the remaining 
galaxies together with a discussion of all the data. 


2. Description of apparatus.—The observations were made at frequencies of 
158 Mc/s and 237 Mc/s, using a receiver based on that described by Ryle and 
Vonberg (6). The parameters of this receiver are substantially the same as 
those of the equipment used in our earlier observations of M31 (1). | 


* In this context a normal galaxy, following Baade and Minkowski (g), is taken to be a galaxy 
which appears to be normal when examined photographically. i : 
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The aerial system used was the 250 ft steerable paraboloid (4). The primary 
feeds at the two frequencies were similar and consisted of a simple dipole and 
reflector connected to the receiver through 85 ft of low-loss coaxial cable. The 
shape of the aerial beam was measured in detail by observing the intense radio 
sources in Cygnus and Cassiopeia; as expected, it was found to be elliptical 
with the major axis parallel to the dipole. The values obtained for the total width 
of the beam between half-power points are given in Table I. The same 
observations were also used to determine the direction of the axis of the beam 
relative to the axis of the paraboloid. At both frequencies the misalignment 
between the two axes was about 10’ and the appropriate corrections have been 
applied to all the observations. 


TABLE I 
The beam-width of the aerial 


Frequency (Mc/s) Beam-width* in E plane Beam-width* in H plane 
158 a°-2 1°°6 
237 1°°5 ee 


* The beam-width is the total width between directions of half-power gain. 


3. Calibration of the intensity scale-—In a previous paper (2) we defined a 
scale of radio magnitude, using NGC 5194/5 as a comparison source, which was 
intended to make the radio and photographic magnitudes of normal spiral galaxies 
numerically equal. On that scale a source with a flux density of S,;, w. m.~* (c/s)~* 
at 158 Mc/s has a radio magnitude m, defined by the equation 


m, = — 53°45 — 2°5 log Sy5¢. (1) 
- For the purpose of the present work it is more convenient to use the discrete 
source 14N5A as a standard of comparison instead of NGC 5194/5. This source 
has a small angular diameter (7) and is circumpolar at Jodrell Bank; the 
equipment was calibrated by observing it at least once during each observing 
period. Preliminary measurements show that it has a flux density of 
75 x 10-6 w. m.~*(c/s)~! at 158 Mc/s which, from equation (1), corresponds to 
a radio magnitude of +6-86. Whitfield (8) finds that it has a spectral index of 
—o-6, 

It is desirable, when observations of sources at different frequencies are to 
be compared, that their radio magnitudes should be defined in such a way that 
they do not vary greatly with frequency; in particular it is convenient that the 
magnitude of the comparison source should be independent of frequency. For 
this reason we have assumed arbitrarily that the radio magnitude of 14N5A is 
independent of frequency and has the value m,= + 6-86 derived at 158 Mc/s. 

To facilitate comparison of our results with those of southern observers we 
have compared the intensity of 14N5A with that of 12N1A in Virgo. The 
12N1A source is the more intense by 2-77 mag. at 158 Mc/s and hence at this 
frequency its radio magnitude is + 4-09. 


4. Observations of M31 (NGC 224) 
4.1. Method of observation.—The observations were made during the late 
summer and autumn of 1958. The region surveyed was bounded by declinations 
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N 34° and N 47° and right ascensions 00" 10™ and o1" 10™, The beam of the 
250 ft paraboloid was scanned over the region at one degree per minute along 
lines of constant declination, the scans being separated by 0-5° in the central 
parts of the survey and 0-75° in the outer parts. These scans were then linked 
together by scanning along lines of constant right ascension at oo" 10™, 
oo" 40™ and o1r"10™, Finally, to make certain of the extent of the nebula, 
the beam was scanned along the major and minor axes. Each scan was repeated 
several times to minimize the effects of interference and receiver drift; all the 
measurements were made at both 158 Mc/s and 237 Mc/s. Some sample scans 
are shown in Fig. 1. 


INTENSITY 


Se ee Se Be ee ee EE SS Se 
34° 36° 3° 40° 42° 44° 46° 48° 
DEGREES (DECLINATION) 


INTENSITY 








Fic. 1.—Sample scans across M31. (a) A scan along the major axis at 158 Me/s. 
(b) A scan along the minor axis at 237 Mc/s. 
Ordinates : intensity in arbitrary units. 
Abscissae : declination in degrees. 


4.2. The radio isophotes of M31.—The records have been converted into 
contours which show the distribution of brightness over the region; the results 
are shown in Figs. 2, 3 and 4. The contours represent lines of equal observed 
brightness and they have been plotted above an arbitrary zero. , 
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The outline of the visible nebula is indicated in Fig. 4 by the broken line 
which, following de Vaucouleurs (9), encloses 95 per cent of the total luminosity. 
It can be seen clearly that the radio emission is much more widely distributed 
than the light, although the shape of the contours is obviously distorted by 
localized sources and irregularities in the general background radiation. 

The localized sources are more prominent on the 237 Mc/s contours, as 
would be expected with the narrower beam. ‘The source at 00%33™, N 39° 
is particularly prominent, and is probably the source originally identified with 
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Fic. 2.—Isophotes of M31 observed at 158 Mc/s. The contours are marked in arbitrary units 
of brightness temperature. Epoch 1950. 
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Fic. 3.—Isophotes of M31 observed at 237 Mc/s. The contours are marked in arbitrary units 
of brightness temperature. Epoch 1950. 
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M31 by Ryle, Smith and Elsmore (10). It is now known to have an angular 
diameter which is considerably less than that of the visible nebulae and this 
fact, together with the discrepancy between the positions of the two objects, 
makes it unlikely that they are physically associated. 
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Fic. 4.—Isophotes of M31 at 158 Mc/s with the background radiation removed. The contours 
are marked in arbitrary units of brightness temperature. The broken line shows the outline of the 
visible nebula. Epoch 1950. 





The small irregularity which appears on the 237 Mc/s contours at about 
one degree from the centre of the nebula may be a discrete source with a radio 
magnitude of about + 9-6; at the distance of M31 this corresponds to a source 
with a total radio emission about 2-5 mag. greater than that of Cassiopeia A. 
The possibility cannot be excluded that it is a discrete source in M31. 

In an attempt to show the radiation from the nebula more clearly the 
background brightness has been removed from the contours observed at 158 Mc/s. 
The background was extrapolated smoothly over the region of the nebula and 
contours of brightness measured above this somewhat arbitrary baseline are 
shown in Fig. 4. At 158 Mc/s the outlines of the nebula are clearly defined 
on the original records and the distorting effects of discrete sources appear to 
be comparatively small; it is therefore likely that the derived contours in Fig. 4 
give a reasonably accurate picture of the radio emission associated bone’ = 
nebula when smoothed with a 2° beam 
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An attempt was also made to remove the background radiation from the 
contours observed at 237 Mc/s. ‘The derived contours proved to be similar 
to those obtained at 158 Mc/s but are not so reliable on account of their lower 
signal-to-noise ratio. 


4.3. The integrated radio magnitude of M31.—The contours at both frequencies 
have been analysed to find the integrated radio magnitude of M31. At 158 Mc/s 
the radio magnitude was found to be m,= + 5-7, which includes all the radiation 
within the 1 per cent contour. This value also includes the localized source at 
00" 33™, N 39°; this source has a radio magnitude of + 9-0 at 237 Mc/s and, if 
it has a spectrum similar to 14N5A, it must contribute about 5-5 per cent of the 
integrated emission from M31 at 158 Mc/s. Subtracting the contribution from 
this source the integrated magnitude of M31 at 158 Mc/s is m,= +5°8. 

At 237 Mc/s it was possible to measure the radiation with reasonable certainty 
only out to the 5 per cent contour. The integrated magnitude out to this level 
is +6-1. Extrapolation to the 1 per cent contour shows that the additional 
radiation is negligible for the present purpose, and the integrated magnitude of 
M31 at 237 Mc/s may therefore be taken as +6-1. A comparison of this result 
with the value of +5-8 at 158 Mc/s suggests that the spectrum of M31 may 
perhaps be steeper than that of 14N5A. The measurements of Seeger, 
Westerhout and Conway (11) show that at 400 Mc/s the integrated magnitude 
of M31 is 0-6 mag. brighter than 14N5A and therefore support this suggestion. 
Assuming the spectral index of 14N5A to be —o-6 (8), the spectral index of the 
integrated emission from M31 averaged over the range 158 to 400 Mc/s is — 1-1. 

It follows from the values given above that the total power emitted by M31 
at 158 Mc/s is 7-8 x 107 w. (c/s)~'ster—". It has been assumed in calculating 
this power that the distance is 630 kpc and that the flux density of the comparison 
source 14N5A is 75 x 10-** w.m.~* (c/s)~! at 158 Mc/s. 

4.4. The general distribution of emission from M31.—An inspection of the 
individual scans over the centre of M31, especially those at 237 Mc/s, suggests 
immediately that there are two superimposed distributions with very different 
widths. The first distribution appears to be two or three times as broad as the 
visible nebula, while the second distribution is much narrower than the first 
and by comparison is highly concentrated towards the centre of the nebula. 
A sample scan along the minor axis at 237 Mc/s is shown in Fig. 1. 

We propose to treat these two distributions as distinct, and to discuss them 
separately. On the basis of the present results alone such a division is necessarily 
rather arbitrary; however, it is supported by a study of the Galaxy. Surveys 
of the non-thermal radiation from the Galaxy have shown that there are at least 
two principal components of the emission: (i) an extended corona which is 
almost spherical and is responsible for the majority of the emission (12); and 
(ii) a second component which is concentrated towards the galactic plane and 
which appears to be roughly coextensive with the visible parts of the Galaxy. 
The second component has been called the disk by Mills (13). We shall therefore 
tentatively identify the two distributions which appear to exist in M31 with 
those found in the Galaxy, and we shall refer to the broad component as the 
corona and to the narrow component as the disk. 


4.5. The disk component of M31.—The disk component has been separated 
ftom the corona by extrapolating the broad distribution over the centre of the 
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nebula as illustrated by the broken line in Fig. 5. The width of the disk 
component obtained in this way is clearly subject to considerable uncertainty ; 
nevertheless we have concluded that the 237 Mc/s records do not show any 
significant broadening of the beam, and it follows that the total width between 
half-power points cannot exceed about 1°. 


OBSERVED CURVE 


INTENSITY 








4 
NORTH FOLLOWING DEGREES SOUTH PRECEDING 


Fic. 5.—Jllustrating the method of separating the disk component from the corona. The solid line 
represents the observed distribution along the major axis of M31 at 237 Me/s. 
Ordinates : intensity. 
Abscissae : degrees measured along major axis. 


Some further evidence on the width of this narrow distribution has been 
provided by Palmer and his colleagues (14) who have surveyed M31 at 160 Mc/s 


using an interferometer with a baseline of 220A. The records show that there 
is no source near the optical centre of M31 with an intensity greater than about 
0°05 of that of 14N5A. It follows that the maximum width between half-power 
points of the disk component must be at least 20’. 

A comparison of the integrated emission from the disk component with that 
from the corona shows that the disk is responsible for about 10 per cent of the 
total emission from the nebula. 


4.6. The corona of M3x 

(i) The apparent distribution of intensity over the corona.—Since the 
contribution of the disk component is only 10 per cent of the total radiation, 
the 158 Mc/s measurements in Fig. 4 are effectively contours of the corona 
smoothed with a 2° beam. The irregularities in the outline of the corona, which 
are particularly marked on the outside contour, are probably not significant 
and may represent minor errors in the separation of the corona from the 
background. 

The general shape of the contours is c¢lliptical and the observed widths of 
the 1 per cent contour are about 10° (major axis) x 7° (minor axis). The position 
angle of the major axis is about 44° and, taking into account the uncertainty in 
the radio isophotes, this is in reasonable agreement with the value of 37°-7 given 
by de Vaucouleurs (9) for the visible nebula. 

The observed results were analysed in an attempt to arrive at a closer 
approximation to the true distribution of the coronal component. The 
radiation due to the disk was subtracted from the observed distribution, and 
corrections were made for the smoothing effects of the aerial beam. Theanalysis 
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showed that the observations can be reproduced closely if the true isophotes of 
the corona are elliptical with gaussian brightness distributions along both axes. 
The total apparent widths of these distributions are 4°-1 (major axis) x 2°-6 (minor 
axis) between half-brightness points, and 10°-6 (major axis) x 6°-7 (minor axis) 
between 1 per cent contours. The apparent axial ratio of the corona is therefore 
0-6, the axial ratio being defined as the ratio of the minor axis to the major axis. 

These results have been compared with those obtained by Seeger, Westerhout 
and Conway (11) who used a 2° beam at 400 Mc/s. The two sets of observations 
are in reasonable agreement. One obvious discrepancy is that the 400 Mc/s 
contours show some curious structure at the extreme north and south of the 
corona. The results presented in the present paper suggest that this structure 
is not associated with M31, but is due to extraneous features of the background 
radiation. 

(ii) The distribution of emissivity in the corona.—We shall discuss the emissivity 
of the corona of M31 on the assumption that the absorption in the nebula at 
158 Mc/s can be neglected, since calculations for our own Galaxy (15) indicate 
that such an assumption is justified. We shall also assume that the corona is an 
oblate spheroid whose principal axes coincide with those of the visible nebula. 
On these assumptions it can be shown (16) that p, the true axial ratio of the 
ellipsoid, is given by 


»_ (b/a)?—sin®i - 
Pe aol (2) 


where i is the inclination of the nebula to the line of sight and b/a is the apparent 
axial ratio of the corona. Taking 7=14°-5 and b/a=o-6 the true axial ratio of 
the corona is 0°57. 

It was stated in the previous section that the apparent variation of brightness 
along the principal axes of the corona may be represented by gaussian 
distributions. It is simple to show that these distributions are exhibited by 
a model of the corona in which the surfaces of equal emissivity are concentric 
ellipsoids of fixed axial ratio and in which the emissivity follows a gaussian decrease 
with distance from the centre. In this model the emissivity e(a,z) at a point 
with coordinates a, z with respect to the centre of the corona is given by 


2 2 
(a, z) =€)exp (- = — ss (3) 
where €, is the emissivity at the centre, a is measured in the plane of the nebula 
and z is measured perpendicular to the plane. 

Assuming the distance of M31 to be 630 kpc and taking the widths of the 
apparent distributions given in the previous section, it can be shown that 
o,=19°3 and o,=11°0, where a and z are measured in kiloparsecs. It follows 
that the total widths of the corona between half-emissivity points are 45 kpc (major 
axis) x 26 kpc (minor axis), and that the total widths between the 1 per cent 
contours are 117 kpc (major axis) x 67 kpc (minor axis). 

The emissivity of the corona was found by integrating the radiation 
from the model and equating the result to the observed total power from 
the nebula. The emissivity at the centre of the corona was found to be 
0°12 x 10° w. (c/s)~! ster! pe~* at 158 Mc/s. This result is in reasonable agreement 
with the value given by Baldwin (17) of 0-24-0°35 x 10° w. (c/s)~'ster—! pe~ at 
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81-5 Mc/s for the mean emissivity of the corona of M31. He assumed in deriving 
this value that the corona is spherical with a diameter of 36 kpc. 

It is interesting to note that, taking the dimensions of the disk component 
suggested in Section 4.5, the emissivity in the disk must be significantly greater 
than that in the corona. 

5. Observations of M33 (NGC 598) 

5.1. The region around M33.—A detailed survey of M33 was carried out at 
158 Mc/s, the region covered being bounded by declinations N 26° and N 36° 
and right ascensions or®15™ and o1®55™. The results of the surveys are 
plotted as isophotes in Fig. 6. The outline of the visible nebula has been 
superimposed as a broken line. 
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Fic. 6.—Isophotes of M33 observed at 158 Mc/s. The contours are marked in arbitrary units 
of brightness temperature. Epoch 1950. 


The most prominent feature of the map is the source at approximately N 33°, 
o1"35™, which has an intensity of about 60 x 10-** w.m.~* (c/s) at 158 Mc/s. 
Some observations made on this source at 237 Mc/s show that it has an angular 
diameter which is considerably less than 1°. Its position is so close to M33 
that it complicates the observations of the radiation from the nebula, and in 
an early interferometer survey made by Ryle, Smith and Elsmore (ro) it was 
itself identified with the nebula. 

The isophotes also show an extended region of low intensity whose centre 
appears to lie somewhat to the south of the nebula. However, the contour 
system is distorted by severe gradients in the background radiation and an 
inspection of the individual records shows that in fact this extended region 
consists of two sources, One of these sources, which has a small angular diameter 
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and an intensity of 13 x 10-** w.m.~* (c/s)~}, lies at about o1"25™, N29°. The 
other source, which has a large angular diameter, is centred at or" 32™, N 30° 15’. 
These latter coordinates agree closely with the optical centre of M33 and we 
therefore propose to identify this extended source with the nebula. 

After correcting for the effects of the confusing sources, it was found that 
at 158 Mc/s the apparent widths of the extended source associated with the 
nebula were 2°-45 in right ascension and 2°-8 in declination between points of 
half-intensity. When allowance is made for the smoothing produced by the 
aerial beam these figures give the corresponding widths of the radio nebula as 
1°*5 (R.A.) x 2°:1 (Dec.). These figures have been checked by making some 
scans through the centre of the visible nebula at 237 Mc/s. At this higher 
frequency the confusion effects were considerably reduced and it was therefore 
possible to measure the widths of the source with more certainty. It was found 
that the apparent widths were 1°-9 (R.A.) x 2°-5 (Dec.) which correspond to 
true widths of 1°-4x2°:1. These figures are in excellent agreement with those 
at 158 Mc/s and will be taken as the true size of the nebula. 

The measurements imply that the source is elliptical with an apparent axial 
ratio of about 0-66; although the ellipticity is small, direct comparison of 
individual records shows that it is real. The observations are not sufficiently 
accurate to give the orientation of the radio source, but show that the position 
angle must lie between 0° and 40°. This range includes the position angle of the 
visible nebula, which is 20°. 

The integrated magnitude of M33 was estimated from the observed peak 
intensity and width of the radio source. It was found to be +7-8 at 158 Mc/s 


and + 7-7 at 237 Mc/s. Assuming the distance to be 630 kpc the total emission 
from the nebula at 158 Mc/s is therefore 1-3 x 107° w. (c/s)~'ster—!. 


5.2. The corona of M33.—The radio source associated with M33 is more 
extensive than the visible nebula and, by analogy with the Galaxy and M31, 
this suggests that in M33 the majority of the radio emission is also generated in 
an extensive corona which envelops the visible nebula. If it is assumed that the 
corona is an oblate spheroid orientated in the same manner as the visible nebula, 
then it follows from equation (2) that its true axial ratio is 0-45. Taking the 
distance of M33 as 630 kpc the total widths of the corona between half-emissivity 
points are 23 kpc (major axis) x 11 kpc (minor axis). 

The emissivity of the corona was estimated by assuming that, as in the case 
of M31, the emissivity follows a gaussian decrease with distance from the centre. 
The value of the emissivity at the centre of the corona was found to be 
015 x 10° w. (c/s)~' ster! pe-* at 158 Mc/s. 

6. Comparison of M31 and M33 

6.1. Comparison of the radio and photographic magnitudes.—The photographic 

and radio magnitudes (at 158 Mc/s) for the two nebulae are given in Table IT. 


Tasie II 
Comparison of the radio and photographic magnitudes of M31 and M33 


Photographic Radio 
Object Type magnitude magnitude 
(18) 
M31 (NGC 224) +453 | 
M33 (NGC 598) +62 +78 
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It can be seen from Table II that, if it is assumed that the distances of the two 
nebulae are the same, then the total radio emission from M31 is greater by two 
magnitudes than the total radio emission from M33. However, the ratio of the 
total radio emission to the total light emission (m,—m,) shows no significant 
difference although the two nebulae differ in type. This subject will be discussed 
more fully in Part II of this paper. 

The measurements of M33 are not sufficiently accurate and cover too small 
a frequency range to decide if there is a significant difference between the spectral 
indices of the two nebulae. 


6.2. Comparison of the corona of M31 with the corona of M33.—The results 
presented here show that in the case of both M31 and M33 the majority of the 
radio emission arises in an extensive corona which envelops the visible nebula. 
The relevant data on the two coronae are given in Table III. The widths 
quoted are the total widths between the 1 per cent contours. 


Taare III 


The coronae of M31 and M33 
Emissivity at centre, 
Nebula Size of Axial ratio w. (c/s)~* ster? pe~* 
corona of corona at 158 Mc/s 


M31 (NGC 224) 10°°6 x 6°-1 0°57 0°12 x 10° 
M33 (NGC 598) 5°°4X2°"4 0°45 o-1g5 x 10° 


The maximum extension of the corona of M31 is about twice the maximum 
extension of M33. However, within the limits of experimental error, the axial 
ratios and the emissivities are equal. 

The visual sizes of both M31 and M33 have been measured by Holmberg (18), 
who obtained a value for M31 of 197’ along the major axis, and 83’ for M33. 
These values may be compared in each case with the size of the major axis of 
the corona; the maximum extent of the corona of M31 is greater than the 
maximum extent of the visible nebula by a factor of 3-2 while the corresponding 
factor for M33 is 3:9. To compare the axial ratio of the corona with the axial 
ratio of the visible nebula it is first necessary to take into account the effects of 
inclination. Examination of several edge-on spirals (19) indicates that, neglecting 
the central bulge, the minor and major axes of Sb and Sc type spirals are in the 
ratio of about 1:12 and therefore they have an axial ratio of about o-08:1. This 
is much smaller than the axial ratio of the corona of M31 and M33 given in 
Table III. However, although the distributions of the radio emission perpen- 
dicular to the planes of M31 and M33 are much broader than the distributions 
of light the ratio of the widths appears to be much the same for each nebula. 
It may, therefore, be concluded that the sizes of the coronae of M31 and M33 are 
approximately proportional to the sizes of the visible nebulae, the maximum 
extent of the coronae being about three times greater than the maximum extent 
of the visible nebulae. 


7. Comparison with the Galaxy.—It has been shown by Mills (13) that about 
go per cent of the total radiation from the Galaxy arises in the corona and about 
10 per cent in the disk component; this result is in agreement with the relative 
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proportions of the two components in M31. He also concluded that the corona 
of the Galaxy is spheroidal with a ratio of minor to major axis of 0-69 and that the 
central emissivity is 4°3 x 10° w. (c/s)~'ster~! pce~* at a frequency of 85 Mc/s. 
In an earlier study of the corona of the Galaxy at 81-5 Mc/s, Baldwin (12) 
showed that the axial ratio was greater than o-5, and he estimated the mean 
emissivity as 1-8 x 10° w. (c/s)~'ster~! pe~* in good agreement with the results 
of Mills. 

The axial ratio of the corona of the Galaxy therefore appears to be 
similar to the axial ratio of the corona of M31. The central emissivity of 
4°3 x 10° w. (c/s)~! ster—* pe~* at 85 Mc/s, which corresponds to an emissivity 
of about 2-3 x 10° w. (c/s)! ster~! pe~* at 158 Mc/s, is greater by a factor of 
about 16 than the corresponding emissivities of M31 and M33. _ As the total 
radio emission from the Galaxy appears to be comparable to the total emission 
from M31, it follows that the corona of the Galaxy is probably smaller than that 
of M31 by a factor of about 2-5 and that it is more nearly equal in size to the 
corona of M33. These conclusions are in reasonable agreement with the results 
of Baldwin (17); he found the emissivity of the corona of the Galaxy to be 
greater by a factor of 6 than that of M31, and his figures also show that the 
corona of the Galaxy is significantly smaller than that of M31. 
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STEADY-STATE COSMOLOGY 
TREATED ACCORDING TO GENERAL RELATIVITY 


W. Davidson 
(Received 1958 November 17) 


Summary 


The interpretation of steady-state cosmology according to general 
relativity, initiated by McCrea (x), is extended in this paper by the 
investigation of external and internal metrics for the static superposition 
of a spherically symmetric mass concentration, or rarefaction, in a 
“* steady-state ’’ medium. 

It is found that a constant total density p, and negative ager Po 
(=—pe) appear to constitute a natural state of relative stability for the 
steady-state medium, and the analysis suggests that this state may be 
realized between the galaxies of a steady-state universe. From a discussion 
of the world-lines of the medium in the neighbourhood of a static disturbance 
it is deduced that such a disturbance has the effect of countering cosmic 
repulsion, or intensifying it, according as it represents a concentration of mass, 
or a rarefaction, relative to the natural state. A mechanism for the creation 
of matter is presented by which the rate of creation is seen to be strongly 
dependent on local gravitational influence, being positive where the medium 
is expanding, and negative (annihilation) where it is contracting. 

Physical limits of size and density are found for a static disturbance 
of the natural state. It is shown that in a disturbed medium the sign of 
the gravitational field depends on the local state of the medium, and in 
a rarefaction the density of inertia, in the Newtonian sense, is found to 
be negative. A theoretical basis for the relative stability of the natural 
state is indicated. 





1. Introduction.—It has been shown by McCrea (x) that the continuous 
creation of matter required by steady-state cosmology, first proposed by Bondi 
and Gold (2), and by Hoyle (3), can be based on the existing field equations of 
general relativity, without introducing, ad hoc, extra terms into these equations, 
In the theory of continuous creation advanced by Hoyle an extra term is added 
to the equations of general relativity, distinct from the energy tensor, so that 
energy is no longer conserved but created out of nothing at a rate determined 
by the extra term. 

McCrea showed that by absorbing the extra term proposed by Hoyle into the 
energy tensor it could be interpreted as a cosmological stress, which in Hoyle’s 
model was taken to be zero. The metric of space-time in the “smoothed out”’ 
steady-state model is established on the basis of Bondi and Gold’s ‘‘ perfect 
cosmological principle’? to be that of the de Sitter universe. Accordingly, 
the application by McCrea of the field equations of general relativity to this 
metric showed that the stress must be uniform, constant, and negative, being 
the negative of the constant energy density in the model. The matter created 
in the model, at the constant rate necessary to maintain the steady state, was 
then interpreted by McCrea as the mass equivalent of the work done by the 
negative stress in the expansion of the model. 
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In this paper we extend the interpretation of steady-state cosmology according 
to general relativity by investigating external and internal metrics for the static 
superposition of a spherically symmetric mass concentration, or rarefaction, on 
a smoothed out steady-state medium. In this way it is found possible to throw 
more light on zero-point stress and energy, to indicate possible factors in the 
mechanism of steady-state cosmology, and to reveal new phenomena which are 
consistent with general relativity, and perhaps with present quantum theory, 
but which have no analogues in classical mechanics. 

2. The smoothed out steady-state medium according to general relativity.—The 
space-time metrics of those expanding or contracting cosmological models which 
are spatially isotropic at every point, and therefore spatially homogeneous, are 
now well known to have the form (4, 5) 
dr? + r? d@* +r? sin? 0 dd? 

; (2.1) 

(1+ hr/4) } 
where R(t) is a disposable function of cosmic time ¢ and k may take the values 
I, 0, or —1. Such a model serves to represent a universe in which a spatial 
isotropy which it possesses in the large becomes extended to the small, by means 
of an ideal smoothing out of its local irregularities arising from mass concentrations. 

The steady-state metric is a special case of (2.1), with R(t) =e! (a, constant) 
and k=o, thus (2) 

dS? = dt® — e*'%(dr? +r? d@* + r* sin? 0 dd?). (2.2) 
(It will be convenient in this paper to let suffix o refer to quantities associated 
with the smoothed out steady-state model of this metric.) It is well known 
that (2.2) is also the metric of the de Sitter universe in the form obtained by 
Lemaitre (6), and by Robertson (7). This was related to the useful earlier form 


(8) 





ast = dt*— R(t) { 


dS? = (1 — F?/a,*) dr? — dl?/(1 — P/a,*) — ? d#* — [* sin? 0 dg? (2.3) 

according to the transformation 
p= —_ tart faglog (1-2 Jay"). (2.4) 

V1-P/a,? 


The field equations of general relativity applied to the general metric (2.1) 
yield for the uniform, and in general time-dependent, proper energy density p 
and proper pressure p the equations (9) 


Kp= 3R?/R® +3k/R? 
~Kp=2R/R+ R*/R?+k/R? |’ (2-5) 


where a dot indicates differentiation with respect to ¢t and « is the Einstein 
gravitational constant. 

In these equations the cosmological constant has been omitted since it is 
superfluous in the steady-state cosmology with which we are concerned in this 
paper. This is because, for the metric (2.2), p and —p are themselves equal 
constants. Thus, following McCrea (1), we put R(t)=e!%, k=o, and find 
that according to general relativity the density and pressure in a smoothed out 
steady-state medium must be taken to be 


(2.6) 


Po= —Po= 


me 
Ka, 














. 
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3. The exterior metric for a static mass concentration or rarefaction in a 
steady-state medium.—lIt is clear that a static exterior metric for a spherically 
symmetric mass concentration in a cosmological medium, which at great 
distance is in the physical state of the smoothed out steady-state universe 
considered in Section 2, will go over to a form equivalent to (2.3) at sufficient 
distance, and go over to this form at all points when the mass of the concentration 
tends to zero. On the other hand, in the neighbourhood of the disturbance 
it will approximate to the well-known Schwarzschild solution. We shall 
therefore adopt the procedure of writing down a well-known metric which satisfies 
these conditions and then interpreting the associated physical state according 
to general relativity. 

It is clear that a metric satisfying the conditions described is 


dS? = (: -*- 5] da ar /(x- =. =) —P.d6*— Pein? dg, (3.1) 


where m is a constant. It will be noted that the metric (3.1) is familiar as that 
for a Schwarzschild particle embedded in the empty but curved world which 
results from admitting the cosmological constant A into the field equations, 
(10). The metrics become formally the same if we put 1/a,” = A/3. 

However, as stated earlier, we abandon the cosmological constant in this 
paper and interpret the curvature of space-time for large / as arising from the 
smoothed out steady-state medium whose density and pressure are given by 
(2.6). A new problem therefore arises of discussing according to general 
relativity how the physical state of the medium is modified by the permanent 
presence of the body embedded in it. The author is aware that the possibility 
of this interpretation of the metric (3.1) has been recognized before (e.g. Pirani 
(11)), but puts forward as new the present investigation of its implications. 

From the equations relating the components of the energy tensor Ty to the 
coefficients of a static metric pf spherical symmetry (12) it is a straightforward 
calculation to show that, for the case (3.1), 

T?= a 5, (3.2) 
where 5” represents the Kronecker symbol. The Greek indices may take the 
values 1, 2, 3, 4 associated with the coordinates /, 0, ¢, + respectively. 

We shall make the reasonable assumption that the stress in the exterior 
fluid can be attributed to the random motion of primary particles, material and 
otherwise, and is therefore everywhere isotropic relative to observers in locally 
inertial frames which are instantaneously moving with the fluid. That is, the 
medium constitutes a perfect fluid characterized by a proper energy density p 
and pressure p. The 7¥ are therefore related to the 4-vector of velocity of the 
fluid, in its mean motion, by the relation 


T= Baal? +p) a (3-3) 
By spherical symmetry we can put d0/dS =d¢/dS =o, so that T2= T= —p. 
Hence, by (3.2), 
p=- 3, ; 
Kay . : (3-4) 
=o, by (2.6) ad 





312 W. Davidson, Steady-state cosmology Vol. 119 


Furthermore, by (3.2) we can substitute 7,4= —p,) into (3.3), so that we find 
Bule+p)(dr/dS)?=0, and since gy,#0, dr/dS#o0, we must have p+p=o. 
Whence 

3 ) 
wy" (3.5) 
=po, by (2.6) J 


It transpires, therefore, that at all points for which the metric (3.1) is valid, 
that is, outside the postulated mass concentration, the density and pressure 
take the constant values associated with the smoothed out steady-state universe 
of metric (2.3). The presence of a stable mass concentration in the medium 
does not evidently affect the maintenance of this density and pressure. Thus, 
according to general relativity the total density py and negative pressure py appear 
as a natural state of relative stability for the steady-state medium. A possible 
mechanism for this stability which is suggested by general relativity will be 
indicated in Section 11. 

Since we have postulated the existence of a stable mass concentration in the 
medium it is necessary to suppose that the relative stability of the medium can 
be upset by sufficiently powerful disturbing factors leading to local gravitational 
condensations. It will be shown in Section 11 that according to general relativity 
there is a critical state of disturbance of the medium at which stable mass 
concentrations may in fact occur. 

If m <o in the metric (3.1) the disturbance of the medium has to be interpreted 
as a rarefaction, and not a concentration, as will be shown in Section 8. Reasons 
will be given in Section 11, however, for believing that such rarefactions cannot 
be permanent. Nevertheless, its inclusion here as a case of stability will lead 
to conclusions which may reasonably be assumed to have some relevance to 
slowly varying disturbances involving rarefactions of the medium. 

If one now considers the actual expanding universe, that is of galaxies immersed 
in an intergalactic medium of an unknown nature, the notion of a steady state in 
the sense just described might be realized by supposing that the intergalactic 
medium should, despite its general expansion, maintain a natural state 
characterized by a constant density p, and constant stress py. If exceptionally 
disturbed, due perhaps to the random motions of galaxies relative to the medium 
as envisaged by Sciama (13), the medium would contract locally to form new 
galaxies. 

The density p,) given by (2.6) may be evaluated approximately by identifying 
, as the reciprocal of the Hubble constant for the model. In this way Hoyle (3), 
and more recently McCrea (14), have arrived at the conclusion that in view of 
the estimated amount of matter in galaxies, the great bulk of the matter in a 
steady-state universe must reside in the intergalactic medium. ‘This fact lends 
support to the above interpretation of our results. That is, the constant density py 
and negative pressure py are not necessarily to be regarded as, in a somewhat 
obscure sense, ‘‘average’’ characteristics derived for a steady-state universe 
by smoothing it out, since it is possible that they may actually be realized in the 
intergalactic medium of such a universe. 

The physical origin of the negative pressure in such a universe can at the 
present time be only a matter of conjecture. McCrea (1) did not propose any 
physical basis for zero-point stress, He pointed out, however, that a zero-point 


p= 
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stress, and perhaps an associated zero-point energy, other than what were usually 
assumed arbitrarily on the basis of classical mechanics, were consistent with the 
field equations of general relativity and also with the virtual zero-point properties 
accorded to space by the current quantum theory of fields. 

In this paper (Section 10) we shall find that under certain conditions our 
steady-state medium possesses inertial characteristics which are similar to those 
that have been described by Massey (15) in relation to negative energy states 
in present quantum theory. In addition, as shown in Section 9, the negative 
pressure in our steady-state medium must be associated with the generally 
repulsive gravitational field that exists in it. Indeed, if the expansion observed 
in the actual universe is in fact accelerating, then general relativity would tell 
us, ruling out the unsatisfactory cosmological constant, that there must exist 
both anti-gravity and negative stress-energy. 

4. The world-lines of the exterior medium.—It is inherent in the representation 
of the undisturbed steady-state medium by the metric (2.2) that the fluid elements 
are treated as at rest in the coordinate system r,0,¢. It follows that the 
medium is in motion in the coordinate system /, #,¢ of (2.3). In fact the 
4-velocity of the medium in this system will be, by (2.3) and (2.4), 


di d6 dp dr\ [il I 
(S-55-95' 95) ~ (re) (4-1) 
We proceed to find the corresponding 4-velocity for the motion of the medium 
concerned in (3.1), noting that we have not now got a reference system 
corresponding to (2.2) in which the medium is at rest. 
The motion is naturally required to be purely radial, i.e. d0/dS=d¢/dS=o, 


Also, since we are explicitly concerned with a local disturbance of the steady 
state, the motion must agree with (4.1) for large J, i.e. 


di 1 


ae a (/ large), (4-2) 


and agree exactly when mo. 

It should now be observed that, so far as the energy tensor is concerned 
(equation (3.3)), the motion may be quite arbitrary, for this tensor does not 
depend on the motion when p+p=o. It should be noted, too, that the usual 
proof in general relativity that the world-lines of the elements of an isotropic 
fluid are geodesics if the pressure is constant fails when p+p=o0. (Thus, in 
general, the equating to zero of the covariant divergence of (3.3) yields, for 
the case when p is a constant, the relation (p + p)A’“., =0, in which *= dx«/dS. 
Unless p+p=o0 one may therefore conclude that the geodesic equation A’A*,, =o 
is satisfied.) Nevertheless, even in the present exceptional case it is reasonable 
to expect that in any local inertial frame the medium will in fact maintain the 
law of inertia for a fluid of constant pressure, and therefore follow a geodesic 
in a general reference frame. A geodesic 4-velocity for the medium will therefore 
be postulated. This postulate is of course satisfied in the undisturbed medium 
(m=0). 

It may easily be shown that the relevant geodesic equations are 


d?] m Td 
a At ap (+3) 


21, ee rey 
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From (4.3) we see that the presence of a mass concentration (m> 0) has the effect 
of countering cosmic repulsion in the neighbouring medium, while a rarefaction 
(m <0) intensifies the repulsion. 

The integrals of these equations are 


MO OS ge 
ds =a aie sis "ae (4.5) 


dr sd 
a8 ~~ am’ (+6) 
fae 
where « is a constant of the motion. The graphs of dl/dS as given by (4.5), 
in the range />0 for the cases m>0, m=0, m<o with varying «, are shown in 
Figs. 1, 2, and 3 respectively. We comment on these curves as follows: 


(i) m>o. Three curves are shown which 
indicate the three different types possible. 
. Curves of type I correspond to 


ke 2S 


for which (di/dS)* has a minimum non-zero 
value. This will occur where d*//dS? vanishes 
which, by (4.3), is at /=ly where 


ly* = ma,?. (4-7) 
The surface /=/, forms a neutral surface 
in the medium where we can regard the 
attraction of the mass m as just balancing 
cosmic repulsion. Curve II is unique, 
corresponding to 


(-3(0)) 


for which di/dS vanishes at the neutral surface 
l=ly. Type III curves correspond to 


{ m\ 2/3) 2 
a< ~3(2) } 4 
FIG.2-m=0 a 


These curves have two branches in between 
which, for a range of / containing /=ly, 
dl/dS is imaginary. 














(ii) m=o. This is the case of the un- 
disturbed medium. From the present point 
of view of possible geodesic motions the curves 
divide themselves into three types similar 
to (i), according as 121. The case a=1 
gives dl/dS = + l/ay. 


(iii) m<o. In this case the acceleration 
NG 3~ see d*l/dS* is always directed away from the 
disturbance. Three typical curves are shown. 
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As a->oo the curves tend to coincide with di/dS= +i/a).. In general di/dS is 
imaginary for sufficiently small /. 

We now consider which of these curves may represent the motion of the 
medium in the various cases. Equation (4.5) indicates that |dl/dS|+l/a, when 
1 becomes large, for any finite «, so that condition (4.2) is satisfied if we take for 
dl/dS the positive root of (4.5) when / is large. The fact that dl/dS-l/a, 
independently of « and m shows that the natural large scale motion of galaxies 
and medium in a steady-state universe can be described by di/dS=I/a,. The 
interest of the present discussion is in the small scale effects of disturbances in 
the intergalactic medium. 

It is naturally required that the motion of the medium should be continuous 
with J. Hence, when m>o we must take the upper branch of curves I and the 
right hand branch (upper half) of curves III (or those portions of them which 
correspond to the exterior of the mass concentration). Both these curves would 
require the medium to be moving outwards everywhere, even near the disturbance. 
This might be possible if m is very small. 

If m is sufficiently large, as in the case of a galaxy, the unique curve II would 
seem to be the most reasonable possibility. In this case the medium would be at 
rest in unstable equilibrium on the neutral surface |=ly, accelerating towards the 
mass concentration for 1<ly, and accelerating outwards for 1> ly. 

In case (iii), when m<o, the medium must everywhere move away from the 
rarefaction. There seems nothing to choose between the various curves for 
this case and indeed they may all be possible since, as we show in Section 11, 
a rarefaction is likely to be unstable. 

5. The rate of creation of matter outside a static disturbance in the medium.—In 
this section we consider the important question of how the rate of creation of 
matter associated with steady-state cosmology is affected by the proximity of 
a mass concentration or rarefaction in the medium. 

For this purpose we divide the proper energy density p of general relativity 
into two parts. One, denoted by p.., will incorporate the contribution from all 
known kinds of ordinary energy, including the energy equivalent of the mass 
of fundamental particles and also their associated kinetic energy due to random 
velocities, and the energy of radiation. The other part, p_, will represent 
‘‘zero-point’’ energy incorporating any other energy which may possibly 
contribute to p in general relativity. Thus 


p=p.+p- (5-1) 
Similarly the pressure p may be written 
p=p.+P-- (5.2) 

The proper density p and pressure p of a perfect fluid are invariants in general 
relativity. We shall assume that the properly measured quantities p,, p_, P,, 
and p_ are also invariants. This means, for instance, that all would 
infer, if it were possible, the same ratio of plus to minus energy in the rest mass 
of a given fluid element. 

According to general relativity total energy is conserved. Creation of ordinary 
energy must decrease the absolute value of zero-point energy, and vice versa 
for annihilation of ordinary energy. The work done by stress must of course 
be taken into account in the conservation process. For a perfect fluid this is 
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expressed by the generally covariant equation of conservation of energy, viz. 
(pA”).,, + pA’, =0, (5.3) 
where A“ is the unit tangent to the world-lines of the fluid elements. 
In virtue of (5.1) and (5.2) we may write (5.3) in the form 


(p4A"), ,+PA?,= —{(pr’),, +P Al}. (5-4) 
The left-hand side is the invariant representing the rate of creation, FE, of 
ordinary energy in unit proper volume. Similarly, the right-hand side is — E_, 
where E_ is the rate of ‘‘creation’’ of zero-point energy or the rate of 
annihilation of ordinary energy. 

For the steady-state medium p and p are constant, equal to py, and py 
respectively. We have found in Section 3 that according to general relativity 
this density and pressure are maintained in the medium, despite the influence 
of local static disturbances on its motion. 

We now make the reasonable assumption that in this steady state the balance 
between what may be termed “‘ positive’’ and ‘‘ negative’’ states of energy will 
depend only on the total density py and pressure py, and consequently p, and p., 
p_ and p_, will also remain constant independent of local influences on the 
motion. 

Accordingly, for E., in the steady state we can write 


E,=(p.+p.),,, (5-5) 
on I 0 ae 
= s+ Ps) eae gr’), (5-6) 
where g is the determinant of the metric coefficients. Substituting from (3.1), 
(4.5), and (4.6), we derive 


, 1d/, dl 
E,=(p.+p eg? as): 
_ FS 2 2/ 
-3(p,+p.){ Md Hate Sa (5.7) 
a V 1+ (a? —1)a,?/P? + 2ma,2/P) 

It is evident that this rate of creation of material and radiation energy is 
influenced by the gravitational effect of the disturbance on the motion of the 
medium. From (5.5) we see that the local sign of E, will be the sign of 2’,, 
since p, +p. will be positive. Now the invariant A", measures the proper rate 
of expansion of the medium. Accordingly, we conclude that where the medium 
is expanding locally matter is being created, and where it is contracting matter is 
being annihilated. 

By letting />oo, or equivalently «1, mo, we get the value of E, for the 
undisturbed steady-state medium, namely (3/a,)(p,+p,). If p, is negligible 
compared with p_,, as is likely to be the case, then this is approximately 3/a, times 
the density of material energy. This agrees with the results of Hoyle (3) and 
McCrea (1), as interpreted by them. 

As an example of the effect of a static disturbance in the medium we suppose 
that m is large and positive as in the case of a galaxy. We have indicated in 
Section 4 that in such a case the critical value for «, namely «= {1 — 3(m/a,)?*}"2, 
is the most likely. Substitution of this value into (5.7) yields, on simplification, 


ly/l—1y?/ 
E,=34(p,+ {ae x 8 
a) P ) \ VI +2ly/l (5 
where /, is the / coordinate of the neutral surface introduced in Section 4. 
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The behaviour of the disturbance factor is easily investigated. We omit the 
details and conclude as follows. Due to the effect of the mass concentration the 
rate of creation in the medium decreases steadily as | decreases, from the undisturbed 
steady-state value at great distance to 1/+/3 of this value at the neutral surface, 
decreases further for 1</y until, at /={2/(1 + +4/5)}ly=0°62ly, annihilation sets 
in for smaller /. 

We conclude, therefore, that, according to our present analysis based on 
general relativity, the creation of matter in a steady-state universe is by no means 
a statistically uniform and random process in space, but on the contrary is highly 
sensitive to local gravitational in e. 

6. The interior solution for a static disturbance of uniform density.—Our main 
concern in obtaining an interior solution of the field equations for a spherically 
symmetric disturbance in the steady-state medium will be with a sphere of fluid 
in which the density is not greatly different from pp, the constant density of the 
exterior medium. It is only in this case that taking explicit account of the 
cosmological background can be expected to lead to significant information 
about the behaviour of the fluid inside such a sphere; and it is with regions of this 
order of density that any theory of the mechanism of steady-state cosmology 
must be concerned. 

For densities p of this order it might be most natural for p to vary continuously 
in the interior, acquiring the value py at the boundary, and also for the 4-velocity 
of the fluid to be continuous through the interior into the steady-state medium. 
Because it is mathematically more tractable, however, yet valuable to illustrate 
physical principles that will be involved in quasi-static regions also, we postulate 
that the proper density in the sphere is uniform and constant. This ensures a 
static interior metric and a fixed boundary for the disturbance. 

Let p be the constant proper density and p the pressure (which will be a 
variable determined from the derived metric) of a perfect fluid occupying the 
region o </</,, where lis the radial coordinate. Outside this region the conditions 
of the steady-state medium obtain. The energy tensor components T,” will be 
given by (3.3), and their relation to a static metric of spherical symmetry of the 
form 


dS? =e dr? — dP —P d® —P sin? 6 dg? (6.1) 
is expressed by the equations referred to in Section 3 (12). From these 
equations, in a manner completely analogous to the derivation of the well-known 
Schwarzschild interior solution (16), we find that the interior metric can be 
written 


dS? =(A— BV 1—Pja*)? dr — dP/(1 — P/a*)— P.d@®—Psin®0d¢*, (6.2) 
where A and B are constants, and a is related to p by 


3 
==, (6.3) 


while the pressure p is the variable 

Bis Ne! 6 

. a(S aoe e+ 
The constants A and B have to be found from the boundary conditions at 

/=1,. These we shall base simply on continuity of pressure, and continuity of 

the metric components g,,, at /=/,. 
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Putting p=p, when /=/,, where by (2.6) py= —3/(xay"), we find 





I I 
A-Bv1-l{/a@ 2A\@ za) 
Continuity of the metric at /=/,, demands that, by (3.1), 


ee 1? 
Mie STs oe a a Se es 
(A-BvV 1-1,7/a*) =1 ar". 


Whence, by (6.5) and (6.6) 
pee pita 2 em ge 
A= +%(1-@/a,2)V 1-2/4 : 6.) 
B= + }(1—3a*/a,°) } 
The signs correspond in the order given, so that there is no loss of generality 
in taking them both positive, in view of (6.2). 

It is to be noted that our interior solution, which has been joined to the 
exterior metric (3.1) by the boundary conditions stated, satisfies the boundary 
conditions established by O’Brien and Synge (17) for discontinuous mass 
distributions in general relativity. For our case of spherical symmetry these 
conditions require that 

8 ap» ze, Fe Salt ~kete 
should be continuous at /=/,. Here Greek letters represent indices 1, 2, 3, 4 
while Latin letters may represent indices 2, 3, 4 only, corresponding to coordinates 
6, ¢, 7 respectively. It is easily verified that our solution satisfies these 
requirements. 

We point out, also, that the fluid in the sphere is at rest in the reference frame. 
This follows from the fact that for the metric (6.1) the energy tensor components 
T;, T, T vanish identically. Hence, by (3.3) we have, for example, 

dl dr 
e’(p +p) ds dS we, 
and since p+p#o in general, it follows that d//dS=o. Thus the fluid is 
unambiguously at rest in this case. There is, consequently, a discontinuity of 
velocity at /=/,, in general, but this is evidently not inconsistent with the 
O’Brien—-Synge boundary conditions. It is to be observed that mo energy is 
gained or lost by the sphere since T,} also vanishes in the exterior medium due to 
the relation py + py=0. 

7. Limits of density and size for a static disturbance in the steady-state medium.— 
To avoid singularities of the interior metric (6.2), which would correspond 
to unacceptable physical singularities, the two expressions 1—//a* and 
A—Bv1-—F/a@ must be of constant sign in the range o</</,. For instance 
we see from (6.4) that if A— BV 1—//a* vanishes for some / then the pressure 
becomes infinite there. We can assume that the physical situation would adjust 
itself to prevent this. We shall restrict our considerations here to the case a? > 0, 
which by (6.3) means p>o. 

For 1—/?/a® to be of constant sign in the range o</</,, we require simply 


Lja<t. (7-1) 
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In the case of A— BV'1—P/a? we note that at /=/, this expression is 
A ba BV 1 —1?/a*= V1 —12/a* >0o, 

using (6.7). Hence we require A — BV 1 — P/a* to be always > 0 in the specified 
range. That is, ifo</<h, 

(1 — aa!) V7 — Ta H(1- 3a4/a,*) VTP] >. (7.2) 
It is easily established that as a* varies between o and oo (p varying between 
and ©) the most stringent condition which ensures the satisfaction of (7.2) is 


l,Ja< V8/9, (7-3) 


this lowest upper limit for /,/a occurring when p>oo. 

Thus, the condition for the absence of singularities varies between the 
inequalities (7.1) and (7.3) according to the value of p. These inequalities may 
be expressed in terms of measurable quantities if we note that the proper radius 5 
of the sphere is given by 





ied Too =asin-(J,/a). (7.4) 


The physical limits vary, Pa between 


8< [3sin [8 
Kp 9 


and (7-5) 


8< 25: 
Kp 2 


according to the value of p. 

This result indicates a physical upper limit to the size of a static body of given 
density, or an upper limit to the density of a body of given size, in the steady-state 
universe. 

It is interesting to calculate on this basis the upper limit of size for an object 
of the mean density of an average galaxy. Taking this density to be 10-** gmcm~* 
we find 5 < 10’ light years approximately. This would seem to be very generous 
according to observational data, but it may place an upper bound to growth with 
age by capture of smaller condensations. 

8. Sign of m as a criterion for a concentration or rarefaction.—The constant m 
in the case of the Schwarzschild solution for a body of spherical symmetry in 
empty space can be identified as the gravitational mass of the body (in relativistic 
units). But here, when the steady-state background is taken into account, this 
matter needs review. 

By (6.6) we find 


(8.1) 


so that 
m2o according as p2 pp. (8.2) 


Thus, the parameter m of a static disturbance is positive for a mass concentration 
relative to the undisturbed steady-state medium, and negative for a mass rarefaction. 

Since m vanishes when p= pp, the form of the equation (4.3) for the geodesic 
motion of the exterior medium indicates that m may be regarded as the excess 
of gravitational mass, taken algebraically, over that which would occupy the 
same coordinate volume in the undisturbed steady-state medium. 





320 W. Davidson, Steady-state cosmology Vol. 119 


9. Dependence of the sign of the gravitational field on the local state of the 
medium.—The fact that the pressure, given by (6.4), in the interior of the 
disturbance is in general non-uniform indicates that a gradient of pressure is 
in general necessary to keep the fluid in the equilibrium that we have postulated. 
Stated otherwise, there is a gravitational ‘‘ force ’’ in the sense of Whittaker (18), 
acting at each point of the medium. The measure of this force, relative to the 
static fluid, is the initial acceleration d?//dS* of a test particle released from rest 
at coordinate / and ‘‘allowed’’ to follow a geodesic. For this one easily finds 

ee BV 1 (9.1) 
dS*  @(A—-BV1—Pja*) 

We observe that, since A—BV1—P/a? is always positive in the interior 
(Section 7), the sign of d?//dS* depends only on the sign of B, and not on/. Thus, 
by (6.7), 





bas 20 according as p¥$ 3p». (9.2) 


It follows that the gravitational field in a static disturbance is repulsive if p< 3po 
and attractive if p>3po. The sign of the field depends, therefore, on the state of the 
medium. 

It may be mentioned that the sign of the gravitational field as defined rather 
differently by Synge (19) depends on exactly the same criterion as found above. 
Based on his definition of gravitational ‘‘force’’ Synge obtains a local form of 
Gauss’ theorem in general relativity. This states that relative to the world lines 
of matter the local gravitational mass in an isotropic medium is p+3p. For our 
case, referring to (6.3) and (6.4), we find that 


pa 6BV 1-P/a (9.3) 
xa*(A — BV1—P/ja®) 
which takes the sign of B, with the same conclusion as before. 
10. Positive and negative inertia.—The gradient of pressure in the static 
interior, referred to in Section 9g, is found from (6.4) to be 
. — 248! (10.1) 
dl xa 1 —Pja(A—-BV1—Pja*y? 
The sign of this expression throughout the medium depends only on the sign 
of AB. The following cases arise: 








(i) p> 3P0 (a* < $ay") 
In this case dp/dl<o, which is a phenomenon familiar to our experience of 
mass concentrations. 


(ii) p=3p9 (a*=}ay") 

The pressure is now constant at the value p = py= — 3/(Kap*). 
(iii) po<p<3Po (449° <a” <a,") 

For this case dp/dl>o so that, since p=p, <0 at /=/,, the pressure is always 
negative its least value being at /=o. 
(iv) p=py (a?=a,?) 

The sphere is now identical with a portion of the undisturbed steady-state 
medium as regards its uniform pressure p, and density py. The argument given 
at the end of Section 6 establishing that the interior fluid is at rest in the 
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coordinate frame now breaks down since py+p)=0; the fluid must in fact be 
expanding in its natural geodesic motion under the repulsive gravitational field 
that exists when p< 3p, (equation (9.2)), there being no pressure gradient to 
counteract it. 
(v) p<po (a*> a?) 

Here dp/dl<o. ‘The pressure is always negative in the interior but decreases 
with / in absolute value. 


It may be surprising that when p decreases through the value p, the pressure 
gradient necessary to maintain equilibrium changes sign from positive to negative, 
despite the fact that the gravitational ‘‘force’’ at each value of coordinate / 
varies continuously with p, without changing sign. In fact, by referring to (9.1), 
we find that the initial acceleration d?//dS* of a geodesic test particle, released 
from rest at coordinate /, varies continuously with p when p decreases through the 


value pp, taking the value 
ind Fg (10.2) 
dS*} 2, Gye 
when p=ppo 


The explanation of this phenomenon can be expressed by saying that the 
density of inertia, in the Newtonian sense, changes sign from positive to negative 
when p decreases through the value po, so that to produce the same effect the motivating 
force of pressure must be reversed. This is seen as follows. 

According to the mechanics of general relativity the acceleration of the fluid 
relative to the rest frame of the test particle mentioned above must be attributed 
to the pressure gradient in that frame. The equations of mechanics for a perfect 
fluid are expressed by equating to zero the covariant divergence of T,” given by 
(3.3). In a local system of coordinates which is instantaneously co-moving 
with the fluid, but geodesic—that is “‘freely falling’’ unaffected by any stress 
gradient existing in the fluid, three of the resulting equations yield the momentum 
equation 

Vp + (p+p)dq/dt =o, (10.3) 
where q is the velocity of the fluid, dq/dt its acceleration, and Vp the gradient 
of pressure. The fourth equation is the energy equation 


dp/dt + (p+p)divq=o, (10.4) 
div q being the spatial divergence of q. Equation (10.4) is the local expression 
of (5.3). 

The momentum equation shows that the proper density of inertia, in the 
Newtonian sense, for a stressed isotropic medium is p+p; this takes account 
of the contribution to momentum of the stress in a moving medium. The energy 
equation takes account of the work done by this stress. 

These results are well known in various forms, but their application to a 
medium of negative p of the same order as p has not generally been envisaged. 
For our interior solution we find, by (6.3) and (6.4), that 


2A 
p+ P= TA_BVi Pa) (19.5) 
which, since A-— BV1—P/a*>o, takes the sign of A. Thus 
p+pzo according as pZ po. 





(10.6) 
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From (10.3) we see how the sign of p+p determines the sign of the stress 
gradient necessary to induce in the medium a given acceleration dq/dt relative 
to a geodesic reference frame. We understand from (10.6), therefore, the change 
in sign of stress gradient when p passes through the value pp, despite only 
infinitesimal change of the repulsive gravitational field indicated by (10.2). 

Thus, like the sign of the gravitational field, the sense of inertial effects in the 
steady-state universe depends on the local state of the medium, according to general 
relativity. 

11. The stability of the steady-state medium.—The analysis in Section 3 
suggested that the state p=p,, =p, may be a natural one of relative stability 
realized in the intergalactic medium of a steady-state universe. In this section 
we put forward, tentatively, a mechanism for this stability based on the mechanics 
of general relativity. 

For this purpose we make use of the results that have been established in 
this paper for static disturbances in the steady-state medium. We make the 
basic inference from the interior solution that, because of its static character, 
there must be a dynamic equilibrium set up between “‘ positive ’’ and “‘ negative ”’ 
states of energy of the primary particles of the fluid. (We recall that we are 
considering only small disturbances of the steady-state medium.) Consequently, 
the derived pressure distribution must be associated with this dynamic 
equilibrium for the different values of p considered. 

On the basis of the consideration that the expansion of any limited region 
of the intergalactic medium will be extremely slow we make the fundamental 
hypothesis that even when the medium is non-static the pressure distribution 
will be adjusted to the prevailing density to secure dynamic equilibrium between 
positive and negative states. In particular, in a spherically symmetric disturbance 
of the medium in which the average density is p at a given local time ¢ the same 
relation between p and p will hold at this moment as was found in the case of 
a static disturbance of the same size and density. It follows according to this 
hypothesis that when a disturbance develops the relative acceleration of the 
fluid elements will be reduced to zero, and the disturbance will simply expand 
uniformly. Also, from this point of view it is the local expansion or contraction 
of the medium that upsets the dynamic equilibrium between positive and 
negative states, causing a net creation or annihilation of matter to restore this 
equilibrium. 

Suppose then that in a limited region of the intergalactic medium a slight 
disturbance of the steady-state conditions arises so that the average value of p 
changes from the value pp. We now notice from (10.4) that, on account of (10.6), 
the effect of the expansion tends to restore the steady-state conditions, i.e. p>p,, 
p-Po. For since divq>o then we shall have 


dp 20 according as p+pSo, 
dt ; (11.1) 


i.e. according as pS py, by (10.6) 


> 


J 
A basis for the relative stability of the steady state is thus apparent. 

Consider now a disturbance of the natural state in a region where the medium 
is contracting. We have seen in Section 5 that inside the neutral surface /=/ly 
associated with a large mass such as a galaxy there will be, when the natural state 
P=Ppo P=pPo prevails, annihilation of matter for /<0-62ly due to the local 
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contraction of the medium (A,”=divq<o) in that region. Suppose now that 
natural state conditions do not obtain in this region, which might happen if the 
galaxy was moving relative to the medium leading to a concentration of mass 
in its wake (Sciama (13)). By (10.6) p+p>o for a mass concentration so that, 
if divq<o, (10.4) indicates that p would tend to increase further. If the 
point is reached when p> 3pp we notice from (9.2), or (9.3), that the interior 
gravitational field becomes attractive, so that even without the influence of the 
passing galaxy the mass concentration is stable and a new galaxy would form. 

On the other hand the gravitational field in a rarefaction is always repulsive, 
and the rarefaction is unstable by the analysis associated with (11.1). 

The question of the stability of smoothed out expanding universes has been 
investigated on a Newtonian basis by Bonnor (20). Employing the methods of 
Newtonian cosmology Bonnor shows that, for small perturbations from the 
smoothed out density, the classical formula of Jeans giving the critical condition 
for gravitational condensations is applicable to the expanding universes in which 
matter is conserved, but not to the steady-state model. Bonnor finds in fact 
that condensations could not form in the steady-state model from small 
perturbations. 

This conclusion by Bonnor is of course in agreement with what we have found 
here to be the basis of the relative stability of the steady-state medium. It does 
not preclude the possibility, however, of condensations forming from large 
disturbances of the medium under the influence of existing condensations in 
the manner we have described. 

Bonnor’s analysis based on Jeans’ formula involves the notion of condensations 
built up by waves of compression arising by statistical variations of density and 
his discussion therefore involves a second-order equation of motion. Our own 
approach is appropriate to the change of a given portion of the fluid as a whole, 
ignoring wave motion on a smaller scale and requiring only the continuity equation 
and the relation between p and p determined by our hypothesis of dynamic 
equilibrium during slow motions. 
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